KONINKLIJKE NEDERLANDSE AKADEMIE VAN WETENSCHAPPEN 


PROCEEDINGS 


S ERE ae 


PHYSICAL SCIENCES 


j 
rs 
“ 
oe 


VOLUME LXIV - No. 1 


NORTH-HOLLAND PUBLISHING COMPANY - AMSTERDAM - 1961 


JOURNAL DE PHYSIQUE 
_ 10, Rue Vauquella eo 


The complete Proceedings consist of three Series, viz. : 


Series A: MATHEMATICAL SCIENCES 
Series B: PuHysicaL SCIENCES 


SERIES C: BroLtoGgicAL AND MErpIcAL SCIENCES 


Articles for these Series cannot be accepted unless 
formally communicated for publication by one of the 
members of the Royal Neth. Academy of Sciences. 


BIOCHEMISTRY 


SILICON TETRACHLORIDE TREATED PAPER IN THE PAPER 
CHROMATOGRAPHY OF PHOSPHATIDES 


IVa. ORIGIN OF THE GRADIENT AND FACTORS INFLUENCING THE 
COURSE OF THE GRADIENT CURVE 


BY 


H. G. BUNGENBERG DE JONG anp J. Tu. HOOGEVEEN 


(Communicated at the meeting of April 30, 1960) 


1. Introduction 


In Parts II and III of this series we have investigated the mechanism 
of the chromatography of egg phosphatides on non-impregnated paper 
and on SiCly-treated paper, using di-isobutylketone—acetic acid —H2O 
(50:25:5) as mobile phase and the large slit-feeding apparatus as 
chromatographic device [1] [2]. 

We came to the conclusion that in both cases the mechanism must 
be considered to consist of adsorption chromatography combined with 
gradient elution. The sequence and degree of resolution of the components 
are different on the two papers, due to the adsorbent being different 
(cellulose; impregnated cellulose). The gradient curve for acetic acid is, 
however, practically the same for non-impregnated paper and for HCl-free 
SiCl,-treated paper. For the aims expressed in the title we may, therefore, 
use both non-impregnated paper and SiCl,-treated paper. We begin 
(section 3) with investigating the influence of the available air-space in 
the chromatographic device on the separation of spots on the course 
of the gradient curve. The results will allow us to draw some general 
conclusions whereby in principle the influence may be predicted of a 
number of other factors. 


2. Methods 


In the present study we use—when not stated otherwise—Schleicher and Schill 
paper 2043b in sheets or strips which are cut in such a way that the direction of 
ascension of the mobile phase coincides with the direction of the smaller suction rate. 

Part of the experiments are performed with non-impregnated paper, part with 
acid SiCly-treated paper of different age and neutralized SiCly-treated paper. In 
preparing the latter we follow the directions given in Part I [3]. 

The phosphatide preparation used throughout this investigation is the ““Reference 
Mixture”, a mixture of column fractions containing lecithin (L); cephalin (C); 
sphingomyelin (8); lysolecithin (LL) and lysocephalin (LC) in proportions suitable 
for chromatography 1). By means of a self-filling capillary pipette (of about 5 mm% 


1) We thank Dr. G. J. M. HOOGHWINKEL for putting a solution of this reference 


mixture dissolved in chloroform-methanol (4:1) to our disposal. 
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capacity) the phosphatide solution (in } % or 1 °% concentration) is brought up 
on the starting points. As chromatographic devices various types of apparatus are 
used especially constructed for our aims, which are described in the experimental 


sections. 
Chromatography was performed in a dark room at 20° (when not stated otherwise) 
with di-isobutyl ketone—acetic acid—Hg2O in the proportions of 50:25:5 or 40:25:5. 
After drying the chromatogram is stained with the Acid Fuchsin—UO2(NOs)e2- 
0.01 N HCl staining solution described in former communications [4] [5]. After 
drying once more, the spots are cireumlined with a pencil. The method used in 
determining the gradient in the mobile phase has been described in detail in Part ITI, 


and will not be repeated here [2]. 


In a few experiments washed SiCl4-treated paper has been used. This kind of paper, 
not yet used in the preceding parts of this series, gives better chromatograms 
(smaller and better separated spots) than acid paper which has hung sufficiently long 
in the air to have lost practically all HCl. 

Sheets of 17 x 28 cm are prepared as follows. We first prepare SiCly-treated 
paper (pulling the paper in 10 see through a 2 % solution of SiCly in CCl4) then 
hang it 30 min in air and subsequently wash each sheet twice during 5 min with 
250ml distilled water. The sheets are then hung in air and can be used three days 
later. As far as we know it does not change its properties with time. 


3. Haxperiments with three chromatographic devices differing widely as to 
available air space 


a) Apparatus 


In a previous communication it was observed that with the large 
slit-feeding apparatus (air space =5.9 1) the spots lie farther apart on the 
chromatogram than with the small slit-feeding apparatus (air space = 0.661), 
though in both cases the front was allowed to ascend to the same 
height [6]. It was concluded that evaporation from the wetted strip may 
play a role in the chromatography. 

From the above observation one should not conclude that increase of 
the air space leads to a larger mutual distance of the spots by virtue of 
the fact that in the large apparatus paper sheets of 17 em width are used 
and in the small apparatus paper strips of only 4.5 cm width. At the same 
front height the evaporating surface of the wetted paper is in the large 
apparatus 3.8 times larger than in the small apparatus. As two factors 
are different (volume of the air space and area of the wetted paper), it 
remains unknown why the mutual distance of the spots on the chromato- 
grams is different in the two kinds of apparatus. 

Four our present aims we shall use three forms of apparatus with very 
different air spaces, in which equally broad strips (4 em) are used. The 
front will be allowed to rise to the same height above the immersion line 
(21 cm) and with starting points on the same distance from the immersion 
line (4.5 cm). These devices will be called A, B and GC, which letters will 
also be used in the figures and tables in this section. A short charac- 
terisation will now follow: 


A = apparatus given in fig. 1. Total air space as used in the experiments 
= 130 ml, of which 60 ml in the upper important part. Details are 
given below (see small print). 

B = the small slit-feeding apparatus. described in a former publication [6]. 
It consists of a glass cylinder of 660 ml capacity, closed by a cork 
which bears the slit-feeding tube. 


C = Erlemeyer flask of 4000 ml capacity, closed by cork bearing the 
slit-feeding tube of a small slit-feeding apparatus. 


Apparatus A, consisting of three loose parts I, I and III, is placed with I in a 
wooden block a attached to an iron plate b (base of a support stand), to give the 
whole assemblage a firm stand. 

Part I is a rectangular glass tray with 180 ml volume. 

Part II consists of a plate of plate glass, provided with a slit (47 mm long and 
0.8 mm wide) to which is attached the “house” of plate glass; inner dimensions 
of crosssection = 5 mm xX 47 mm). 

The total air space in this ‘‘house”’ from slit to upper end is only 60 ml. Part I, 
built up of plates and strips and blocks of glass is sealed together with “Devcon, 
wire resistant”’, (black in the figure) while care is taken that as little as possible 
of this cement stands in direct contact with the inner air space of the house 1). 


Compare the crosssection Ila. 


1) In our first attemps to build an apparatus of type A, Part II and ILI were 
made wholy of Perspex. This material proved to be quite inadequate for our pur- 
poses, the inner surface of the “house” swelling and liquefying in the vapours of the 
mobile phase. As a result it was decided that the apparatus should be built of glass, 
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Part II] is a glass piece, fitting loosely in the upper opening of Part Il. It has a 
slit and a hole. The upper end of the paper strip IV with perforation hole is put in 
the slit so that the holes coincide, after which a short glass pin is put in the holes. 
Before fastening the strip, it has been provided with phosphatide spots on the 
starting points on line st and has been pushed with its upper end through the slit 
of the groundplate of Part Il, until it comes out above the opening of the house. 
After fastening to Part II, the latter is placed into the opening of Part I. 

To prevent that the paper strip makes contact with the walls of the glasshouse 
it must be provided with a weight V at its lower end. For this purpose the lower 
end of the strip has four incisions ¢ as shown in the figure. By folding the paper as 
shown to the right in fig. 1, we can now push in the glass-tube V with sealed-in piece 
of iron rod (volume 10 ml; weight 53 g). After having controlled that the strip hangs 
stretched in the glasshouse of II, nowhere touching the walls, Part II is placed as 
a lid on Part I. The tray has been filled beforehand with so much of the mobile 
phase (100 ml), that after immersion of the weighted lower part of the strip, the 
meniscus coincides with the pencil line 7 on the paper. After the front has risen to 
the pencil line fr, Part II is lifted up and with a pair of scissors we cut off the lower 
part, lift up Part III with which the paperstrip is removed out of the glasshouse 
of Part II. 


b) Particulars on the experiments 


With the above apparatus A, B and C chromatograms have been made 
on non-impregnated paper, on practically HCl-free SiCl4-treated paper 
and on neutralised paper, using as mobile phase di-isobutylketone — acetic 
acid — H,O0 =40:25:5 and allowing the front to rise 21 em above the 
immersion line. With the two first-named papers and the same mobile 
phase chromatograms were made without applying the Reference mixture. 
They served for the determination of the acetic acid gradient in the mobile 
phase wetting the paper. On the strips pencil lines were drawn beforehand 
at 3; 6; 9; 12; 15; 18 and 2] cm above the immersion line. The front 
was allowed to ascend to the 21 cm line. Then the strip was rapidly cut 
with a pair of scissors along the pencil lines. We obtained seven zones 
in which we determined the mean acetic acid concentration of the mobile 
phase wetting the zone. The method is described in details in Part TIT 
of this series. 


c) Results with non-impregnated paper 
0) GRADIENT DETERMINATIONS 


Two independent series of gradient determinations have been made, 
to get an idea of the possible errors. Compare Table I. As the analysed 
paper zones (3 <4 cm) are much smaller here than in the gradient deter- 
minations in Part IIT (4.517) one would beforehand expect greater 
errors than in Part IIT (where an error of at least 1 weight percent 
acetic acid was present). ' 

When we compare the results of the series a and b in the Table this 
is confirmed. We find here six times a difference of about 2 weight %, 
once of 3.5 % (C IIT) and once even of 6.1 % (BIV) 


o 


When for each series apart the data of Table I are plotted two figure 
are obtained which resemble one another much. Thus despite the errors 
the conclusion is not endangered that the descend of the gradient curve 


is steepest in apparatus A, less steep in apparatus B and still less steep 
in apparatus C. 


Ds) 


TABLE I 
Gradient determinations in non-impregnated paper 
ee ee eee eee 


Mean acetic acid concentration (weight %) and in parentheses 
| location of centre of zones (in % of fri) in two independent series 
Zone | of experiments a and b 
Apparatus A Apparatus B Apparatus C 
a b- 5 ‘T b a b 
a Oo 

I 43.4 43.6 44.2 45.2 41.5 40.7 
(7.0) | (7.2) (7.1) (7.0) (7.3) (6.9) 
Mt | 42.5 40.5 44.3 44.7 34.2 36.2 
(21.3) (2125) (21.6) (21.0) (21.9) (ZR) 
INE 40.1 7 41.0 43.1 30.8 34.6 
(35.8) (3525) (36.2) (35.4) (36.4) (35.6) 
IV 42D 4 31.9 38.0 30.4 30.7 
| (50.2) (49.5) (50.4) (50.0) (50.5) (49.9) 
AYE oe 26.8 24.2 26.3 26.3 28.3 
(64.6) (64.0) (64.7) (64.3) (64.6) (64.2) 
VI 62 ie 8a 19.0 18.4 24.5 24.4 
(79.0) (78.5) (79.0) (78.6) (78.6) (78.5) 
Vil Pero" apie. 0:8 Wen 18.4 24.1 24.7 
(93.1) (92.9) (93.1) | (92.9) (92.8) (92.8) 


* extract unreliable by leak of the flask. 


f CHROMATOGRAMS 


The chromatograms obtained from the “Reference Mixture” (5 mm? 
of a 4 or a 1 % solution have been applied on the starting points) are 
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given in fig. 2. The Ry,-values 1) of the spots are given in the next survey 
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Apparatus | Lee | S Ties Le 
A 0.65 0.65 0.58 
B 0.695 0.645 0.48 
C 0.728 0.645 0.365 


1) The notion Ry, = distance base of the spot to immersion line (7) divided by 
Q . . . . 

fr; distance of front to immersion line, as a characteristic number for describing 

the position of a spot on a chromatogram, has been introduced in former publications 


[7] [1]. 
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y) COMPARISON OF CHROMATOGRAMS AND GRADIENT CURVES 

We read from Table I that in apparatus A (with the smallest air space), 
the gradient curve proceeds first nearly horizontally, and then descends 
very steeply. The corresponding chromatogram A in fig. 2 shows two 
spots only, which lie close together. The Ry, values of these spots lie in 
the tract of % fr; values where the steep descent of the gradient curve A 
is situated. 

The gradient curve in apparatus B proceeds downwards much less 
steeply (see Table I) and the spreading of the spots on the chromato- 
gram B in fig. 2 is greater now, besides we obtain already three spots. 

With apparatus C, the slope of the gradient curve is still smaller 
(see Table I) and the spread of the three spots is still much greater 
than with apparatus B. Compare chromatogram C in fig. 2. Thus there 
is a very distinct correlation between spread of the spots and the slope 
of the gradient curve. 

A further discussion of the results will follow in section 4. 


d) Results with practically HCl-free SiCl4-treated paper 
a) gradient determinations 

SiCly-treated paper was used after hanging 91 hrs in air. In this 
time span most of the HCI present in the paper directly after preparation 
evaporates. The very small quantity still present can be completely 


TABLE II 
Gradient determination in practically HCl-free SiCly-treated paper 
| Mean acetic acid concentration (weight %) and in parentheses 
Zone | location of centre of zones (in % of fr;) id seal 
Apparatus A 1 Apparatus B Apparatus C 

Nee al 41.7 43.5 37.2 
(7.2) (Tak) (7.2) 

It | 41.3 42.8 32.5 
(21.6) (21.7) (21.6) 

iM | 41.9 40.7 Bie 
(36.1) (36.2) (36.1) 

TAY! 41.2 35.4 29.1 
(50.5) (50.2) (50.2) 

Vv 26.4 25.1 27.8 
(64.8) (64.2) (63.7) 

VI 8.4 Hg) 24.0 
(DEN) (78.5) (78.0) 

VIL 2.1 18.1 23.6 
(93.1) (92.9) (92.7) 


8 


neglected, compared with the acetic acid present '). Thus the titrations 
of the extracts of the zones refer practically only to acetic acid. 
The results of the gradient determinations are given in Table IT. 


B) Chromatograms 
The chromatograms obtained from the Reference mixture are given 
in fig. 4. The Ry,-values of the spots are given in the following survey: 


Apparatus | C | L | LC +5 LL 
| | 
A | 0.66 0.62 0.605 0.58 
B 0.67 0.57 0.54 0.498 
Cc | 0.66 0.47 0.40 0.325 


y) Comparison of Chromatograms and Gradient curves 

We read from Table II that the apparatus with smallest air space (A) 
gives a gradient curve with very steep descending branch. The slope of 
the descent of the gradient curve is less in apparatus B and still smaller 
in apparatus C (with the largest air-space). 

The chromatograms in fig. 4, show once more a similar correlation 
between slope of the gradient curve and separation of spots as was present 
with non-impregnated paper: the separation of spots is smallest when 
the slope of descending branch of the gradient is steepest. The separation 
increases as the slope of the gradient curve diminishes. 

These results confirm the picture arrived to in Parts Il and III of 
this series of the mechanism of the chromatography on SiCl4-treated 
paper, namely adsorption chromatography combined with gradient elution. 

Further particulars are discussed in section 4. 


e) Results with neutralized SiCly-treated paper 


We found already in Part III, that in the determination of the acetic 
acid gradient difficulties arise by the NH4Cl present in neutralized paper 
[2]. For this reason we have not made gradient determinations here. The 
chromatograms obtained are given in fig. 3. We see here in principle the 
same influence of the increase of the air-space in the chromatographic 
apparatus as in the case of non-impregnated paper (fig. 2) and SiCl,- 
treated paper (fig. 4). With increase of the air-space the spread of the 
spots increases. So this is an indication that for neutralized papers too 
will hold that with increasing air-space the slope of the gradient curve 
diminishes. 


Next survey gives the Ry, value of the spots. 


Apparatus C L | LC +5 | aa 
| | 
| 
A 0.65 0.59 0.56 0.49 
B 0.65 0.545 0.51 0.46 
Cc 0.64 0.45 0.39 | 0.295 


1) For decrease of HCl content with time see fig. 6 in Part ILI of this series [2]. 


i) 


4. Origin of gradient and role of air s pace in the chromatographic apparatus 


a) Origin of gradient 


We will start our discussion from the combined results of the three 
sets of gradient determinations given in the Tables I and II 1). In Table III 
we give for each apparatus A, B and C the mean acetic acid concentration 
in the seven zones and the mean location of the centre of these zones 
expressed in % of fr; (=distance front to immersion line). The data of 
the table have been plotted in fig. 5. 

In this figure, three striking points draw our attention: 

1. In the apparatus with the smallest air space (A), the gradient curve 
has a very steep descent, in the apparatus with larger (B) and very 
large air space (C) the descent is less and still less steep. 

2. The gradient curves A, B and C intersect one another. 

Taking into account the errors still present (though diminished by 

taking the mean of three sets of determination) it seems possible that 
the intersections points of the three gradient curves coincide. 

Point 3, which suggests the existence of a turning point on the gradient 

curve will not be considered here in more detail. 

Point 1 is important for the question of the origin of the gradient. 

One could for instance suppose that when no air space was present in 
the chromatographic apparatus, the mobile phase would have up to the 
front everywhere the same composition as the mobile phase which is 


sf 


TABLE III 
Gradient determinations in Apparatus A, B and C (mean of the three series 
in Tables I and II) 


Mean acetic acid concentration (weight %) in zones I-VII and 


| 


Apparatus in parentheses location of the centres of zones (in % of fr;)* 
Ss feats Vey yoo Vere vis 

IN 42.9 | 414 | 41.2 42.0 26.6 7.8 12 
(7.1) | (21.5) | (35.8) | (50.1) | (64.4) | (78.9) | (93.0) 


B 44.3 43.9 | 41.6 Bole) 25.2 19.1 17.9 
(7.1) | (21.4) | (85.9) | (50.2) | (64.4) | (78.7) | (93.0) 


C 39.8 34.3 32.2 30.1 27.5 24.3 24.1 
(7.1) | (21.5) | (36.0) | (50.2) | (64.2) | (78.4) | (92.8) 


* fr; = distance front to immersion line. 


1) In Part III of this series, it appeared that the gradient curves on non-impreg- 
nated paper and on HCl-free SiCl4-treated paper run in the same apparatus did not 
differ more than the experimental error (1 or 2 weight % acetic acid). In the Tables 
I and II of the present communication, the experimental errors are greater (2, 
sometimes more than 2 weight %, see Table I), but the same differences also are 
found between the results in Table II and each of the series in Table I. To minimize 
the errors and to obtain, therefore, a better picture of the gradient curves in apparatus 
A, B and C, it is indicated to start from the means of the results in table I and II. 
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sucked up at the immersion line. When an air space is present a gradient 
could then arise by preferential evaporation of acetic acid from the 
upper part of the wetted strip, the lower end always being fed by the 
unaltered mobile phase. 

When this supposition applied one would expect that this preferential 
loss of acetic acid would be greater as the available volume of the 
container is greater. This then would lead to an ever steeper gradient, 
so that at last the acetic acid concentration of the mobile phase in the 
upper part would become smaller. 

When we consider fig. 5, we see that just the reverse actually occurs. 
With the largest air space the steepness of the acetic acid gradient curve 
is smallest, it increases in steepness as the available air space is smaller. 
There is every reason to suppose, that when an apparatus could be realised 
with air space zero, the descent of the gradient curve will be still steeper 
than is shown in apparatus A }), 

In this case, the paper will behave as an adsorption column, which 
preferentially adsorbs acetic acid from the mobile phase. 


50 acetic acid 


40 


30 


20 


10 


20 40 60 80 100 


Fig. 5. Gradient curves obtained with apparatus A, B and C (mean of the three 
sets of gradient determinations in section 3). 


1) Attempts have been made to realise chromatography in paper immersed in 
mercury, but because of undesired wetting effects of the paper by mercury (especially 
of SiCli-treated paper) and other complications, we have not succeeded in constructing 
an apparatus in which the air space approximates zero. 


1] 
b) Role of the air space in the chromatographic apparatus 


The partial vapour pressures of the components of the mobile phase 
will in the case considered above be very different at both sides of the 
small zone in which the very steep gradient is present. If there is no 
air space, no possibility for interchange will exist. As soon as an air space 
is present such an interchange becomes possible by diffusion through the 
air space. Thus acetic acid evaporating from the lower part of the paper 
may diffuse to the upper part of the paper and may be taken up by the 
mobile phase present there. Vice versa, di-isobutyl ketone evaporates 
from the upper part of the paper, diffuses (through the air space) in 
downward direction and is taken up by the liquid wetting the lower part 
of the paper. As a result the gradient curve for acetic acid will diminish 
in slope, whereby the new gradient curve will intersect the very steep 
original gradient curve (see point 2 in paragraph a). 

We may further expect that more of the components can be transported 
through the air space as its volume increases not so much as a direct 
consequence of the enlarged volume itself, but because in general the 
area of the cross section of the air space perpendicular to the paper will 
be greater in a larger vessel. Thus with a larger volume of the container 
a gradient curve will be obtained with smaller slope. This is just what 
we found experimentally, compare in fig. 5 the curves A, B and C. 


c) Resolution of spots on the chromatogram and air space 
in the chromatographic apparatus . 


In the preceding Part I and II of this series we came to the conclusion 
that the chromatography of phosphatides with di-isobutyl ketone 
—acetic acid—H2O mixtures isto be considered as adsorption chromato- 
graphy combined with gradient elution. For each spot a characteristic 
composition of the mobile phase would exist at which that spot is at its 
base eluted with maximal rate (compared to the rate of ascension of the 
front). 

One must thus expect that when the slope of the gradient curve is 
diminished, the relative distance of the spots on the chromatograms will 
increase. This is just what has been found in section 3. 

In fig. 6 we have plotted the R,,-values given in the surveys of section 3 
for the three apparatuses A, B and C. We see that on the three kinds 
of papers used the mutual distance of the spots increases from A over 
IB. to: Ce 

We first discuss the chromatograms on non-impregnated paper. In A 
the upper spot has not yet been dissolved in an 5 and L+C spot, this 
takes place in B and C. 

An interesting point is that with increase of the volume the L+C spot 
moves in upward direction, the S spot remains on its place and the 


LL+LC spot moves downwards. 


100 


90 not impregnated Si Cl, treated SiCl, treated 
(neutralized) 


Fig. 6. Ry, values of the chromatograms of fig.s 2-4. 


When we consider the gradient curves A, B and C in fig. 5, this becomes 
clear at once. The intersection of the three curves lies about at the same 
fr; value namely at about 64 %. The Ry, value of the S spot has about 
the same value, namely 0.65 (A); 0.64° (B) and 0.64° (C). Hence the 
decrease of the slope of the gradient curve from apparatus A over B to C 
should not alter the position of the S spot on non-impregnated paper. 

The L+C spot though not visibly separated must lie in A somewhat 
higher than fr;=0.64 and therefore will shift in upward direction with 
decrease of the slope of the gradient curve (A — B — C). For the LL+ LC 
spot applies the reverse. In A its Ry-value lies lower than 0.64; with 
decrease of the slope of the gradient curve (A — B ->C) it should only 
move downward on the chromatogram. 

Now turning to the SiCly-treated papers we find that here no spot 
moves upwards at increase of the volume of the chromatographic apparatus 
(see fig. 6), The upper spot, which here is C (cephalin) remains, however, 
in its place, the three other spots move downwards with decreasing slope 
of the gradient curve (A + B — C). 

As C lies on both kinds of SiCly-treated paper at about the same Ry, 
value as § on non-impregnated paper, C must lie on or very near to the 
intersection point of the three gradient curves. Hence it will not move. 
The three lower spots of course must move downwards with increasing 
volume of the chromatographic apparatus. 


d) Concluding remarks 


The air space in a chromatographic apparatus is not an inert entity, 
but takes an essential part in the resolution of spots. This knowledge 
allows us to foresee, or at least to explain, the influence of a number of 
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other factors on chromatography, which will be investigated in the 
sections 6-10. 


ae oe 
5. The mechanism of elution and the sequence of spots on the chromatograms 


Already in paragraph a of section 4 we came to the conclusion that 
if no air space were present, the paper would behave as an adsorption 
column, which preferentially adsorbs acetic acid from the mobile phase. 
This sheds light to the mechanism of the resolution of the spots on 
the chromatograms. The elution process may be considered as an 
adsorption displacement by acetic acid from the mobile phase. The 
concentration of acetic acid in the mobile phase required for eluting 
a phosphatide spot at its base will be higher as the adsorption intensity 
of the phosphatide spot will be higher. Thus the lower a spot is situated 
on the chromatogram the stronger it is adsorbed. So the sequences for 
increasing adsorption intensity are: 


L+C<S<LbL+LC (non-impregnated papers) 
C<L<LC+S8<LL (SiCl,-treated papers) 


Taking into consideration that the phosphatides are adsorbed from a 
hydrophobic medium (di-isobutylketone) on a polar adsorbent (cellulose 
or silica-covered cellulose), they will be actively adsorbed by means of 
hydrophilic groups in their molecules. Small differences in lengths of the 
hydrocarbon chains should be relatively less important. So we have to 
consider the ampho-ion part of the molecule (cholinephosphate; colamine- 
phosphate) and eventually present OH groups. 

Non-impregnated paper—cellulose being a neutral substance —makes 
evidently practically no distinction between the ionized colamine group 
and the quaternary ammonium group of choline. This follows from the 
fact that L and C form a composite spot, equally so LL+LC. That the 
last-named spots lie so much lower on the chromatogram has two causes; 
first, they contain an OH group, and further one long hydrocarbon chain 
is no longer present. In sphingomyelin an OH group is present too, but 
it still possesses two long hydrocarbon chains. Thus compared to lecithin 
it should be adsorbed stronger, but not as strong as lysolecithin. 

The above leads to the sequence of the three spots actually found in 
the chromatograms on non-impregnated paper: 


(upper) (L+C)—S—(LL+ LC) (lower) 


On SiCly-treated paper the adsorbent is silica, which envelopes the 
cellulose fibres [2]. Silica being an acidic substance it does discriminate 
between the ionized primary amine group (of colamine) and quarternary 
ammoniumgroup (of choline). The affinity of the latter ionized group to 
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the adsorbent is evidently greater than the former). The outcome is 
a resolution of the composite spots which we found on non-impregnated 
paper. Thus the L spot lies lower than the C spot, similarly the LL spot 
lower than the LC spot. For the rest the influence of the presence of an 
OH-group (sphingomyelin), and the influence of introduction of an 
OH-group accompanied by loss of a long hydrocarbon chain (lyso products) 
remains the same. 

The above explains the actually found sequence of the spots on 
SiCly-treated paper: 


(upper) C—L—(LC +8) —LL (lower) 


The S and LC spots come so near to one another that they can only 
with difficulty be resolved; as a rule they form a composite spot. 


6. Comparison of ascending and descending chromatography by means of 
a slit-feeding apparatus 


a. Apparatus 

The apparatus used is given in fig. 7. It consists of a glass-eylinder I (8.3 cm 
inner diameter and 29 em high, volume about 1.6 1), and a plate of glass II, with 
two holes, in which with corks are fastened a glass rod with hook (a) and the filling 
tube (b) which gives entrance to the slit-feeding tube (c), having a slit 1.5 mm 
wide and 70 mm long. For ascending chromatography the filling tube (b) is slid 
downwards, and the paper strip is fastened to the hook of rod a by means of a 
perforation hole. The coverplate II is shifted to a position such that a coincides 
with the axis of cylinder I. The lower end of the strip dips 5 mm into the slit of 
the slit-tube. 

For descending chromatography rod a is not used. The filling tube b is drawn 
up into the position as shown in the figure. At right angles to the middle of the 
slit-feeding tube a glass rod d is sealed, ending in a hook. The paper strip hangs 
on this hook by means of a perforation hole. It is stretched down by a glass rod e 
fastened in two slits in the free hanging end of the strip. The other end dips 5 mm 
into the slit-tube. The starting points lie 5 em from the immersion line, that is 
1.5 cm past the perforation hole. To make the ascending chromatography comparable 
here too the strip had a perforation hole, 1.5 em below the starting line, while the 
distance of the starting line to the slit was also 5 em. 


b) Results 
a) Non-impregnated paper 

Fig. 8 gives two chromatograms on non-impregnated Schleicher and 
Schill 2043b of $° spots (5 mm3) of the Reference mixture. As mobile 
phase was used di-isobutyl ketone —acetic acid—H20=50:25:5 and the 
front was allowed to rise (or descend) to 22 cm from the immersion line. 
The time needed to reach this was in descending chromatography shorter 
(216 min) than in escending chromatography (256 min). It is observed 


1) In former investigations it has already appeared that the affinity of quater- 
nary ammoniumions to SiQs is greater than of the cations of primary amines. 
Compare in Colloid Science I, Elsevier, Amsterdam (1949) chapter LX, p. 304, fig. 29. 


Fig. 7. Slit-feeding apparatus used to compare descending and ascending chromato- 
graphy. For ascending chromatography the slit-feeding unit is lowered. 


Figs. 8-9. Comparison of ascending (ase.) and descending (desc.) chromatography 
of the Reference mixture on non-impregnated paper (fig. 8) and on practically 
HCl-free SiCly-treated paper (fig. 9) respectively. 
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that with descending chromatography the three spots stand wider apart 
than with ascending chromatography. It must be concluded that with 
descending chromatography the gradient is less steep than with ascending 


chromatography. 


B) HCl-free SiCly-treated paper 

Here too Schleicher and Schiill 2043b is used, but now 12 days after 
preparation with 2% SiCly in CCl. 

Fig. 9 gives the two chromatograms obtained. The results are similar 
to those obtained with non-impregnated paper: The four spots stand 
farther apart with descending chromatography, so the gradient is less 
steep here than with ascending chromatography. 

The time necessary to reach the front line (23 cm from the immersion 
line) is shorter in descending chromatography (222 min.) than in ascending 
chromatography (282 min). 


c. Discussion 


In ascending chromatography the capillary forces move the front in 
a direction against gravity, in descending chromatography in the same 
direction as gravity. So there is no problem why the time for the front to 
reach a same distance from the immersion line is shorter in descending 
chromatography than in ascending chromatography. 

If diffusion were the only means of transport through the vapour 
space, one would from the shorter time available expect that in descending 
chromatography the gradient curve is steeper than in ascending chromato- 
graphy, that is the spots would stand nearer to one another. 

As just the reverse is the case, there must be a second transport 
mechanism which in the case of descending chromatography is much 
more effective in lessening the slope of the gradient curve than in the 
case of ascending chromatography. 

The wetted part of the paper strip will be covered on both sides by a 
layer of saturated vapour of the mobile phase. This layer is denser than 
the vapour in the container, and will move downward along the strip. 
Given this downward movement, there must be an upward movement in 
another place of the container. So a rotation in the vapour space at-both 
sides of the strip (or sheet in the case of a large slit-feeding apparatus) 
sets in, moving downwards along the strip and upwards along the glass 
walls of the container. 

The direction of rotation is the same both in ascending chromatography 
as in descending chromatography, but the efficiency of transport of a 
given component of the mobile phase is different. Let us choose for our 
argument acetic acid. In ascending chromatography the rotation current 
in its downward movement takes up acetic acid mainly from the slit-side 
part of the wetted strip. This acetic acid taken up can be given off to 
the mobile phase at the front side of the wetted strip only after having 
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completed a half cycle in upward direction along the glass wall. Mean- 
while it may have lost part of its acetic acid by diffusion. Thus the transport 
of acetic acid is not so efficient. In descending chromatography the reverse 
applies. Here the rotation current descending along the wetted strip 
takes up acetic acid from the slit side part of the strip and gives it off 
at the front side part of the strip while further descending. Here the 
transport of acetic acid is much more efficient. For the other main 
component of the mobile phase (di-isobutyl ketone) just the reverse holds, 
so that we obtain: 


— ee 


Transport of named component | ascending descending 
towards part of the strip | chromatography chromatography 
acetic acid, towards frontside | less efficient more efficient 
| | 
di-isobutyl ketone, towards slitside | more efficient less efficient 


From the greater spread of the spots obtained with descending chromato- 
graphy than with ascending chromatography we must conclude that 
the efficiency of the transport of acetic acid is a much more important 
factor in lessening the slope of the gradient curve than the efficiency of 
the transport of di-isobutylketone. We will not try to explain this in 
detail (possibly the lower boiling point of acetic acid plays a part). In 
any case the results show that we have in general (so also in ascending 
chromatography) to reckon with two transport mechanisms in the vapour 
space, namely diffusion and rotation. 

From what follows below it is probable that in descending chromato- 
graphy the rotation current has a labile character. 

In our experience it is much more difficult in descending chromato- 
graphy to obtain regular and reasonably reproducible chromatograms 
than in ascending chromatography. For no obvious reason, perturbations 
may very easily occur. We found often that the three spots on one row 
(especially the lower spots) stand on an oblique line, though the front 
was horizontal. In fig. 8 and 9 the best chromatograms have been 
reproduced. That in fig. 9 shows still an irregularity, the middle lower 
spot lying higher than the adjacent one. For routine work we therefore 


cannot recommend descending chromatography. 
(To be continued) 
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7. Comparison of chromatograms run in the direction of smaller and larger 
ascension rate 

In Part I-III of this series and the preceding sections of the present 
Part we have always cut the Schleicher and Schiill 2043b in such a way 
that in running the chromatogram the mobile phase ascends in the 
direction of the smaller ascension rate. It was found empirically that in 
this way a somewhat greater separation of spots is obtained than when 
the chromatogram is run in the direction of the larger suction rate. 

We give here some examples with SiCly-treated Ederol 208 paper 
(small slit-feeding apparatus, di-isobutyl ketone—acetid acid —H,0= 
50:25:5, ascension height of the front =18 cm above the immersion line). 
This paper shows a marked difference in ascension rate into the two 
directions of the paper and the above-mentioned effect is clearly present 
here. Compare fig. 10 in which the pair A and B refers to acid SiCl,- 
treated paper of 21 hr age, C and D to neutralized SiCl4-treated paper and 
E and F to washed out SiCly-treated paper 1). 

In the survey below we give for each chromatogram the ascension time 
in minutes, the Ry,-values (mean of $ % and 1 % R) of the upper spot (C) 
and the lower spot (LL). As a measure of separation we take the differences 
of the Ry, values of the C and LL spot. 


chromatogram | A B zm. a | E F 
ascension time (min) | 275 207 265 | 202 266 — 200 
Ry, of C spot 0.639 0.63 0.63 0.625 0.625 | 0.615 
Ry, of LL spot 0.35 0.385 0.35 0.38 0.385 | 0.395 
separation 0.28 0,246 0.28 0.245 || 0.24 | 0.22 
kind of paper | acid (21 hr age) neutralized | washed out 


1) We have also compared two strips of acid SiCl4-treated paper of } hr. age. 
Here the distinct effect is not present. But as such newly prepared SiCly-treated 
paper rapidly loses HCl to the air, the difference in ascension time (more than 1 
hour!) of the strips into the two directions, is no longer short compared to the 
ascension time. Hence the papers are no longer comparable. 


= F 


Fig. 10. Comparison of chromatograms of the Reference mixture run in the 

direction of smaller and larger ascension rates. Arrows give the direction of larger 

ascension rate. A and B on SiCla-treated paper, 21 hours after preparation; C and D 
on neutralized paper; E and F on washed SiCly-treated paper (Ederol 208). 


The results show the signs of a change in the slope of the gradient 
curve. As the Ry, values of C for the two chromatograms belonging to 
each pair are within the errors equal, the turning point of the gradient 
curve lies at or very near to the C spot (as was also the case with SiCl,- 
treated Schleicher and Schill 2043b, compare fig. 6). 

For that chromatogram of each pair on which the Ry, value of the LL 
spot lies lower, the gradient curve must proceed with less slope. 

A simple explanation offers itself from the distinct correlation between 
ascension time and separation in the above survey. For each pair the 
separation is greater as the ascension time is longer. With a longer ascension 
time there is more time available for transport through the vapour space, 
with as a result a greater decrease of the slope of the gradient curve, 
hence a greater separation of the spots. 

In this simple explanation the anisotropy of the texture of the paper 
does not play a role itself. Whether this anisotropy may really be neglected 
must be left undecided here. 


8. Influence of stirring in the vapour space on chromatography 
a) Large slit-feeding apparatus with stirrers 


The apparatus used in this section is given in fig. 11. It has practically 
the same dimensions as the large slit-feeding apparatus described in a 
former communication [6]. The manipulations and precautions to be 
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10 cm 


Fig. 11. Large slit-feeding apparatus with stirrers fastened in front and back 
wall of the glass-container. 


taken are exactly the same. The apparatus only differs from the formerly 
described one by the possession of two stirrers a (glass strips 1.4 x 19 em) 
fastened to an Akulon axis, turning on Akulon bearings, fastened in the 
front and back glass wall of the container. Each stirrer can be driven by 
a separate synchronic motor standing about 30 em beside the chromato- 
graphic apparatus '). The synchronic motors are so placed that the two 
stirrers turn in opposite directions. Chain transmission ensures a constant 
speed of stirring. The stirrers make 50 revolutions per minute. After 
the paper sheet is brought in its place, the stirrers are set in motion. 
By viewing through the apparatus (with 4 Volt-incandescent lamp behind 
the back wall) one can make sure that the stirrers pass one another 
in an approximately vertical position. We now fill in the mobile phase 
and chromatography begins. Despite stirring the front rises horizontally 


out of the slit tube and remains practically horizontal during the 
experiment. 


') The synchronic motor should not stand too near the chromatographic 


apparatus, to avoid disturbances from radiation of the motor which becomes hand- 
warm while functioning. 
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So with the above apparatus we may compare chromatography without 
stirring (both stirrers are then placed vertically), with one stirrer turning, 
and with both turning. 


b) Chromatography of the Reference mixture on non- 
impregnated paper with di-isobutyl ketone —actie acid — 
H2,0=40:25:5 and 50:25:5 respectively 


We have run with each mobile phase three chromatograms to a front 
height of 18 cm above the immersion line, one without stirring, one with 
only the front stirrer moving and the third with both stirrers moving. 
The chromatograms obtained with one stirrer turning and two stirrers 
turning hardly differed. It therefore may suffice to give in the figures 12 
and 13 only the chromatograms obtained without stirring and with two 
stirrers turning. 

One observes a great influence of stirring; the spots le much farther 
apart than without stirring. The next survey gives the R,, values of the 
three spots. 


Mobile phase | stirrers turning | L+—C S LL + LC 
= 0.69 0.63° 0.42 
AQT 2555 I | 0.70 0.63 0.32 
2 | 0.72 0.635 0.31 
= 0.67 0.62 0.40 
Gras erd) i! 0.60 0.52° Or 
2 0.61 0.52 0.19% 


Fig. 12. Influence of stirring on chromatography of the Reference mixture with 


di-isobutyl ketone —acetic acid —H.O = 40:25:5 on non-impregnated paper. A 


stirrers at rest, B the two stirrers turning. 


Fig. 13. Influence of stirring on chromatography of the Reference mixture with 
di-isobutyl ketone—acetic acid—H2O0 = 50:25:5 on non-impregnated paper. A 
stirrers at rest, B the two stirrers turning. 


c) Chromatography of the Reference mixture on washed 
SiCl4-treated paper with di-isobutyl ketone —acetic acid 
— H,O=50:25:5 


The chromatograms (front height 18 cm above the immersion line) 
are given in fig. 14. The Ry,-values of base and upper margin of the spots 
are given in the next survey. 

The results are in principle the same as in the preceding paragraph: 
With stirring a greater mutual distance of the spots is obtained than 
without stirring. 


stirrers turning | C L LC+S5S LL 


d) Discussion 


In section 4 the greater spread of the spots which is obtained by 
increasing the volume of the chromatographic container (see section 3) 
has been discussed. It was assumed that the greater space allowed a more 
effective transport by diffusion of acetic acid evaporated from the lower 
part of the wetted paper and taken up by the higher part of the wetted 
paper, similarly of di-isobutyl ketone in the reverse direction. Thus the 
gradient curve becomes less steep and as a consequence the distance of 


the spots on the chromatogram increases with increasing volume of the 
container, 


Fig. 14. Influence of stirring on chromatography of the Reference mixture with 
di-isobutyl ketone—acetic acid—H2O0 = 50:25:5 on washed SiCly-treated paper. 
A stirrers at rest; B the two stirrers turning. 


Accordingly, in an apparatus with vapour space, diffusion is the limiting 
factor in the chain of events (evaporation — transport through the vapour 
space by diffusion—condensation) which will determine the degree of 
diminished slope of the gradient curve compared to the very steep slope 
in an apparatus without air space (though such an apparatus could not 
be realized). 

The slow transport through the vapour space which takes place in a 
container without stirring is probably not only due to diffusion. In section 
7 we have obtained indications that a slow transport by convection in 
the vapour space also plays a role. 

So it must be expected that a greater decrease of the slope of the 
gradient curve will occur with stirring than without stirring because 
diffusion and rotation as limiting factors are cancelled out by stirring. 
Consequently the spots will stand farther apart with stirring than without 
stirring. This is indeed what has been found in this section. Compare 
figs. 12, 13 and 14. 

Before we may safely accept the above given explanations we shall 
check in the next section that indeed the gradient curve proceeds with 
smaller slope while stirring in the vapour space, than without stirring. 


9. Influence of stirring in the vapour space on the gradient curve 
a) Chromatograms 

For the experiments in this section SiCl,-treated paper has been used, 
which has been exposed sufficiently long to the air (90 hours or longer) 
to have lost practically all HCl. 
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Two chromatograms on paper sheets 17x28 cm have been run with 
di-isobutyl ketone — acetic acid —H.O =50:25:5, allowing the front to 
rise to 22.5 em above the immersion line, the one without stirring in the 
vapour space, the second with the two stirrers turning. The chromato- 


grams are given in fig. 15. 


Clea ereverere 


Fig. 15. Influence of stirring on chromatography of the reference mixture with 
di-isobutyl ketone—acetice acid—H20 = 50:25:5 on old SiCly-treated paper (practi- 
cally HCl-free). A stirrers at rest; B the two stirrers turning. 


The following survey gives the Ry, values: 


without stirring 


with stirring 0.38 0.285 0.19 


b) Gradient determinations 


For these determinations we used paper sheets, on which besides the 
immersion line and the future front line (22.5 em above the immersion 
line) were drawn further horizontal pencil lines on 4.5 em mutual distance. 
After rise of the front to the front line, the paper was cut rapidly into 
five zones, each of 174.5 em, which were analysed as given in detail 
in Part HI of this series [2]. The immersion zone (17 x 0.5 em) is of course 
cut off and disregarded. 
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The results have been given in Table IV. The data of the table have 
been used in fig. 16 to draw the gradient curves. 


TABLE IV 
Influence of stirring in the vapour space on the acetic acid gradient in practically 


HCl-free SiCly-treated paper 


ee 


Without stirring Stirring with two stirrers turning 
- | 
position wet strip | acetic | acetic | position wet strip | acetic | acetic 
of zone; | —dry strip | acid acid of zone; | —dry strip | acid acid 
Op OL ire | mg | mg % Thy GOH Ae; mg mg ve 
0 — 19.9 811 311 388 OFS os 782 252 322 
19.9— 40.2 746 249 334 |19.7— 39.5 711 206 208 
40.2— 60.2 686 174 254 | 39.5— 60.1 736 202 274 
60.2— 80.3 533 ee 21° | 60.1— 80.0 614 170 27? 
80.3-100 313 64 207 | 80.0-100 380 114 30° 
ve acetic acid 40-1 pci’ 
304 Sis a 


a 
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Fig. 16. Influence of stirring on the acetic acid gradient in the mobile phase 
wetting the paper (same mobile phase, ascension height of the front and kind of 
SiCly-treated paper as in fig. 15). A without stirring; B with two stirrers turning. 


c) Discussion 
x) Change in slope of the gradient curve as a result of stirring in the vapour space 


When we first consider the chromatograms given in fig. 15, one sees 
that with the paper used a similar strong effect of stirring on the spread 
of spots is present, as we met in the preceding section. 

It has been argued in preceding sections that the great influence of 
stirring in the vapour space must be due to an elimination of diffusion 
and rotation as limiting factors as a result of which the gradient curve 
should proceed less steeply. 

The gradient determinations in this section confirm this argument. 
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The gradient curve has indeed a much smaller slope than without stirring. 

A peculiarity is that the transport through the vapour space has been 
so efficacious that the gradient curve with stirring proceeds nearly 
horizontally to the right of the intersection point with the gradient curve 
without stirring, or even (when this is no experimental error) goes here 
through a flat minimum. Compare fig. 16, B. 

That this remarkable shape of the gradient curve may arise, must be 
due to the fact that the liquid wetting the paper decreases strongly from 
immersion line to front (see columns 2 and 6 in Table IV). When (as a 
result of stirring) the vapour in the whole vapour space has nearly the 
same composition, the acetic acid concentration may increase with 
greater rate just in the last zone below the front, where the least amount 
of liquid is present. 


B) Position of Ry-values on the two gradient curves 


The results in the preceding parts II and III of this series led to a 
picture of the mechanism of the chromatography, namely adsorption- 
chromatography combined with gradient elution. The base of a spot 
should lie at a particular composition of the mobile phase (as an appro- 
ximate measure of which we take the acetic acid concentration) where 
it is eluted, while higher up the phosphatide is readsorbed. 

When by external means we are able to change the gradient curve, it 
should appear that the Ry, values of the spots lie on the gradient curves 
at equal acetic acid concentration. 

In section 3 we found that the slope of the gradient curve decreases 
with increase of the volume of the chromatographic apparatus. The 
gradient determinations in this case, however, are subject to great expe- 
rimental errors (the zones analysed being only 4x34 em) to allow a 
control of the above discussed position of the Ry, values on the gradient 
curves, 

It would seem that we are here in a more favourable condition. The 
error will be less since the zones analysed are 17x44 em, and gradient 
curves do not proceed so steeply). 

In fig. 17, we have redrawn from fig. 16 the gradient curve without 
stirring (left diagram) and with stirring (right diagram). On each curve, 
points have been indicated with abscissa corresponding to the Ry, values 
given in the survey of paragraph a). 

We see that in the case of points LL, LC+S8 and L the spots on the 
chromatograms are situated within experimental error at the same 
acetic acid concentrations. Compare the three horizontal broken lines 
which pass the points on the two gradient curves at distances of at most 
one weight % acetic acid (the experimental error). 

It is not possible to draw such a horizontal line through the C-spots. 
Must we conclude that the mechanism proposed in Parts If and III of 
this series for chromatography without stirring, only holds for LL, 


acetic acid 


Slo 


acetic acid 
S/o 


| 
20 40 60 80 100 20 40 60 80 100 


Fig. 17. Position of the R;;-values of the spots on the two gradient curves of fig. 16. 
A without stirring; B with two stirrers turning. 


LC+S8 and L, but does not apply for C? This conclusion is not valid, for 
with stirring the gradient curve has its deepest point at an acetic acid 
concentration which lies higher than the elution concentration of C 
(given by the point on the left gradient curve). Hence with stirring the 
C spot on the chromatogram is no longer adsorbed on the paper, but 
moves in an eluted state. 

Here the question arises as to the mechanism in which chromatography 
of the C spot occurs in the case of stirring. Another question also arises: 
should we regard the mechanism of chromatography with stirring of 
the other spots (LL, LC +L) still as adsorption chromatography combined 
with gradient elution ? 

We shall not go into these questions here. They will be dealt with experi- 
mentally in the next following part of this series. 


10. Influence of temperature 


a) Experimental 


The experiments have been made on non-impregnated Ederol 208 paper in the 
large slit-feeding apparatus provided with stirrers (see fig. 11). At each temperature 
investigated (4°, 20°, 25°C) two chromatograms (ascension height of the front 
18 cm) have been run with di-isobutyl ketone—acetic acid—H20 = 40:25:5, of 
the Reference mixture (5 mm 4} % solution); one chromatogram without stirring, 
the second chromatogram with both stirrers turning. The chromatograms obtained 


are given in fig. 18 and fig. 19. 


6b) Results 


The results of the experiments without stirring (fig. 18) show clearly that with 
increasing temperature the mutual distance of the three spots increases. Hence we 
conclude that with increasing temperature the slope of the acetic acid gradient 
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decreases. A remarkable detail is that the separation of the two upper spots 
(L + C; 8) which has taken place going from 4° to 20°, does not increase (even 
becomes less) going from 20° to 25°. Whereas the S-spot remains on its place from 
4° to 20°, it shifts in upward direction going from 20° to 25°. 

The results of the experiments with stirring (fig. 19) show in principle also an 
increase of spreading of the spots with increasing temperature, at least when the 
distance of the upper spot (L + C) to the lower spot (LL + LC) is considered. 

The effect of increase of the temperature is less pronounced than in the experi- 
ments without stirring. This is plausible as with stirring the gradient already has 
a smaller slope and it cannot be made very much smaller by another factor. 

In fig. 19 we see a similar detail as was mentioned above for the chromatograms 
obtained without stirring: 

Going from 4° to 20° the 8S spot remains on its place, the L + C spot moves 
somewhat upwards and the LL + LC spot somewhat downwards. So these are 
signs of a regular decrease of the slope of the gradient. 

Going from 20° to 25° all spots move in upward direction. It seems therefore 
that the gradient curve shifts in upward direction. 


c) Discussion 
a) Complexity of the temperature influence 


When it is tried to argue which will be the influence of increased temperature 
on the mutual spreading of the spots, that is on the change in the gradient, we 
meet with the difficulty that so many factors are involved. We may divide them 
into four groups. For each factor we shall design in the survey below the influence 
of increased temperature with + when the gradient by this factor alone would 
become less steep, that is the spreading of spots becomes greater. The sign | will 
be used in the reverse case. 


I. Ascension time 
a) surface tension (driving force of ascension) decrease with temperature, 
hence by this factor alone the ascension will become slower and more 
time will be available for the transport mechanisms, so +. 
b) viscosity (resistance for ascension) decreases with temperature hence 
available time for transport is shorter, so |. 


Il. Partial vapour pressures 


c) vapour pressure increases with temperature, hence all other factors being 
constant, the amount transported through the vapour space is greater, 
hence +. 


II. Transport mechanisms 


d) diffusion in the vapour phase increases with temperature, so in the same 
time more transport, hence +. 

e) driving force of the rotation will increase with temperature as the down- 
ward moving gas layer near the paper (section 6) will increase in density 
as a result of increased vapour tension, Hence more transport, so +. 


f) viscosity of gases increases with temperature, hence slows down the con- 
vection by rotation in the vapour space. Hence less transport, so ae 


IV. Adsorption and elution 


In the above groups I, If and III it was assumed that in the range of 
temperature investigated, adsorption and elution are not affected by increase — 
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of the temperature. When this is not true it is difficult to argue the effect 
of temperature on both and from this the influence the spreading of the spots 
is difficult to foresee, and this will not be tried here. 


B) Temperature influence between 4° and 20° 


When in fig. 18 the chromatograms at 4° and at 20° are compared, there is no 
indication that we must seriously take into account factors belonging to group IV 
of the above survey. The S spot remains practically on its place, as was also the 
case with non-impregnated Schleicher and Schiill paper and the same mobile phase 
(see fig. 6). Hence the turning point of the gradient remains on its place going from 
4° to 20°. So we shall only consider the groups I, II and III of the above survey. 

As to group I, it appears from the time in minutes written on each chromatogram, 
that from 4° to 20° the ascension time distinctly decreases, so the influence of 
temperature on factor 6 preponderated over that in a. 

So when in the shorter time available a greater mutual distance of spots is obtained 
at 20° than at 4°, in the experiments without stirring (fig. 18) it follows that the net 
result of increased temperature on the factors mentioned in groups II and III is 
a distinct decrease of the slope of the gradient curve. 

In the experiments with stirring (fig. 19) the slow transport mechanisms are 
removed. So the factors mentioned in group III need not to be considered for 
explaining the increase of mutual the distances of the spots. So group II remains, 
the increase of the partial vapour pressures with temperature. 

In the experiments without stirring, the transport mechanisms are the limiting 
factors, which is evident from the great influence of stirring. It follows that the 
net effect of increasing temperature on factors d, e and f leads to a decrease of the 
slope of the gradient curve, and as factor f would in principle lead to the reverse 
result, the increased rate of diffusion (d) and of rotation of the vapour mass by 
factor e, preponderate. 

Resuming we may say that without stirring the increased distance of the spots 
with increase of the temperature results from increase of vapour tension and of the 
increased rate of the transport mechanisms. 


y) Temperature influence between 20° and 25° 


It is not so easy to explain the changes which take place in the chromatograms 
developed at temperatures from 20° to 25°. It seems likely that in this temperature 
interval factors involved in group IV can no longer be neglected, since there are 
signs that the gradient curve shifts in upward direction (compare paragraph 6). 


11. Summary 


1. With specially constructed chromatographic apparatuses differing 
widely in available air space the resolution of the spots of the Reference 
mixture (C=cephalin, L=lecithin, LC and LL=their lysoproducts and 
S=sphingomyelin) has been investigated. Non-impregnated and SiCly- 
treated paper were used, with as a mobile phase di-isobutyl ketone — acetic 
acid — H20. The acetic acid gradient in the mobile phase wetting the paper 
has been investigated. 

2. Resolution of spots is small with a very small available air space; 
the gradient curve in this case is very steep. 

With increasing air space the resolution of spots increases and the 
slope of the gradient curve diminishes. The results support the conclusions 
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reached in Part III of this series, i.e. the location of a chromatographic 
spot indicates a particular acetic acid concentration on the gradient 
curve, at which the spot is eluted at its base with maximal rate. 

3. In the not realizable case of an apparatus without air space, the 
gradient curve will consist of three parts; a horizontal part from immersion 
line to a certain fraction of the distance front to immersion line, an 
extremely steep descent and a second horizontal part up to the front and 
with acetic acid concentration zero. The origin of the gradient is attributed 
to adsorption of acetic acid to the paper. Elution of spots occurs by 
adsorption displacement. The sequence of the spots on non-impregnated 
paper and on SiCly-treated paper can be explained from this point of view. 

4. The air space in a chromatographic apparatus allows acetic acid 
transport by diffusion through the vapour phase from the lower part 
of the wetted strip to the upper part. Similarly di-isobutyl ketone transport 
occurs from the upper part of the strip to the lower part. So the slope of 
the descending part of the gradient curve is lowered and as a consequence 
the mutual distance of the spots on the chromatogram increases. The 
effect is greater as diffusion is more effective, for instance in an apparatus 
with larger air-space. 

5. With a specially constructed slit-feeding apparatus descending and 
ascending chromatography have been compared. 

The results (greater separation of spots in descending chromatography) 
point to an additional transport mechanism in the vapour phase. A 
downward movement of the vapour occurs along the wetted strip and 
upward along the glass wall of the container. In descending chromato- 
graphy the transport of acetic acid from the slit-side part to the frontside 
part of the wetted strip by means of this rotation current of the vapour 
is much more efficient than in ascending chromatography. Hence the 
slope of the gradient curve is smaller here than in ascending chromato- 
graphy. In practice descending chromatography is not recommended. 
It is difficult to obtain acceptable and reproducible chromatograms. 

6. The recognition that the vapour space in a chromatographic 
apparatus partly determines the degree of separation of spots allows to 
foresee, or helps to explain, the influence of other factors. 

7. Witha given chromatographic paper, separation is somewhat greater 
when the mobile phase ascends in the direction of the slower suction rate 
than in the direction of faster suction rate. The ascension time is longer 
in the first mentioned case, so here is more time for diminishing the 
slope of the gradient curve. 

8. One should expect that stirring in the vapour space will increase 
the resolution of spots. Stirring will eliminate as a limiting factor the 
slow transport mechanisms in the vapour space (diffusion; rotation). 
Hence the slope of the gradient curve will be diminished in the case of 
stirring. 

These expectations have been confirmed by means of a specially con- 
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structed, large slit-feeding apparatus, provided with stirrers. The experi- 
ments show indeed an increased resolution and a strongly dimished slope 
of the gradient curve. 

An unforeseen but easily explained additional detail is that the gradient 
curve with stirring shows a shallow minimum, the implications of which 
have been discussed for the mechanism of chromatography with stirring. 

9. Itis difficult to predict the influence of temperature on the separation 
of spots in general. The net result will be composed of the effect of 
temperature on all partial processes involved in chromatography. 

Without stirring the separation of spots increases in the range 4°-20°. 
It is argued that this is due to the preponderating influence of increased 
partial vapour tensions, increased rate of diffusion, and increased intensity 
of rotation. 

10. With stirring the separation of spots also increases with temperature. 
The limiting transport mechanisms are eliminated by stirring so the effect 
must be due to the preponderating influence of increased partial vapour 
tensions. 

11. Our experiments suggest that in the range 20°-25° the influence 
of temperature on adsorption and elution begins to become apparent. 


Department of Medical Chemistry, 
University of Leyden 
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PHYSICAL CHEMISTRY 


GRAPHITIC OXIDE 
VI. SOME PHYSICO-CHEMICAL PROPERTIES 


BY 


J. H. DE BOER anp A. B. C. VAN DOORN 


(Communicated at the meeting of September 24, 1960) 


1. Introduction 

In the preceding paper (communication V) [1] we discussed the swelling 
of graphitic oxide by the adsorption of water. From this discussion 
resulted the remark that perhaps the réntgenographically determined 
c-spacings might be wrong in so far as these spacings should not represent 
the actual situation and should only be average values. 

In order to investigate this suggestion we determined the density of 
anhydrous graphitic oxide and of some water containing samples using 
benzene as imbibition liquid. For it is known that graphitic oxide does 
not swell in benzene and in this manner we can thus measure the changes 
in volume of graphitic oxide caused by the penetration of water between 
the hexagon layers. 

Moreover we determined the density of anhydrous graphitic oxide in 
a number of other imbibition liquids. 

The final conclusions from these experiments and from those described 
in the preceding paper will be discussed. 


2. The density of graphitic oxide in benzene 

The density of graphitic oxide was determined in a pycnometer (about 
4 ml) in the following way. About 200 mg graphitic oxide were immersed 
by 1 ml of dry benzene and afterwards the sample was de-aerated by 
careful evacuation at a pressure of about 18 mm Hg. Thereupon the 
pyenometer was filled with benzene and placed in a thermostat at 25° C, 

Figure 1 gives the change of the specific volume of graphitic oxide (Vs) 
in relation with the water content of the samples. It is seen that the 
specific volume reaches a minimum value at a water content of 4.2 %; 
afterwards the specific volume increases linearly with increasing water 
content. 

As a matter of fact the specific volume is calculated per gram material, 
so the amount of graphitic oxide decreases with increasing water content. 
It is, however, much more convenient to relate the volume of 1,000 mg 
anhydrous graphitic oxide (Vs) plus the volume of the adsorbed water 
(Va) with the water content. By way of example we recalculate two points 
of figure 1 in this manner: 


Vg (cm) 
0,60 


0.55 
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Fig. 1. Vs is specific volume in benzene of graphitic oxide containing known 
amounts of water. % HzO has been calculated on anhydrous graphitic oxide. 


One sample graphitic oxide contains 1.6 °%, HzO (calculated on dry 
material) and the specific volume is 0.481 ml. One gram of this product 
thus contains 984.3 mg graphitic oxide and 15.7 mg H2O. Thus 1,000 mg 


anhydrous graphitic oxide and 16.0 mg H2O together have a volume of 


1,016 - ae 
1,000 x 0.481=0.489 ml. 


The second sample contains 4.3 °%, water and the specific volume is 
0.468 ml. One gram of this sample consists of 958.8 mg anhydrous 
graphitic oxide and 41.2 mg water. Thus 1,000 mg anhydrous graphitic 


1,000 eon 
1,023 <x 0.468=0.488 ml. 


In this manner all experimental points of figure 1 are recalculated and 
are shown in figure 2 (curve A). We see now that graphitic oxide indeed 
is capable of adsorbing 4.2 °, water without any swelling. There must 
thus be free space in the graphitic oxide lattice. In figure 2 is also drawn 
curve B which relates the water content with the sum of the volumes of 
1,000 mg anhydric graphitic oxide (0.489) and of the amount of free 
water (density=1), corresponding with the amount of adsorbed water. 

The slope of curve B is not equal to that of the ascending part of 
curve A. However, this is caused by the fact that graphitic oxide, 
containing more than 10% water, looses part of the water which is 
taken up by the dry benzene. This was experimentally verified in the 
following way. After the determination of the density of graphitic oxide, 
the upper layer of benzene was withdrawn with a pipette, and the water 
content and the density of this benzene were both determined. The 


oxide plus 43 mg water together have a volume 
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specific volume of the graphitic oxide was now again calculated, using 
the changed density of the benzene; moreover the water content of the 
graphitic oxide was corrected for the loss of water. This specific volume 
of graphitic oxide, calculated per 1,000 mg anhydrous material, is given 
in figure 2 by curve C. We see that curves B and C in this figure are 
completely parallel. We therefore must conclude that the density of the 
adsorbed water is indeed equal to 1. 


Vg +V (cm?/g) 
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Fig. 2. Vs + Va is volume of 1 gram of graphitic oxide plus the volume of the 
amount of water adsorbed. ° H2O has been calculated on anhydrous graphitic oxide. 


3. The density of graphitic oxide in several imbibition media 

The technique of these determinations is the same as already described 
in the preceding section. 

Preliminary experiments were carried out with different graphitic 
oxides viz. 3H’, 6H’ and with 3X. No difference in specific volume between 
these samples was observed using water, benzene or methanol as imbibition 
liquid. Also in this respect the different graphitic oxides all show the 
same behaviour. 


In table I are given the specific volumes of graphitic oxides in a number 
of imbibition media. 
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TABLE I 
———————— 
Imbibition medium | Specific volume ml/g 

ee eee 
methanol 0.383 

1—4 dioxane 0.399 
ethanol 0.400 
pyridine 0.412 

water 0.423 
nitrobenzene 0.446 

helium (corrected) 0.465 
reciprocal réntgen-density 0.466 
benzene 0.488 
carbontetrachloride 0.515 


The specific volume of graphitic oxide in helium?) is 0.475 ml/g. It 
was recently shown [2] that the true specific volume of an adsorbent is 
best approximated by correcting the specific volume in helium for the 
distance to which the helium atoms approach the surface of the adsorbent. 
If we take for this distance the radius of the helium atom, calculated 
from the b-term of the Van der Waals equation, i.e. 1.3 A, then this 
sterical correction has the value: S x 1.3x 10-8 ml. Here S is the outer 
surface of the adsorbent; in our case S=74.8 m2/g, so the correction has 
a value of 74.8 x 104 x 1.3 x 10-8=0.010 ml/g. The corrected specific volume 
in helium becomes thus: 0.475—0.010=0.465 ml/g. This value agrees 
very well with the reciprocal X-ray density, i.e. 0.466 ml/g. The latter 
was calculated as follows: one gram C7O4He, having a unit weight 150, 


contains EUS ee 
150 


corresponding with 3.5 hexagons, is 3.5 5.25=18.4 A2 whereas the 
distance between two hexagon layers is 6.3 A. Thus the calculated 
specific volume becomes: 


units C7O,H»e; the surface of one unit C7O4He, 


8.02 x10” 18.4 x 6.3 x 10-%4= 0.466 ml/g. 
150 

The agreement between the calculated and the experimentally observed 
value (in helium) of the specific volume seems to indicate that graphitic 
oxide has retained the graphite character in the hexagon layers during 
the oxidation. On the other side this graphite character is improbable 
in view of the hydroaromatical character which we must assume in order 
to understand the flexibility of the hexagon layers. In a future paper 
we hope to return to this discrepancy. 

The difference between the specific volume in water and in helium 
(corrected) is 0.042 ml. This is in perfect agreement with the fact (known 
from the réntgenographic investigation (communication V) [1] and shown 
above) that graphitic oxide can adsorb the first 4.2 % of water without 


1) Our thanks may be expressed here to Mr. P. ZWIETERING who determined 


the specific volume of graphitic oxide in helium for us. 
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any swelling. This water, which penetrates into the lattice, is apparently 
taken up by free space between the hexagon layers. Gases, including 
He and liquids which do not penetrate into the lattice, cannot be taken 
up by this space. 

2emarkable is furthermore the very low specific volume of graphitic 
oxide in methanol and in dioxane. The dioxane molecules, which are 
much larger than the water molecules, evidently fill a larger space between 
the hexagon layers than the water molecules. The reason of this fact is 
perhaps that the water molecules are adsorbed on specific centres. One 
could imagine that these specific centres are to be identified with the 
carbon-oxygen complexes. 


4. Discussion of the results 

There is a discrepancy between the X-ray data of the swelling in water 
(communication V [1]) and the specific volume data of the swelling in 
the present communication. According to the present data 4.2 % of 
water is taken up by graphitic oxide without any swelling, whilst any 
more water, taken up by swelling, adds its own volume (specific volume 
of this water being 1.0) to that of the graphitic oxide. The X-ray data 
agree with the lack of swelling by the first 4.2 ° of water. From 4.2 % 
to about 16 % of water, however, X-ray data indicate only a linear 
swelling of the distance between hexagon layers from 6.3 A to 6.8 A. 
The absorption of 12 % of water, therefore, results only in an increase 
of this distance of 0.5 A, whilst it should be 2.1 A. This latter figure is, 
indeed, found on further swelling, as already indicated in commu- 
nication V [1]. 

It has already been stated in communication V that 16.2 °% of water 
corresponds with the so-called unimolecular layer of the B.E.T.-theory; 
it is the Vm value. This Vm value, however, is too small for a unimolecular 
layer of water molecules, covering the total surface of all hexagon layers. 
As already stated above these water molecules may be bound to specific, 
hydrophylic centres. 

From the conclusion of the investigation of the oxidative properties 
of graphitic oxide (communication IV [3]) we know that the bruto 
formula C7O,H» can be written as: 


{C7(= 0)HjO3H = {C7(OH)}O3H 


Furthermore it appeared that the hydrogen atom of the O3H combination 
is bound to an oxygen atom in either an OH or an OOH group. In this 
discussion we will tentatively call this combination an “OH” group. 
It seems justified to identify the Vm value of the water-sorption isotherm 
with that quantity of water required to hydrate all hydrophilic, say 
“OH”, groups. The absorption of water may thus take place on two 
different groups viz. the enolic OH group and the above mentioned 
“OH” group. The swelling of graphitic oxide takes place moreover in 
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two steps. The first 4.2 % H2O (ie. 0.36 mol H2O per unit C7O4H2) are 
taken up without any swelling; afterwards the adsorption of water is 
accompanied by swelling. When 16.2% H2O is adsorbed (i.e. the 
Vm-value) there is exactly 1 mol HsO more adsorbed since the first 
4.2 % H2O; for 16.2 % HsO corresponds with 1.35 mol H,0. 

This, however, does not mean that really every “OH” group is hydrated, 
for the Vm-value corresponds with a statistical monolayer. As a matter 
of fact those “OH” groups which are situated near the edges of the 
lamellae are hydrated more easily than the “OH” groups in the middle 
of the lamellae. However, if the first “OH” groups (i.e. the most accessible) 
are hydrated, these hydrated “OH” groups can serve as second adsorption 
centres. This is the more probable if one keeps in mind that the attraction 
forces between the hexagon layers are essentially Van der Waals forces. 
These are, however, inversely proportional to r3 if r is the distance 
between the layers. As this distance is increased by the adsorption of 
the water molecules, the energy required for the second molecule of water 
may be smaller. It is true the real adsorption energy will be smaller for 
the second molecule but the total adsorption energy (i.e. real adsorption 
energy minus swelling energy) may be more favourable than if the second 
water molecule penetrates further into the lattice in order to hydrate 
the “OH” groups lying further inwards. These same arguments hold 
for the adsorption in the third and higher layers. 

A similar process of gradual penetration from the outside to the inside 
of the layer lattice is found with Mg-vermicullite. WaLKrr [4] observed 
that the penetration of water causes an increase of the layer-distance 
from 11.6 A to 13.8 A. This increase starts from the outside and is 
followed by a second increase from 13.8 A to 14.4 A, even before the 
first increase has reached the centers of the platelets. The gradual swelling 
here proceeds in distinct steps, whilst in graphite oxide it proceeds 
continuously. This latter behaviour is caused by the great flexibility of 
the hexagon layers (see communication V [1]) in contrast to the rigidity 
of the layer of the Mg-vermicullite lattice. One may visualize that in the 
case of graphitic oxide the penetrating water causes the layers to bend, 


as is shown schematically in figure 3. 
This process of bending of the lamellae kas eae 
proceeds so far until it is energetically more ee Se ae 


favourable for the water molecules to penetrate 


further into the lattice rather than to separate = ~~ 


the layers more from each other. Se ee 


Fig. 3. The swelling of a particle of oS 
graphitic oxide by the uptake of water. 
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Schematically this process of swelling is given in figure 3. It is evident 
that we obtain a continuous spectrum of c-spacings. Réntgenographically 
we only determine an average of the mostly occurring ¢-spacings. The 
extreme values of the c-spacing, which occur only in relatively small 
quantities but which, however, correspond with a large part of the 
adsorbed water, are not determined réntgenographically and thus the 
e-spacing, determined in this manner, is too small. 

This process of adsorbing and swelling goes on in the manner described 
till graphitic oxide has a water content of 16.2 °%%; because from this 
point onwards the swelling in the direction perpendicular to the laminae 
increases linearly with the water content, and one has to conclude that 
from this stage onwards swelling and further penetrating into the layer 
lattice proceed equally well. For the average increase of the c-spacing 
from 16.2 to 32.4 °%% H2O (i.e. the forming of the second statistical layer 
of adsorbed water molecules) is 2.1 A which is a value also found with 
the swelling of montmorillonite [5, 6]. 

The Vm-value from the water vapour isotherm is thus probably best 
identified with the number of hydrophilic centres in graphitic oxide. On 
these centres the water molecules are bound by localized adsorption. 
This explains the large difference between the “‘water-surface” and the 
theoretically calculated lamellar surface of graphitic oxide viz. 585.1 and 
738.4 m? per gram respectively (see our preceding paper [1)). 

From the density in helium and from the specific surface we can 
calculate the number of hexagon layers per crystallite of graphitic oxide. 
Suppose there are n particles per gram; the volume of 1 gram is: 
n x arzh = 0.465 « 1024 A’, The surface per gram is: n x 22r2= 0.70 x 1020 A. 


From this it follows h= 130 A corresponding with 130 is about 20 hexagon 


6.3 
layers. 
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PALEONTOLOGY 


DICERORHINUS HEMITOECHUS FALC. IN THE NETHERLANDS 
BY 


H. LOOSE 


(Communicated by Prof. I. M. vAN DER VLERK at the meeting of September 24, 1960) 


During sand dredging operations in the Zwarte Water, near Westerveld 
(province of Overijssel), a rhinoceros skull was found in May ’58. 

It had the general form of Coelodonta antiquitatis, with some important 
differential characteristics. These all pointed to Dicerorhinus hemitoechus 
Falc., a species not yet recorded from the Netherlands. 

Of this well preserved skull (a few, not very important, fragments and 
the mandible are missing) a summary description, followed by some 
general remarks on the species, is given here. 


(See fig. 2-5). The nostrils extend farther backwards than in C. anti- 
quitatis. While in that species the posterior border lies approximately 
between P? and P4, in this skull the nostrils end above the posterior rim of P*. 

The rugosity of the nasal horn base is considerable. The frontal horn 
base is seen in lateral view by a slight curve of the frontals only. 

While C. antiquitatis is known for its complete bony nasal septum, 
in this skull the bony septum is incomplete in the middle. There are 
traces of a cartilaginous part where the bony septum ends. 

The zygomatic arch shows a double bend. The maxillar part extends 
in a straight line to the rear; the jugal turns up and backward at an 
angle of 30°; at the glenoid fossa, the arch runs downward toward the 
median sagittal plane. 

The articular plane on the anterior rim of the glenoid fossa is 11.5 em 
wide, measured perpendicular to the m.s. plane. In two skulls of C. anti- 
quitatis values of 9.5 and 10 cm were found. 

Behind the fossa there is a well developed postglenoidal processus, 
forming the anterior rim of the auditory meatus (a triangle pointing 
downward with one angle), the paroccipital processus providing the 
posterior rim. 

The parietals are deeply constricted by the cristae temporales (minimum 
width 3.4 cm; in skulls of C. antiquitatis 7.5-12 cm). 

The occipital crest is not as wide as the base of the skull. This causes 
the occipital view of the skull to be a trapezium. In C. antiquitatis this 


part is a rectangle. 
The occipital crest of this skull is 16.5 cm wide, in C. antiquitatis about 


YB} Cita. 
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In the occipital crest of C. antiquitatis a median fossa is an exception, 
sometimes there is a median protuberance. The present skull has a marked 
median fossa. 

The foramen magnum has the shape of an ellipse, with its shorter axis 
in the median sagittal plane. In C. antiquitatis the foramen is higher, 
with the shape of an ovoid or triangle. This is caused in part by a difference 
in the orientation of the foramen magnum. In the skull from the Zwarte 
Water the angle between the vertical axis of the foramen magnum and 
the roof of the mouth is 119°. In ©. antiquitatis (ZEUNER, 1934, 19 ind.) 
the maximum was 107°, the median 95°, the minimum 83°. 


In 1822 Cuvrer described a rhinoceros skull without bony nasal septum 
under the name Rh. leptorhinus. De CuristoL thought (an opinion entirely 
founded on reproductions of drawings of the type specimen) that a 
complete bony septum had been present. Therefore he described the 
skull without septum, found in 1835 near Montpellier (dép. Hérault), 
under the name Rh. megarhinus. 

In 1846 Owen described part of a skull from Clacton-on-Sea (Essex) 
with septum as identical with Rh. leptorhinus Cuv., thereby endorsing the 
judgement of Dr CurisTou. Finally FALCONER in 1860 created the species 
Rh. hemitoechus for the Clacton skull and some others with an incomplete 
bony septum, after showing the impossibility of maintaining the name 
Rh. leptorhinus, given by CuvieR to a skull without septum. 

The species, described by JAGER originally in 1839 as Rh. kirchbergensis 
and subsequently in 1841 as Rh. Merckii, is rejected by FALCONER as 
founded on insufficient data (2M, M38, 1 lower molar) and provisionally 
identified as belonging to Rh. megarhinus. 

Now confusion is complete. Only one species remains well defined, 
despite frequent changes in the name of the genus: Coelodonta antiquiratis 
Blum. Much has been published on the other pleistocene species in the 
century after FaLconER, with usually every author taking the opposite 
view of his predecessor. 

In this battle of systematics, Rh. leptorhinus and Rh. megarhinus 
appear as rather vague and loosely defined names of doubtful value. Most 
of the heavy fighting has been (and in some countries still is) over the 
question: are hemitoechus and kirchbergensis separate species? 

Complete, well preserved skulls are rare. The efforts to differentiate the 
species D. etruscus—hemitoechus—kirchbergensis are therefore mostly based 
on differences in dentition. The fact that this is extremely difficult for 
D. hemitoechus and D. kirchbergensis is no proof that these species are 
actually one. When the complete skulls are compared, all evidence is to 
the contrary 1). 

F. Znuner has shown 2) a connection between anatomy of the skull 


1) STAESCHE °4]. 
2) ZEUNER ’34. 
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and ecology in the Rhinocerotidae. To this end certain standard measu- 
rements were treated statistically. 

Owing to lack of material ZeuUNER was unable to apply his method to 
D. hemitoechus. In the present publication. measurements according to 
ZHUNER (and an additional dimension) are given for 8 crania of D. 
hemitoechus (6 in the British Museum, Natural History, the skull from 
the Zwarte Water and the cranium from the Stuttgart Staatl. Museum 
fiir Naturkunde). 


8. Collet 


Fig. 2. Leiden No. 93302 occipital view. 
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y 83 95 107 95 19 110.” Tis > 120 ~ Lig. “S | 82 95.5 109 95.5 4 
m 40 53 80 55 34 28 37 44 37 6 | 4!1 55 66 54 5 
58 —14 —1.5 +27 +1 32 |—21 14 —9 -14.6 6 |—42 17.5 13 —22.5 4 
Z 75 88 120 92.6 12 34 49 60 47.75 8 | 


r 


q 


The measurements for Coelodonta and D. kirchbergensis have been taken from ZEUNER °34. At the ext 


right are the values for the skull from the Zwarte Water 


Abbreviations 


1, st and z have been measured in mm, all others in degrees. All 
measurements, except z, have been taken in the median sagittal plane. 
The m.s.p. is not mentioned in the definitions, to keep them as simple 


as possible. (See also fig. 1). 


| Distance rhinion-basion (length from tip of nasals to lower margin 
of foramen magnum). As the development of the occipital crest 
varies for species, age and sex, the basal length of the skull has 


been taken as standard. 


st Incision in the middle of the occipital crest. 

Oo Angle between opisthion + occipital crest and the parietals. 

i Angle between parietals and hindmost point of the occipitals. 
n Angle between parietals and tangent of hornbases. 


= height opisthion-extended 


roof of 


mouth 


p b= distance aboral end of 
roof of mouth-intersection 


gl s with a. 


po Angle between opisthocranion +opisthion and roof of mouth 
(aborally extended). 

y Angle between basion-+opisthion and roof of mouth. 

m Angle between basion + opisthion and _parietals. 

5 x=m-—o 


Minimal width between the cristae temporales. 


H. LOOSE: Dicerorhinus hemitoechus Fale. in the Netherlands 


Leiden No. 93302 vertical view. 


Fig. 4. Leiden No. 93302 lateral view. 
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Of the 3 species Coelodonta antiquitatis is the largest. D. kirch- 
bergensis is slightly larger than D. hemitoechus. This conclusion is 
confirmed by dental evidence given by Stanscue ’41. 


Mostly absent in Coelodonta, often replaced by a protuberance. 
No such protuberance has been found in the hemitoechus skulls 
examined. According to Durrst 1926, p. 143, protuberances and 
crests indicate tangential insertion of tendon on bone, incisions 
and fossae an insertion perpendicular to the bone. 


Always identical in D. hemitoechus. In Coelodonta and D. kirch- 
bergensis irregularities in the relief of the occipitals sometimes occur. 


No significant differences between the three species. 


The angle y in D. hemitoechus was found to be considerable, even 
when no exact measurements were possible. 


The figures for y are: Leiden 93302 119 
Brit. Mus. 20013 120 

27836 = 

45205 110 

M5113 a 

40946 a 

45206 118 

Stuttgart 16938 116 


STAESCHE ’41 p. 115 had already seen the divergence of y. As he 
had only the Stuttgart skull at his disposal, he could not decide 
between a fortuitous aberration and a typical characteristic. The 
fact that in all skulls of D. hemitoechus the angle y was rather 
large, must decide for the latter. 

D. hemitoechus and Coelodonta antiquitatis were animals of the 
open country, getting their food mainly by grazing, also by 
browsing on small shrubs. In Coelodonta, when grazing, the skull 
was held in a direct line with the cervical vertebrae, all pointing 
in a straight line to the ground (ZEuNER 734). A reconstruction 
of D. hemitoechus should correspond with the white rhino, 
Ceratotherium simum Burch., holding its neck approximately 
horizontal when at rest, the head pointing to the ground. 


Because of the large values for y, m is small. This is seen even 
better with « (=m—o). ZHuNER has measured values between 
+27 and —14 for Coelodonta. In D. hemitoechus all measurements 


were negative. 


In Coelodonta z varies between 72 and 120, in D. hemitoechus 
between 34 and 60. 
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The material of the genus Dicerorhinus used in this study may not 
be very impressing in number. Yet a few conclusions are warranted. 

If we disregard the form of the occiput and the dentition, the best 
characteristic of Dicerorhinus hemitoechus Falc. is y. Also, m would be 
useful. but m-max. is already slightly higher than m-min. for Coelodonta 
or D. kirchbergensis. Both m and y reflect the exceptional angle between 
foramen magnum and rest of the skull. 
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MECHANICS 


ON THE TORSION OF A CYLINDRICAL BAR WITH SLITS 
BY 


J. B. ALBLAS }) 


(Communicated by Prof. C. J. BouwKamp at the meeting of October 29, 1960) 


1. Introduction 


The problem of the torsion of a cylindrical body, which was first 
investigated in 1855 by de Sarnt—-VENANT (ref. 1), has since attracted 
much attention. General methods for the evaluation of the torsional 
rigidity have been developed and an exact or approximate solution is 
known for a great number of cross-sections. It is remarkable that not 
much attention is paid to a class of interesting and important special 
problems, which arise if the cylindrical bar is weakened by one or more 
slits. Although some special cases of this kind have been treated, a general 
theory that bears practical results is not yet available. The energy method 
for the solution of torsion problems can of course be used here, but it 
does not suit this purpose very well. 

The solutions which have been found for the special cases are nearly 
all obtained by means of some conformal mapping. FILON (ref. 2) has 
investigated the problem of the torsion of a cylinder of circular cross- 
section with one or two keyways. The boundaries of the keyways are 
confocal hyperbolas and he gives special attention to the case when the 
hyperbolas degenerate into straight lines. SHePHERD (ref. 3) treats the 
problem of a circular shaft with one, two or four equal and evenly placed 
slits of any depth from the surface towards the axis. WIGGLESWORTH 
(ref. 4) has generalized one of the problems of Shepherd: the problem of 
the circular cross-section with two slits extending from the ends of a 
diameter to any depth along the diameter. 

The first who has attacked the class of slit problems in a different way, 
viz. by means of an integral equation technique, is Nowacki (ref. 5). 
Nowacki gives the solution for the special case of a rectangular cross- 
section. In general his integral equation is hard to solve. 

It is the purpose of this paper to give a general method for the evaluation 
of the reduction in torsional rigidity. To this end we also start from an 
integral equation formulation as Nowacki does. We indicate a solution 
method for the singular integral equation(s), by reducing it (them) to 
Fredholm equation(s). But we denote special attention to a variational 
principle, which we derive from the integral equation(s), and from which 


1) Technological University of Eindhoven. 


48 


we deduce directly an approximate value of the reduction in torsional 
rigidity. This principle is similar to RAyLErcH’s principle (ref. 6) or to 
Levine and Scuwrinaer’s principle (ref. 7) for diffraction problems. 

In this method we assume that the problem is solved (or approximately 
solved) for the same cross-section with less or no slits. Using trial functions 
we derive from this solution an approximate value of the reduction factor. 
As we know the character of the singularity at the end of the slit and 
the value at the boundary, it is easy to find “physically reasonable” trial 
functions, which can be expected in many cases to give a good approxi- 
mation without too much computation. 

As an example we solve the integral equation for a complete radial 
slit in a circular cross-section. We compare the exact solution with the 
result which is obtained by the variational method. The agreement is 
very good. 


2. Mathematical formulation of the problem 


We consider a cylindrical bar with an arbitrary cross-section with one 
or more slits (ef. fig. 1). The only non-vanishing stresses can be derived 
from a stress function F(x, y) according to 


= 
Fig. 1 
OF Ou 
(1) Tea = OG >, Tey = — 0G —, 


where G is the shear modulus and « the twist per unit length. The function 
F(x, y) satisfies 


(2) AF = 


in the domain R of the cross-section, while on the boundary C; of R 
(including the slits), the boundary condition is 


(3) F=0. 
If we have solved the problem (2), (3) we find the torsional rigidity D by 
(4) D=2G ff F(x, y) da dy. 
R 
There are several methods which can be tried to solve the problem (2), (8). 


One of them makes use of a Green’s function formulation. If we know 
the solution G(x, y; &,) of 


eG eG E 
(5) AG= sa + ya = — Indes) o(y—n), 
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with the boundary condition 
(6) G=0, along Cj, 


we find from Green’s formula 


| | (G4F— Fae) dx dy=—2[ [@ de dy + 2nF (5,4) = | (G = 


R R C 


- 


1 


from which we derive 
(7) F(é,n) = 2 | | G(x, y; &,n) dx dy, 
. 


where n is the normal in the outward direction. 
The formulation (5), (6) cannot be considered as a simplification of the 
problem (2), (3). If we replace the boundary condition (6) by 


(8) G=0, along C, 


where C is the boundary of the domain R, excluding the slits (or including 
a smaller number of slits), the problem (5), (8) is obtained, which can 
more easily be solved. In the following we shall assume that the stress 
function solution ((2), (3)) of this problem is known. We now have 
(9) F(é,n) = =| G(x,y; &,7) da dy + =| a ds, 

R L 
where L is the union of the disconnected ares of the excluded slits and s 


is measured along the slits. Introducing the unknown function S(s) by 


q SViext — oF OF 
(10) SOA ea 
where the + and — signs denote the “upper” and “lower” sides of the 


ares, we have from (9), (7), 
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1 /f l fy : “iff 
(11) F(é,7) = Fo(&,7) iis or | aa. E, ”) S (s) ds, 
L 


where F(&, 7) is the stress function of the bar without slits (or with a 
smaller number of slits). If we denote by o the pair of coordinates 
(é,7) on L, we find from (11) for S(s) the integral equation 


(12) f G(s, o) S(s) ds = —2nFo(o), (o on 0): 
L 
We remark that, in general, the formula (12) is an abbreviated notation 


for a system of integral equations. 
With the help of (4) we find for the rigidity 


(13) D=Do+G@ f Fo(s) S(s) ds, 
L 
where Do is the rigidity of the cross-section with a smaller number of 


slits. In the derivation of (13) the formula (7) has been used. For the 
reduction factor D,,, we have from (13) the simple expression 


(14) Dea = —G J Fo(s) S(s) ds. 


4 Series B 
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The problem formulated in this way is divided into two simpler ones: 
(a) the construction of the Green’s function for the contour, 
(b) the exact or approximate solution of the integral equation (12). 

The construction of the Green’s function for a simple contour can be 
reduced to the determination of a function which maps the domain D 
onto the unit circle. If, with the complex variable 
(15) 2=x+ty, 
this function is denoted by 


(16) w=f(z), 


the Green’s function becomes (ref. 8) 


vo | HHA (20) 
17 G(z, 20) = — lo re ; 
ol ii © |i Fea fle) 
where the point z=2o is carried into the origin w =0. We can always write 
(18) G(z, 20) = —log |z—zo| +9(z, 20), 


where g(z, 20) is harmonic in R. 
We also can write 


G(z, 29) = —log (z—20) + A(z, 2), 


and with the help of this expression we can derive from (12) a singular 
integral equation, which can be reduced to a Fredholm integral equation 
of the second kind. We shall not develop this method here. 


3. The variational principle 


The main purpose of our considerations is the calculation of the 
reduction in torsional rigidity. Without solving the complete system of 
integral equations (12), it is possible to give an expression for this reduction 
which can be used for practical computations. We introduce the factor A by 


(19) 42 ee 
and we have from (14) and (12) 


2 F'o(s) S(s) ds - j F(a) S(o) do 


» i 
(20) A- (1/22) f ds Fda G(s, 0) (a) S(o) 
or 

[ J Fo(s) S(s) ds]? 


A 2x) F de Fda G(s, a) Sa) So) * 


The expression (21) is more than a trivial identity, for it is easy to see 


that the value of A is stationary for small variations of the function S(s); 
in formula 


(22) oA 


I| 
= 


al 


The reduction in torsional rigidity is obtained if we can find a function S(s) 
which makes the expression (21) stationary. If we assume for S(s) a 
physically plausible expression S*(s), which contains a number of unknown 
constants, the “‘best’’ value of A is 
[ J Fo(s) S*(s) ds]? 

») — ate L 

(23) oe —(1/2z1) fds f do G(s, a) S*(o) S*(s) ° 
Ts =i 


The boundary values of S(s) are known. We have S(s)=0, if s is on O 


and S(a) ~ if a is an inner endpoint. We always have to take 


V/s—a]’ 
trial functions S*(s) satisfying the boundary conditions. 
We now prove the stationary property. We have 


pln =! ds [ do G(s, 0) S(s) S(a) = 


do G(s, a) S(s) -6 S(o) 


Fo(s) 6S(s) ds | Fo(o) S(a) do + 


L 


: \ 


[S 
Rls 
-} 
as 
< 
YH 
5) SI SE ie 
oO. 
iva) 
—— Jie | S| 


+ [ Fo(c) 6S(o) do- | Fo(s) S(s) ds. 
L 
For the right side of (24) we have 
\ = [ ds. -6S(s ee ae 2 )+ Fo(s)-2 z| 
L 


(25) 


/ sae: OS (6 Jana -G(s, a) S(s) + Fo(a) -2 a} 


By virtue of (12) the terms in braces are zero, so 6 A=0. 


4. Example. Circular cross-section with complete radial slit 


(a) exact solution 
As an example we solve the problem of the circular cross-section with 


radial slit. Aas 
The conformal mapping of the circle with radius a onto the unit circle 


|w|<1 is obtained by 
z 


(26) oO= a 


The representation, in which the point zo is carried into the origin, is 
according to (17) 


Z2—Z0 
(27) ay z a?—Zoz 4 


I 
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from which Green’s function is derived: 


2 Z 20 
5Ye) 1 SS es = 
(28) G(z, Zo) log | a ane 
With 
(29) z=re, 29= pe", 
the Green’s function becomes 
r2o2 


a? +2 279 cos (0—9) 


(30) Gir,9; 0, p) =4 log 


r2-+ 92—2ro@ cos (9—p) 


We take the slit with length 6 along the a-axis. Then, in the integral 
equation (12) we can put 6=gm=0 and we have 


: —" 
(31) | log S(r)dr=—2n-—=%, (b<g<a). 
a ai @ - 
We introduce non-dimensional quantities 
(32) F=rla, d=o/a, b=b/a, 
and equation (31) becomes (dropping the bars) 
6 
1— ° J 9.9 
(33) | log —| S(ar) dr = — “(1 - a? 07), (6<@< 1). 
0 
We proceed by transforming the variables. We introduce x, y, by 
I1—r 1—o 1—z 1l—-y 
4 r= —, Y= =, r= ;¢=——_, 
(3 ) : l+r Y l+o : lia e l+y 


and equation (33) becomes after some operations 


1—2 
S(a—— 
x+y | (« =) es m (1 a2 (=)) 1—b 7, 
T—Y | (1l+a)? pe 2a l+y/ /? (; + 4 is ) ‘ 


We define the function T(x) by 


aR es 1 yf. L—e 

(36) T(2) = qo S(a) 

and we differentiate equation (35) to y. We find 

0 

37 E ~ wde = l—y 

Me 1 J ae Poe Oy 
1+b 

If we proceed with the new variables 

(38) y2=t, P=s 

equation (37) becomes 

al 


(39) Tt Vali ee Bi oe etelese) | 
J (j/s) orig (0<t< (4) 


De 
The equation (39) is a well-known singular integral equation (ref. 9). 
We look for the solution which remains finite for s=0, but becomes 
: P 1—b 
infinite for s=( 
1+b 
of (33) is finite for r=qa and infinite for r=b. 


The solution of (39) is 
i=) 57) aay RS 
a) (EN 
i 1+b 1—j/t dt 


5) 
) . This agrees with the requirement that the solution 


(40) Ty. sine wise 
—————— — 
: ys) lpn | yt (1+ /t)® ts’ 
| (=) a) 
or 
1+» 
+t 
(41) ja eae 4a 8 eee \/ I—b\? yp _d 
CWEnea ORES Gen) ps 
— —s? 9 


We have evaluated the integral (41) for the case b=0, the complete 
radial slit. The formula becomes 


(42) Tie) [ 1-e 


and we find (dimensionally) 


2) sin A G[HAE Dew» GA) 4G 


For the torsional rigidity we find from (13) 


ie Dee | Fy(r) S(r) dr = 


0 
1 


(44) 2 = AG + 4a4G | ie) = = 


“ 


0 


= nat (1— a) 


\ 970? 


a result found already in another way in ref. (4). 
64\ . 
The numerical value of x(1 -z5) is 0.8780557. 


(b) Comparison with the approximate value found by means 
of the variational principle 


A trial function which has the right boundary value at r=a and the 


right singularity at r=0 is 
- l _- 
(45) S*(7)= (5 -7). 


With this function we compute the reduction factor A from (21). 


We have 


(46) A= = a ; i agi = 
ju) feoG9) 
0 0 
We find 
1 
(47) ne 2°) (= — 2) da = 27/20. 
0 
. l— ay ] 2 i ae 1+ y2 
cy ell) t 2 oe 


With (48) the denumerator of (46) is 


t4 


1 
S 2). S.-i ae iea 
\ 3-5 + | {——— log (1 t?) dt 4 


0 


The integrals of (49) are calculated by means of a series expansion of 
the integrand. We find 


1 


r ¢ (1—t) (1—#8) ae. _ We 18 a Eo 
(50) | —— log (1 — t?) T p dt = Pp — Ts log 2-+t 35 96° 
0 
We also find 
0-8) 7, Tae 1 ee 
(51) 4 { —OL log 7 dt=a?— 5 log 2-4-5. 


0 


From these computations we find 


: PR ee pd hid MP 146.04 
(52) A= 5% e767 0.66646 a4, 
from which we derive for the torsional rigidity 
(53) D=Do—0.67 Ga=0.90 Gas, 


The agreement with (44) is good. A better result can be obtained, if 
we add to the trial function S*(7) as given by (45) a polynomial in 7, 
containing a number of constants, which are to be determined with the 
well-known Rayleigh-Ritz technique. 


i 


On 


5. Discussion 


Although we have restricted the paper to the discussion of one simple 
example, it is obvious that the method which has been developed in 
sections 2 and 3 can be used for the calculation of the reduction in 
torsional rigidity for a wide class of cross-sections with slits. There are 
two difficulties, which arise in a more general case: 

(a) it is often rather complicated, if not impossible to find the exact 
Green’s function. 
(b) the integrations are very involved. 

In many cases of practical interest we can construct an approximate 
Green’s function with the help of (17) from a mapping function, which 
gives an approximation for the conformal representation onto the unit 
circle, The theory of this approximation is well developed (ref. 10). 
For the execution of the integrations numerical methods are preferable. 
Especially if one has the disposal of an electronic computer, the method 
given in this paper can be used to obtain rapidly accurate solutions of 
technically important problems. 
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PALEONTOLOGY 


AUTOFLUORESCENCE AND AGE OF SOME FOSSIL POLLEN 
AND SPORES 


BY 


P. VAN GIJZEL !) 


(Communicated by Prof. 1. M. vAN DER VLERK at the meeting of October 29, 1960) 


ABSTRACT 


The term fluorescence palynology is proposed by the author for the application 
of the phenomenon of autofluorescence of fossil pollen and spores in palynology. 
Some conclusions can be drawn from measurements of relative fluorescence intensity. 
Fluorescence palynology concerns firstly the relation between autofluorescence and 
the taxonomic positions of fossil pollen and spores, and secondly the relation 
between the spectral ratio of these objects in UV-light and their geological age. 
This ratio (Q-quotient) in general decreases with increasing geological age. It is 
calculated from measurements of fluorescence colour which shifts from green or 
yellow to orange with increasing age. In principle fluorescence palynology has 
solved the problem of contamination of humie sediments with secondary pollen 
and spores. 


INTRODUCTION 


Since HaiTINGER and others started to study fluorescence, this 
phenomenon is often used in biological, medical and others sciences 
(c.f. BRAuTIGAM & GRABNER, 1949; Goopwiy, 1953 and others). 

However, to this date only very little has been published on fluorescence 
of pollen grains. BERGER (1934) stated this phenomenon of recent pollen 
in his studies of hay-fever procreators. WErzeL (1959), who studied slides 
of organic sediments from Germany, was the only one who observed the 
phenomenon of luminiscence of sporomorphae and other microfossils, but 
gave no special attention to the fluorescence of pollen and spores. 

In this note special attention is given to the use of autofluorescence for 
different palynological purposes, summarized in the term fluorescence 
palynology, which is hereby introduced. This phenomenon and _ its 
applications in chronostratigraphy will be described in this preliminary 
note. 

Autofluorescence or primary fluorescence of matter is based on the 
faculty to transpose light with comparitively short wavelengths (unvisible 
UV-light) into visible light with much longer A, without any activation 
by means of colouring (which would cause secondary fluorescence). 
If fossil pollen and spores are observed with UV-light they show splendid 
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fluorescence colours, varying with age and type. The fluorescence light 
emitted by pollen and spores is not monochromatic. Their colour is 
determined by the intensities at different wavelengths. As photo 
electrical measuring instruments are too inaccurate for measuring these 
intensities, the subjective fluorescence measuring method (‘Methode der 
subjektiven Fluoreszenz-Messungen”’) according to RucH was used. 
My measurements, carried out in Prof. Ruch’s laboratory in Ziirich, 
supported some of my suppositions concerning fluorescence palynology. 

By means of this accurate method, measurements are taken of intensities, 
relative to those of fluorescing uranglass. The values of these intensities 
(in the figs. signed as I,e),) have to be considered as preliminary, because 
they are dependent on the adjustment of the apparatus. Independent of 
this adjustment however, is the spectral ratio Q, being a proportional 
number. 

The relative intensities are measured at the wavelengths of 500, 530, 
560, 600 and 640 mu. The spectral ratio is obtained by dividing the 
relative intensity at 500 my (green) by the one at 640 my (orange). The 
quotient Q forms an important scale for fluorescence colour. In the 
following the importance of these values for palynological purposes will 
be dealt with. 

While preparing the pollen preparations, the use of acid liquids which 
may diminish autofluorescence, must be avoided. 


AUTOFLUORESCENCE AND TAXONOMIC POSITION 


Some differences in colour appear when fossil pollen and spores are 
fluorescing. Firstly there is a difference in colour between pollen and spores. 
While the latter have a characteristic blue colour, the pollen is emitting 
light of longer wavelengths, colouring the grains from pale green to 
orange or red. Secondly, differences in colour are seen between various 
pollen types. The fluorescence colours of some of the most important 
pollen types from the Quaternary of N.W. Europe are given in table I. 

In this table the variation in colour has been given, which is greater 
for the grains of younger deposits (Holocene) than for older grains 
(Pleistocene). Besides, with increasing age an important change in colour 
takes place towards the red part of the spectrum (see next paragraph). 
This variation within one pollen type is related to different species inside 
that type. Investigations on fossil pollen compared with recent species 
under UV-light, are in progress. The relations between autofluorescence 
of pollen and spores and their taxonomic positions can be summarized 
as follows: 

A pollen type or growp of pollen types of a certain age may possess w 
characteristic fluorescence colour. 

Measurements of the relative fluorescence intensity give possibilities 
for the determination of the fluorescence characteristics of different pollen 
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types. In fig. | diagrams are given of some pollen and spores from Upper 
Holocene peat. They show that the spectral ratio Q is a useful criterion 
for an objective determination of fluorescence colours. The measured 
grains are chosen arbitrarily. 

To determine an average value as standard characteristic, it is necessary 
to observe a large amount of grains of the same type. Q-values of Pinus 
silvestris of Holocene age show a clear variation (see fig. 1). Notable is the 
little difference in the values of various sediment types of Holocene age. 


The cause of autofluorescence of pollen and spores is unknown. Probably 
the main chemical component of these grains, known as sporonine 
(ERDTMANN, 1954; ArzELius, 1956) or tryphine (ERDTMANN, 1960) or 
sporopollenine, is responsible for it. Possibly admixtures of other com- 
ponents in minute concentrations may cause the fluorescence. The differ- 


TABLE I 


Fluorescence-colours of some pollen and spores from Quaternary peats and clays 


Fluorescence-colours 


i Upper Holocene Middle-Pleistocene 
Type ~ —— “a rie —— 
; variation in| : variation in 
main-colour | E main-colour | wekeatet 
Pollen : 
Alnus yellowish green yellow yellow | yellowish orange 
Betula yellowish green yellow yellow _ yellowish orange 
Carpinus | yellow not observed yellow pale yellow 
Fagus yellow | from green to | —_ — 
orange 
Picea yellow | not observed _ yellow yellowish brown 
Pinus silvestris pale green yellowish green yellow | from yellowish 
green to yellow- 
| ish orange 
Quercus pale green _ yellowish orange| yellow not observed 
Corylus yellow _ yellowish green | yellowish orange yellow 
Ericaceae yellowish green | from pale green) yellow not observed 
to yellow 
Gramineae orange from yellow to) yellowish brown} not observed 
red 
Cyperaceae orange yellow orange red not observed 
Spores : 
Filices blue blueish green from pale blue} not observed 
to green 
Sphagna blue blue grey pale green not observed 


Observations in UV-light with Leitz-filters 2mm UG1 and 4mm UG1 (filtering visible light). 
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ences in fluorescence colour between pollen and spores are obvious. This 
may be caused by a difference in chemical composition. The high 


resistance of these components, however, 


makes it very difficult to 


determine the chemical composition. 


Trel. 
938 AUTOFLUORESCENCE OF 
SOME SPORES 
0.055 
COLOR: BLUE 
0.050 AGE: SUBATLANTIC 
arp EK 105/L) 
0.045 -- Filices 
Sphagnum 
0.040 
0.035 
0.030 
0.025 
0.020 
0.015 
0.010 
0.005 
LS 500 $30 560 600 640 mu 
I rel. 
0.065 0-260 AUTOFLUORESCENCE OF 
PINUS SILVESTRIS POLLEN 
0.060 a 
ni : COLOR: GREEN 
0.055 % AGE : SUBATLANTIC 
" Samples 
0050 Steet aeeeee EK 105/L (young 
4 Sphagnum peat) 
0.0451 Q_o7g 5 KG 180/BR (humic 
clay) 
0.040 
0.035 
0.030 
0.025 
a020 
0.015 
0.010 
0.005 
— 500 640 mu 


Fig. 


I rel. 

0.060 AUTOFLUORESCENCE OF 
SOME POLLEN GRAINS 

0.055 

0.050 AGE: SUBATLANTIC 
(Sample EK 105/L) 

Q0454 ee eee ee Quercus (orange) 

Ericaceae/yellow) 
0.040 Q= 1.53 ———- faqus (orange) 


D5 500 530 560 600 640mu 
I rel. 
01065 AUTOFLUORESCENCE OF 
SOME BETULACEAE POLLEN 
060 
COLOR: YELLOW 
0.055 AGE: SUBATLANTIC 
(Sample EK 105/L) 
0.050 Alnus 
-- Corylus 
0.045 ——-—-— Betula 
0.040 
0.035 
0.030 
0.025 
0.020 
0.015 
0.010 
0.005 
1— 500 530 560 600 640m 


ANALYSIS: P.v.GUZEL 
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AUTOFLUORESCENCE AND GEOLOGICAL AGE 


Fluorescence palynological investigations of many samples of varying 
age show a distinct relation between the fluorescence colour and geological 
age. This may be formulated as follows: 

Autofluorescence of fossil pollen and spores show with increasing age a 
change in colour of short wavelengths towards colours of longer wavelengths, 
followed by complete invisibility. 

This phenomenon is clearly noticeable on pollen of Pinus silvestris. 
In Holocene clays and peats these grains are fluorescing pale green, In 
similar Pleistocene deposits yellowish green to yellow and in Tertiary 
clays yellowish orange to red-brown. Other forms as Picea, Alnus, Betula 
and Corylus, from the Holocene and Tertiary show the same tendency. 
Spores of Sphagnum from Tertiary deposits are not fluorescent, those of 
Filices faintly pale blue to green. In the fluorescence of Lower Tertiary 
pollen yellow and brownish colours are dominating, while green and blue 
are lacking. Upper Mesozoic pollen and spores are even darker and less 
yellow. On the other hand, Carboniferous spores show a complete 
extinction of visible fluorescence light. The fluorescence of pollen from 
these old deposits is strikingly different from that of Holocene ones 
showing a brilliant colour-gradation. It is notable that the pollen from 
Tertiary browncoals show a revival in colour, in comparison with that 
of the Middle and Lower Pleistocene peats. This may be caused by the 
process of carbonization of peat into browncoal. This chemical reaction 
possibly activates the fluorescence of pollen. 

Fig. 2 shows the spectral composition of some Pinus silvestris grains at 
different ages. The content of green light gradually decreases in favour 
of that of orange light. The colours of the grains in this figure are as 
follows: Upper Holocene pale green (Q= 3.00), Middle Pleistocene pale 
yellow (Q = 1.75), Lower Pleistocene yellow (Q = 1.44) and Middle Oligocene 
redbrown (Q = 0.68). 

Many measurements have still to be taken to determine the exact 
relation between fluorescence and age. 


APPLICATIONS TO STRATIGRAPHICAL PROBLEMS 


1, Fluorescence palynology can be applied when palynological age- 
determination fails due to contamination of humic sediments with older, 
secondary pollen. This, for example, is the case with so-called ‘“potklei” 
(pottery clay) ofthe N. Netherlands (c.f. v. GigzEL, OVERWEEL & VEENSTRA. 
1959). The method of IvpRsEN (1936) could not be used here. The same 
problem arises in investigations of the detailed biostratigraphy of that 
territory, still in progress. In principle the application of fluorescence 
palynology has solved this problem. Of the majority of the pollen it can 
be established which grains are secondary or autochthonous. Some time 
ago I presumed (c.f. v. GiuzEL c.s., 1959) the top of this pottery-clay was 
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deposited during the Saale-glacial, while strong erosion of older sediments 
resulted in a contamination of this clay with large quantities of Pliocene, 
Miocene and older pollen and spores. Measurements of Pinus silvestris’ 
spectral ratio confirm this supposition. Some of the grains from the clay 
have a Q-value which corresponds with a Middle Pleistocene age. On 
the other hand about 80 % of the grains possess a lower spectral ratio, 
according to that of Tertiary clays. These observations are of importance 
for biostratigraphical studies of other contaminated sediments. The greater 
the difference in age of secondary and autochthonous pollen, the clearer 
is the difference in fluorescence and the easier they may be distinguished 
from each other. 


SPECTRAL RATIO AT 

DIFFERENT AGES IN 
Q=3.00 AUTOFLUORESCENCE OF 
PINUS SILVESTRIS POLLEN 


{ rel. 
0.030 


0.025 
——— SUBATLANTIC (KG 180/8R) 
--+------ HOXNIEN (Neede-clay) 
——-—— EARLY PLEISTOCENE 
(Belfeld- clay) 
ra oe MIDDLE OLIGOCENE 
(Septarien-clay) 


0.020 


0.015 


0.010 


0.005 


41—>_~—=300 530 560 600 640 mut 
ANALYSIS: Rv.GUZEL 
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2, It is well known that dating and correlating of some not con- 


taminated Pleistocene sediments is impossible because of the absence of 
zonings in pollen-diagrams (WoLDSTEDT, 1950, 1954, 1958). New points 
of view may be brought forward by the application of fluorescence 


palynology. 

3. Fluorescence palynology may be applied on species determination 
for biostratigraphical purposes, especially of Tertiary pollen and spores, 
which are hard to determine. The principle of this question has already 


been discussed above. 
Perhaps there are possibilities to test the relations of extinct and 


recent living species. 


Extensive investigations on principles of fluorescence palynology and 
on the aforementioned applications are in rapid progress and will be 
published later. Special attention will then be given to age determination 
and time correlation. 
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ZUSAMMENFASSUNG 


Der Verfasser hat Fluoreszenz-Mikroskopie auf die Pollenanalyse angewandt. 
Subjektive Fluoreszenz-Intensitétsmessungen haben folgendes erwiesen: 

1. Es besteht ein Zusammenhang zwischen jedem Pollentyp und_ seiner 
Fluoreszenzfarbe bei gleichem Alter. Wiahrend die Sporen blau bis blaugriin 
fluoreszieren, haben die Pollenkérner immer Fluoreszenzlicht mit héheren Wellen- 
langen. 

2. Im allgemeinen tritt in das primaire Fluoreszenzlicht der Pollen- und Sporen- 
kérner bei zanehmenden Alter eine Aenderung der Intensitat nach héheren Wellen- 
langen auf. Wahrscheinlich wird das Farbeverhiltnis (‘‘spectral ratio”) durch das 
Alter bedingt. 

Die Anwendung der Fluoreszenz-Palynologie hat im Prinzip die Lésung fiir die 
Frage der sekundir Pollen verunreinigter Sedimente gegeben. 
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PHYSICS 


SOME EXPERIMENTS ON HEAT TRANSFER AND MAGNETISM 
BELOW 1°K. I 


BY 


A. R. MIEDEMA 


(Communicated by Prof. C. J. Gorter at the meeting of October 29, 1960) 


SYNOPSIS 


The experiments which will be reported have been performed at temperatures 
below 1° K which were attained by the method of adiabatic demagnetization. 
A short description of the experimental methods will be given in chapter I. 

Experiments on a new type of paramagnetic substances, solidified solutions of 
paramagnetic ions in alcohols, are described in chapter II. At very low temperatures 
one obtains a glassy substance in which the magnetic ions are regularly distributed. 
The interactions between the magnetic moments of the ions and the electrical 
field due to the surrounding alcohol and water molecules appears to be much larger 
than usually found in paramagnetic crystals. In solutions of low concentration, 
irreversible processes occur during demagnetization which indicates that long 
relaxation times exist between magnetic moments and the heat waves. 

The problem of heat transfer at very low temperatures is considered in chapter III. 
Using alcoholic solutions as a cooling substance single crystals have been cooled. 
The temperature dependence of the heat-flow measured between 0.6 and 0.015° IX 
showed that a boundary resistance was not important. At the lowest temperatures 
the heat-flow was determined by the limited thermal conductivity of the single 
crystals. The experimental data are discussed using the theory for heat transfer by 
phonons and the different methods for indirect cooling, which have been used by 
several investigators, will be compared. It is suggested that in indirect cooling 
experiments the best results will be obtained if for the cooling salt a targe number 
of flat plates, cut from a magnetic single crystal, which are glued to copper plates 
are used. This cooling salt can absorb large amounts of heat in a relatively short time. 

The magnetic behaviour of single erystals of MnNH,-, CoNH,-, CuNH,- and 
Cuk-tutton salt is described in chapter IV. The crystals of the first two compounds 
show a very anisotropic magnetic behaviour which is antiferromagnetic in one 
direction and apparently ferromagnetic in a direction perpendicular to this. From 
the two copper tutton salts investigated CuNH,-tutton salt seems to be antiferro- 
magnetic whereas in CuK-tutton salt very high values for the susceptibility have 
been measured, which suggests the occurrence of ferromagnetism. In the latter 
salt a remarkable minimum is found in the specific heat at an entropy of } R In 2. 

In chapter V the results are compared with a molecular field model, which 
accounts for both dipolar and exchange interactions. From X-ray investigations 
one expects each magnetic ion to have exchange interaction with two types of 


magnetic neighbours. The contributions of the different neighbours are estimated 
from the magnetic data. 
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The hyperfine structure interaction in magnetically ordered CoNH,- and MnNH,- 
tutton salt has been investigated by means of the anisotropic distribution of the 
y-radiation emitted by ®°Co and ®4Mn nuclei (chapter VI). The magnitude of the 
y-ray anisotropies and the direction of the nuclear alignment axes differ from the 
values predicted from the magnetic results. 


APPARATUS AND SOME EXPERIMENTAL METHODS 


1. Haxperimental arrangement 


For the method of adiabatic demagnetization one needs a suitable 
paramagnetic substance in a cryostat with liquid helium at as low a 
temperature as possible and mounted in such a way that it can be 
magnetized in heat contact with the liquid helium (isothermally) and 
demagnetized under thermal insulation (adiabatically). The equipment 
used in the Kamerlingh Onnes Laboratory has been described extensively 
several times, ref. [1], [2] and [3]. A short description will be given here. 
The cryostat is located between the poles of an 80 kW electromagnet, 
which gives a field of 22300 oersted in a pole gap of 7.5 em and which 
is homogeneous to within 1 % over a region with a diameter of 5 cm. 

In order to obtain a low starting temperature for the demagnetization 
process the liquid helium is evaporated under reduced pressure. By means 
of a combination of an oil diffusion pump (Booster) and a large capacity 
mechanical pump the temperature of the liquid helium bath can be 
reduced to somewhat below 1° K. The heat contact between the salt 
and the liquid helium bath is made and broken with the help of exchange 
gas; the sample is mounted in a vacuum space and during isothermal 
magnetization the space is filled with helium gas at a pressure of about 
10-3 mm Hg. Thereafter the gas is pumped off with a Hg diffusion pump 
and then the field is removed. The sample being now thermally isolated 
from the surrounding helium bath the demagnetization is practically 
isentropic and the temperature of the salt falls. After demagnetization 
the sample warms up again slowly. The amount of heat which enters the 
sample can be reduced to less than 0.1 erg/s, which makes it possible 
to perform the low temperature measurements nearly isothermally, if the 
heat capacities of the cooled substances are not too small. 


2. The measuring bridge 

The magnetic properties, such as the static susceptibility yo and the 
real and imaginary parts of the oscillating susceptibility y’ and y", have 
been determined by means of an induction bridge of the Hartshorn type. 
It can be used for both ballistic and a.c. measurements. A schematic 
diagram of the bridge method can be seen in fig. 1.1. M1 represents the 
paramagnetic salt surrounded by the primary and secondary windings 
of a mutual inductance. This is connected in series with a variable 
inductance Ms, which is a part of the bridge circuit. If the bridge is used 
for ballistic measurements, the primary circuit is fed from a storage cell; 
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ZC 


Fig. 1.1 The measuring bridge. 


G generator PC primary coil cryostat SB secondary bridge 

A amplifier SC secondary coils cryostat ZC ballistic galvanometer 
Am ammeter PB primary bridge V vibration galvanometer 
8 sample 


for measurements with alternating measuring fields the voltage is taken 
from a low frequency signal generator (200 Hz), the detecting device 
being a ballistic galvanometer (Z.C.) or a vibration galvanometer, 
respectively. The vibration galvanometer is used as a null detector. The 
secondary voltage of M7; can be resolved into two components only one 
of which is in phase with the alternating current and can be compensated 
by M.. The other component, which is due to a.c. losses occurring in 
the sample can be compensated by tapping off a small voltage from a 
potentiometer (71, 72, R) in the primary circuit. Zero deflection of the 
vibration galvanometer can only be obtained if both M2 and Rerr, which 
equals 717r2/(71-+7r2+ R) are set to the correct value. The two components 
of the susceptibility are found from the experimental values of Mz and 
Rett by means of 


(1.1) y' =(—Mo+Me)/¢ and (1.2) 7" = Rets/@q 


where yo = x’ —iz", My is the value of Mz when compensating the empty 
mutual inductance M), q¢ is the efficiency factor when measuring the 
variation of M, and w the angular frequency of the measuring field. In 
our apparatus 7; =12, R is a resistance of 102, 103, 104 or 1052 and re 
is variable from 0—100Q2. The values of Mp and q are found from 
susceptibility measurements between 4 and 1° K. In this temperature 


region most of the investigated paramagnetic substances follow the 
Curie law 


(1.3) 4=0/T 
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The Curie constant C equals NJ(J + 1)g?up?/3k and its value is known 
if the value of g is known from other experiments, which is usually the 
case. The temperature values are calculated from the measured vapour 
pressure of the liquid helium bath, using the tables of Van Dijk and 
Durieux (7'55 scale). The differences between the 7'sg temperature scale 
and the 755 scale are of no practical importance; the differences between 
the two scales are smaller than 0.003° K and the corresponding error in 
our measurements is less than 0.5%. From a graph of the relation 
between M2 and 1/7 the values of Mo and q can be derived easily. 


3. Temperature measurements below 1° K 


It is easy to determine the temperatures as long as Curie’s law is valid, 
since then the susceptibility is inversely proportional to the temperature. 
In general, the quantity C/z is denoted by 7*, called the magnetic 
temperature and equal to the real temperature when the energy differences 
caused by the crystalline interactions are low compared to 7’ and thus 
Curie’s law applies. 

The magnitude of a magnetic field inside a magnetized sample differs 
from the external field due to the demagnetizing field. At high tempe- 
ratures the deviation from Curie’s law caused by the demagnetizing field 
is described by 7=C/(7'—6) with 6= — NC’ (where N is the demagnetizing 
factor and C’ the Curie constant per cm*). A second correction occurs 
because of the interactions between the magnetic ions. If the magnetic 
ions have a cubic arrangement the interaction gives a contribution to 
6 of (47/3)C’, if no exchange interaction is taken into account. For a 
sphere N =47/3 and thus the resulting value of 0=0. 

The 7'* and susceptibility values, yo, 7’ and y”, given in the following 
correspond to the values of dM/dh for a spherical sample, h being the 
measuring field. When the investigated magnetic material had not a 
spherical shape, the experimental values of dM/dh were corrected by 
means of: 


(1.4) 4=Yell + (42/38—N) yy] 


where ze is the experimental value of the susceptibility and yy the 
susceptibility per cm?. 

At the lowest attainable temperatures large deviations from Curie’s 
law may occur. In measuring the absolute temperature it will be essential 
to use the fundamental definition of temperature introduced by Kelvin: 


(1.5) dQ—T Aas 


Performing a large number of demagnetizations from various initial 
conditions it is possible to obtain the relation between entropy, calculated 
from the initial field and temperature values, and a certain quantity 
which will be called the thermometric parameter P; (for instance, 7’*). 
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When a known amount of heat is added to the sample the increase of 
entropy can be derived from the change of Py. 

During the measurements to be described here, the heat is added by an 
alternating magnetic field. The heat absorption in the sample is given by: 


(1.6) dQ/at— (Y%4)horoz" 


ho being the amplitude of the alternating field. The heat may be added 
continuously or in short periods. In the first case the absolute temperatures 
follow from: 


(1.7) T = (dQ/dt)/(dP1/dt)(a8/dP,) 


In the second case the deduced value of 7' is the averaged value during 
the heating period being found by means of the formula: 

te 
(1.8) T= J (Ypho%w x" dt/(Se—S1) 

ty 
When determining absolute temperatures, it is important to use a heat 
input which is large compared to the heat leak: at least 10 times larger. 
A correction for the heat leak cannot be accurate, since it may cause a 
inhomogeneous distribution of the temperature. 

The value of 7” depends strongly on entropy in most cases and thus 
the amplitudes of the alternating fields must vary with entropy in order 
that the required heating be obtained. If the thermometric parameter 
depends on the amplitude of the measuring field, it may not be convenient 
to use the first measuring method, an accurate relation between P; and 
entropy for each field value not always being available. If quantities 
independent of io are measurable, as, for instance, v’ and y” in the case 
of alcoholic solutions of paramagnetic ions, the first method is preferable 
to the second one. 

Another way to add heat to a demagnetized sample in a nearly homo- 
geneous way can be found in the absorption of y-rays. The y-rays cause 
a constant heat input, which can be calibrated, measuring the entropy 
variation at high temperatures, where Curie’s law is valid. 

It is clear that if the relation between entropy and temperature has 
been found, any temperature-sensitive property can be used as a 
secondary thermometer as its variation as a function of entropy can 
easily be determined. Magnetic thermometers, such as y’ and 7", have 
the disadvantage that their measured value depends on the shape of 
the sample and may vary a little for different samples of the same 
paramagnetic substance. Relatively large quantities of a salt are needed 
to get a high measuring accuracy. 

In some cases a preferable thermometer may be the anisotropy of 
gamma-rays emitted by oriented radioactive nuclei. The y-ray intensities 
emitted in a given direction are dependent on temperature. In several 
cases the relation between y-ray anisotropy and temperature has been 
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thoroughly investigated. A y-ray thermometer is shape independent in 
zero field, it does not vary for different samples of the same salt and it 
can also be used if only small amounts of salt are present. The latter 
fact means an important advantage for experiments on heat transfer 
(see chapter III). A y-ray thermometer can only be used for crystals, 
not for alcoholic solutions (chapter II) and for many nuclei it is not 
highly accurate at temperatures above 0.1° K. There are, however, ions 
which have such large h.f.s. interactions that an anisotropic y-ray 
distribution is found even at 1° K for radioactive nuclei (for example 
166Ho [4]). 
(To be continued ) 
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ON SOLIDIFIED SOLUTIONS OF PARAMAGNETIC SALTS 


The possibility of investigating solutions of paramagnetic ions below 
1° K only exists if solvents are used which do not crystallize. Some 
liquids are known which can be supercooled if the temperature is lowered 
below their normal freezing point. These liquids attain a glassy state at 
very low temperatures instead of crystallizing. Examples of such liquids 
are ethanol [1], l-propanol [2] and glycerol [3]. Pure ethanol may 
crystallize between 110 and 160° K if the cooling has not been fast enough, 
whereas with propanol it is difficult to obtain the crystalline state. Gibson 
e.a. describe 1-propanol at 125° K as having a jelly-like consistency which 
contracts at lower temperatures to a hard, transparent glass. 

Solutions of copper chloride in ethanol and propanol have been inves- 
tigated with X-rays at liquid air temperatures. The X-ray patterns showed 
that no crystals of copper chloride were present. !) 

The solubility of inorganic salts in the three above mentioned liquids 
is rather high due to the large electrical dipole moment in the liquids. 
Especially nitrates and chlorides can be dissolved in large amounts; 
sulphates are less soluble in these alcohols. The maximum obtainable 
concentrations at room temperature are given in table II.I for a few 
paramagnetic salts. They are compared with the concentrations of 
magnetic ions in crystals which have been used frequently for adiabatic 
demagnetization experiments. 


TABLE II, I 
Concentrations of paramagnetic ions in some solutions and crystals, used for 
demagnetization experiments, units mole/em? 


LS 


Solubility in 
Salt l-propanol at Crystal Concentration 
300° K 
FeCls- 6H20 17.5 3% 10 CesMga(NOs)12°24H2O 28 x 10-4 
Cr(NOsz)3:9H2O eo kOe CrK(SOq)2-12H20 37 X 10-4 
CuCle:2H20 8.5 x 10-4 CukKo(SO4)2:6H2O 55 x 10-4 


1) I am much indebted to Prof. C. H. Mac Grnuavry and her coworkers who 
performed the X-ray analysis. 
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The apparatus used during the measurements on alcoholic solutions is 
shown in fig. 2.1. The solutions are put in a glass sphere with a diameter 
of about 3 cm. This sphere is connected to a piece of CrK-alum, which 
absorbs the heat coming from the bottom of the apparatus. In preliminary 
experiments it was found that the thermal contact between the sample 
and the wall of its surrounding sphere could become very bad, due to 
the contraction of the solidified aleohol. After reducing the temperature 
of the helium bath to between 4 and 1° K the temperature of the sample 
decreased only very slowly, and after a demagnetization it could take a 
long time to get the sample back to 1° K, even when exchange gas was 
present between the glass sphere and the wall of the apparatus. In order 
to improve thermal contact, helium gas was put inside the glass sphere 
before mounting. ; 
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Big. 2.1 
Apparatus used for experiments with solidified solutions. c 
a alcoholic solution ea 
b filling tube with He gas 
@ 


¢ exchange gas 
d liquid hehum bath 
e CrK-alum guard cylinder 


1. Cr(NO3)3-9H20 in 1-propanol (7.5 x 10-4 mole/cem*) 

Susceptibility measurements were carried out with an a.c. measuring 
bridge (frequency 225 Hz). In the liquid helium temperature range 
Curie’s law was valid within the experimental accuracy. Demagnetizations 
were performed starting from a number of fields and an initial temperature 
of about 1° K. From this the relation between entropy and susceptibility 
has been derived. The results, which reproduced for several measuring 
days, are shown in fig. 2.2. The real part of the susceptibility 7’ increases 
rapidly at entropies below F In 2. The imaginary part 7” becomes important 
at entropies below 0.6 R. The entropies are calculated from the initial 
values of the magnetic fields and temperatures with g = 2 and are corrected 
for the entropy of the lattice heat waves and the entropy change due to 
the electric contributions to the energy splittings. For the first correction 
a value of S,=0.05RT? has been derived from demagnetizations from 
initial temperatures of 1 and 1.5° K. If compared to the values for crystals 
this correction is rather large, as may be seen from the corresponding 
value of the Debye constant 9 which equals 40° K. 

At entropies above R In 2 the experimental points are in good agreement 
with Hebb and Purcell’s formulae, if for the stark splitting parameter 
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Fig. 2.2 Susceptibility versus entropy for a solution of chromic nitrate in propyl 
aleohol (7.5 x 10-4 mole/em?). A curve calculated with 6/k = 0.30° K is given 
in the high entropy region (dashed line). 

O 10x x'/R A 10° x x"/R 


the value 6/k=0.30° K is assumed. The calculated relation is given in 
fig. 2.2 as a dashed line. The value found for 6 roughly equals the value 
obtained with crystalline chromium nitrate, ref. [4] and [5], and it 
differs only slightly from the values found in the chromium alums [6]. 

The relatively high values of 7” below S=0.6 R gave us the opportunity 
to measure the thermodynamic temperature in this entropy region by 
the a.c. heating method. 7” proved to be independent of the amplitude of 
the alternating field for the values between 0.08 and 8 oersted used. 
The energy was added continuously (10 erg/s) or in periods of 10 seconds 
(30 erg/s). Fig. 2.3 shows the relation between entropy and absolute 
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Fig. 2.3 The relation between entropy and absolute temperature of a solution of 
chromic ions in propyl aleohol (7.5 x 10-4 mole/em3), The dashed curve is calculated 
with 6/k = 0.30° K. 
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temperature obtained. The lowest temperature attained was 0.003° K, 
The experimental points measured at entropies below 0.8 R are in 
agreement with the S versus T' relation for entropies above R In 2, calculated 
with the electrical field splitting concluded from the S versus y’ relation 
(dashed line in fig. 2.3). 

In fig. 2.3 the temperatures are plotted on a logarithmic scale, hence 
the derivatives of the S versus T curve are proportional to the specific 
heat. It can be concluded from the figure that the specific heat has a 
pronounced minimum between 0.05 and 0.02° K. At lower temperatures 
it rises roughly in proportion to 1/7. 


2. Cr(NO3)3-9H20 in 1-propanol (2.5 x 10-4 mole/em?) 

The value for the entropy of the heat waves of Sp=0.05 RT? found 
with the former solution suggested that the lattice entropy was small 
enough to reach lower temperatures when demagnetizing less concentrated 
chromium nitrate solutions. Experimentally, however, it was found to 
be hardly possible to reach temperatures below 0.1° K. Figure 2.4 shows 
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Fig. 2.4 Warming up curves obtained with a Cr(NOz)3:9H2O solution (2.5 x 10-4 
mole/em?) after demagnetization from H = 2 x 104 Oe, 7’ = 1° K. The sample 
was kept at H = 1 x 104 Oe, 7 = 0.5° K for a time increasing from curve a to curvef. 


a 0 min. d 6 min. 
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a number of warming-up curves obtained after demagnetization from the 
same initial field and temperature ratio (H/T =2 104 Oe/deg). For the 
curves b—f the magnetic field was reduced to half its initial value at 
first and after a waiting time, varying from 1 minute for curve b to 
31 minutes for curve f, the magnetic field was removed completely. The 
figure shows that after a longer waiting period the warming-up rate is 
significantly smaller. 

The results indicate that irreversible processes occur during demagne- 
tization, which may be ascribed, for instance, to a temperature difference 
between heat waves and electron spins. Demagnetizing to 0.5° K at first 
has the result that a large part of the entropy of the heat waves is removed 
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before the complete demagnetization takes place. The relatively large 
difference existing between the curves e and f (waiting times 15 and 
31 minutes, respectively, suggests that a relaxation time between the 
electron spin system and the heat waves at 0.5° K in a field of 104 Oe is 
of the order of 10 minutes, at least. 

In earlier experiments with a more diluted chromium nitrate solution 
(110-4 mole/em3) no temperatures lower than 0.5° K could be reached 
after demagnetization (H/T =2x 104 Oe/deg). The longest time during 
which the magnetic field was kept at half its initial value was 10 minutes 
in this case. With more concentrated solutions of chromium nitrate, 
ferric chloride and copper chloride no indications for the occurrence of 
irreversible processes have been found in the warming-up curves. 


3. CrCl3-6H2O in 1-propanol (10 x 10-4 mole/em*) 
Two, freshly prepared solutions of chromium chloride of the same 
concentration have been investigated, one solution being prepared at 
room temperature, the other one being warmed to 70° C when dissolving 
the chromium salt. Figure 2.5 gives the relation between entropy and the 
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Fig. 2.5 The relation between entropy and magnetic temperature for three solutions 
of chromic ions in 1-propanol. 


© Cr(NOs3)3-9H2O cone. 7.5 x 10-4 mole/em’, prepared at 20° C 
A CrCl3-6H2O cone. 10 x 10-4 mole/em’, prepared at 70° C 
CL) CrCl3-6H20 conc. 10 xX 10-4 mole/em’, prepared at 20° C 


magnetic temperature in the high entropy region for both chromium 
chloride solutions and for the chromium nitrate solution mentioned 
before. The curves can be characterized by the value of the Stark splitting 
parameter 0, which, as a matter of fact, will be an average value for ions 
in a solution. For the solution prepared at 70° C we derived a value for 
d/k of about 0.5° K, for the solution prepared at room temperature the 
mean value of 6/k is 1.0° K. The S versus T* relation for chromium nitrate 
is in agreement with 5/k=0.30° K. 

An explanation for the different 5 values can be found in the chemical 
behaviour of chromic ions. Two complexes of chromic ions occur in 
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chromium chloride, a violet and a green one, their formulae being 
[Cr(H20)¢]®+ and [Cr(H20),4Clo]+, respectively. According to Orde [7] 
the number of each type of ions present in a solution of chromium chloride 
depends on temperature and concentration. In crystalline chromium 
chloride only the green modification is present. Upon dissolving the salt 
[Cr(H20)4Cl2]* ions will be transformed to [Cr(H20)s]8+ ions, the equili- 
brium time varying from many days at room temperature to several 
minutes at 70° C. 

When preparing a solution at room temperature it may be expected 
that mainly [Cr(H2O),Cls]* ions are present, whereas when preparing a 
chromium chloride solution at 70°C part of the green modification will 
be transformed into the violet one. In the chromium nitrate solutions 
which were always prepared at room temperature, only [Cr(H2O)]8+ ions 
occur. The magnetic measurements with this solution thus represent the 
magnetic behaviour of chromic ions surrounded by six water molecules 
and it is not surprising that the derived value for the crystal field splitting 
(0/k=0.3° K) equals the value found for chromium nitrate crystals. 

If the results obtained with the chloride solution prepared at room 
temperature (6/k=1° K) correspond to the 6 value of [Cr(HzO)4Cly|+ ions, 
one might evaluate from the 6/k value of the warmly prepared solution 
its composition; 6/k=0.5° K corresponds to about 70 % [Cr(H2O)6¢]#+ 
and 30 % [Cr(H2O)4Cle]* ions. 


4. FeCl3-6H20 in 1-propanol (17.5 x 10-4 and 3.5 x 10-4 mole/em?) 

Two solutions of ferric chloride were investigated, the one with the 
lower concentration being prepared by diluting the concentrated solution 
a factor 5 with a propanol-water mixture of the same composition as 
present in the 17.5 10-4 mole/cm? solution. 

Susceptibility measurements showed that in the temperature range 
between 4 and 1°K Curie’s law was not obeyed. The x versus 1/T' plot 
was not a straight line but the slope of the curve decreased about 15 % 
between 4 and 1° K. 

From measurements performed with the same apparatus on a spherical 
sample of CrK-alum with the same volume and nearly the same suscepti- 
bility per cm? as the 17.5 x 10-4 mole/cm® ferric chloride solution, it was 
derived that at 2° K the x value of the ferric ions was not more than 
66 % of the value calculated from the theoretical Curie constant. This 
suggests that the electrical crystal field at the position of a ferric ion 
due to the surrounding water and alcohol molecules, causes splittings of 
the S=5/z ground level which are of the order of several degrees Kelvin, 
These splittings are much larger than those known for hydrated crystals 
containing ferric ions. 

Figure 2.6 shows the relation between both 7 and 7'* and entropy. 
The 7'* values are calculated with the mean slope of the x versus 1/7’ curve 
between 4 and 1° K and the entropy values, which are calculated from 
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Fig. 2.6 The relation between entropy and 7, 7* obtained with solutions of 
FeCls-6H2O in 1-propanol. 

oO 2 scone, 1 10-4 mole/em?® 
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7.5 X 10-4 mole/cm$ Ad FP eons: 3 ~ 10-4 mole/em? 
S=5/y and g=2, have been corrected for the entropy of the heat waves, 
as found with solution 1. The 7’ values have been measured with the 
relaxation heating method. It is clear that the 7’, 7’* and entropy values 
given in fig. 2.6 have no absolute meaning; they can be used, however, 
when comparing the two solutions with different iron concentrations. 
The figure shows that in the low entropy region the values of 7’ and 7'* 
derived for the two solutions differ by practically the same factor as the 
concentrations do (a factor 5). 

The 7 and 7’* values which are reached at entropies smaller than 
R In 2 can be expected to be determined by the magnetic dipolar inter- 
actions between different magnetic ions. VAN ViEck [8] has shown 
that to a first approximation this interaction is characterized by the 
so-called internal magnetic field H;, which can be found from the formula 


(2.1) (Hi)? =2u? > v5.8 
D+q 

#? is the average square of the magnetic moment of one ion, q indicates 
an arbitrary ion and the summation over p has to be extended over all 
other ions in the whole crystal (solution). Our experimental results are 
in agreement with this formula since in a 5 times less concentrated 
solution the distances rpg are altered by a factor 3)/5 and hence H, 
decreases a factor 5. 


The result differs from the relation between H; and concentration 
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measured in crystals diluted by replacing magnetic ions by diamagnetic 
ions. In this case, (see, for instance, ref. [9] for measurements with cobalt 
ammonium tutton salt) H;2 and not H; depends linearly on concentration, 
in agreement with formula (2.1), since now the distances rp, remain 
unchanged but the number of ions over which the summation should be 
extended decreases proportionally with the dilution factor. 


5. CuClz-2H2O in 1-propanol (8.5 x 10-4 mole/cm3) 

In the liquid helium temperature range the susceptibility of a copper 
chloride solution followed Curie’s law within the experimental accuracy. 
Demagnetizations have been performed from a number of fields and the 
resulting relation between S and 7,7'* can be seen in figure 2.7. The 
entropy values are calculated from the conditions before demagnetization 
with g=2.2 and S=¥% and are corrected for the entropy of the heat 
waves with S,=0.05 RT3. 
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Fig. 2.7 CuCl.-2H20 in 1-propanol (8.5 « 10-4 mole/em?). 7 and 7'* as functions 
of entropy. 


The imaginary component of the susceptibility becomes measurable 
below S=0.5 R and reaches its highest value, 7/7’ =0.03, at the lowest 
entropy. 7” proved to be independent of the amplitude of the alternating 
field between 0.08 and 0.8 oersted. Below S=0.45 & the absolute tempe- 
rature has been measured with the relaxation heating method. The 
lowest temperature reached is 0.008° K. . 

Assuming that 7’ equals 7* in the high entropy region, the specific 
heat can be derived from the slope of the curve in fig. 2.7, the temperature 
scale being logarithmic. The specific heat at relatively high temperatures 
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is found to be proportional to 1/7? and can be written as Cy = 0.0020 R/T?. 
At lower temperatures the specific heat has a minimum at an entropy 
of S=0.4 # (fig. 2.8). 
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Fig. 2.8 The specific heat of a solution of CuClz-2H2O in 1-propanol (8.5 « 10-4 
mole/cm’) as a function of temperature. 


Such a minimum in the specific heat may occur if short range order 
effects are important in the copper solution. The distances between copper 
ions will vary irregularly while the positions of the neighbouring oxygen 
atoms will be different for different copper ions. It seems reasonable to 
suppose that one particular copper ion has only one nearest copper 
neighbour, being defined as the copper ion to which it has the strongest 
magnetic interaction. An ordering of copper spins in pairs may cause 
a minimum in the specific heat if ferromagnetic or both ferromagnetic 
and antiferromagnetic interactions are present. 

A system of 2 spins (S= 14) can be described by 4 basic states. For 
completely free spins the 4 states have equal energy and the entropy 
per ion will be k In 2. The exchange interaction causes a splitting into a 
singlet and a triplet state. For interactions with the ferromagnetic sign 
the triplet state has the lowest energy and when going to low temperatures 
an entropy of k(In 4—In 3)/2 will be removed. For interactions with the 
antiferromagnetic sign the singlet lies lowest and when decreasing the 
temperature the total entropy of k ln 2 will be removed. 

Strictly, the splitting into a singlet and a triplet state is only true for 
isotropic g values. In our case the assumption of an isotropic g value 
will be roughly justified since g—2 is small. The value of g has been found 
by comparing the magnetization in a field of 2x 104 Oe with the value 
of the zero field susceptibility at 1° K. Since we are not able to measure 
the magnetization directly in such a large field it was derived from the 
magnetization curve in a longitudinal field after the sample was 
demagnetized from H/T’=2x 104 Oe/deg. At low temperatures the 
magnetization is practically saturated in a field of 600 Oe and without 
much difficulty it was possible to extrapolate the isentropic magnetization 
curve to H=2x 104 Oe. 


The zero field susceptibility at 1° K is proportional to g2 whereas the 
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magnetization in a field of 2 x 104 Oe is roughly proportional to g as it equals 
(2.2) M = Ngup(e* —e~*)/2(e* +e7*) 


where «=gupH/2kT. As a result a g value of 2.20 + 0.10 has been derived 
which nearly equals the mean value for copper ions in several paramagnetic 
crystals (chapter IV). 


6. Ce(NO3)3-9H20 in ethanol (25 x 10-4 mole/em?) 

Preliminary measurements have been performed with an alcoholic 
solution of cerium nitrate. Between 4 and 1° K Curie’s law was valid. 
After demagnetization from H/T =2x 104 Oe/deg a value for 7* of 
0.004° K was measured. In a second experiment with the same sample 
the magnetic properties were completely different. The slope of the 
calibration curve was much smaller and no low temperatures were reached 
after demagnetization. We ascribed this effect to the crystallization of 
ethanol, when the sample warmed up very slowly between the two 
measuring days. 


7. The negative result of an attempt to obtain nuclear polarization 


A solution of ferric chloride (5 x 10-4 mole/em?) in propanol, containing 
about 10 wC of both radioactive 5*?Mn and ®Co has been investigated. 
Polarizing fields between 100 and 1000 Oe were used but no y-ray 
anisotropy larger than 0.5 9 was present. The temperatures being low 
enough, the lack of y-ray anisotropy may suggest that for cobalt and 
manganese ions in an alcoholic solution the h.f.s. interaction is very 
anisotropic. This results in a strong coupling of the nuclear spins to the 
electrical crystal field, which in a solution will cause nuclear alignment 
in a direction which varies from ion to ion. For cobalt ions an anisotropic 
h.f.s. interaction exists in several paramagnetic crystals but for manganese 
ions it is rather surprising. Apparently it was not possible, even in an 
external field of 1000 Oe, to magnetize the manganese electron magnetic 
moments in a direction which differs from their preferential direction, 


as determined by the crystal field. (To be continued) 
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PHYSICS 


SOME EXPERIMENTS ON HEAT TRANSFER AND MAGNETISM 
BELOW 1°. If] 
-_ 
A. R. MIEDEMA 


(Communicated by Prof. C. J. Gorter at the meeting of October 29, 1960) 


ON HEAT TRANSFER 

1. General remarks 

One of the main purposes of demagnetization work is, besides investi- 
gating the thermal and magnetic properties of paramagnetic materials, 
to cool down other substances. In this kind of experiment it is necessary 
to achieve good thermal contact between a material which can be cooled 
by the method of adiabatic demagnetization (the cooling salt) and the 
material under investigation. In such experiments, where heat transfer 
is essential, the following thermal “‘resistances’’ may be important. 


a. The inverse heat transmission coefficient between the electron spins 
of the cooling salt and its heat waves. 

b. The inverse thermal conductivity of the cooling salt. 

c. The contact resistance at the surface between the cooling salt and 
the transfer medium. 

d. The thermal resistance of the transfer medium, which, in case a 
certain distance exists between the cooling salt and the substance 
under investigation, will be called the heat-link. 

e. The boundary resistance between the transfer medium and _ the 
substance under investigation. 

/. The inverse thermal conductivity of the substance itself. 


The cooling salt consists of a magnetically concentrated material in 
which the “resistance” a is expected to be of no importance. For example, 
according to the data of BéLtaEr [1] the energy transfer between the 
electronic spin system and lattice heat waves in a crystal of CrK-alum 
is in the liquid helium region given by @=3 x 1087'-0.9 erg/s cm3, the 
extrapolated value of this heat flow being much higher than the flow 
permitted by other resistances, as will be shown later. 

The thermal conductivity of a lattice is known to be proportional to 
its specific heat and therefore to 7 at very low temperatures. In the 
case of a non-superconducting metal, the heat conduction will be mainly 
due to the free electrons, which results in its proportionality to 7’. For 
this reason if a copper wire is used for the heat-link, its diameter being 
chosen not too small, the thermal resistance d will be less important 
than 6 at the lowest temperatures. 
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A number of physicists have performed experiments on thermal contact 
below 1° K. In order to know the cooling rate of the investigated sub- 
stance, one must measure its temperature. Thus in the first instance, 
materials having properties known as functions of temperature, are to 
be used. These materials will be referred to as the thermometer in the 
following. 

The first metal-to-cooling-salt contact was constructed by Mnnpoza [2]. 
Thin copper sheets were soldered to a copper rod. The spaces between 
the copper sheets were filled with powdered CrK-alum crystals and a 
binding agent of plastic cement. The sample was pressed under a pressure 
of 2000 atmospheres. The cooling salt and the thermometer were prepared 
in the same way, the contact area being 30 cm?. Mendoza obtained data 
on the heat transport between 1 and 0.2° K, the latter temperature 
being the equilibrium temperature when the cooling salt was demagnetized 
from H/T’'=12 103 Oe/deg, the thermometer being not demagnetized 
at all. The results for the heat flow between the two salt pills could be 
represented by the empirical formula: 


(3.1) Q=pA(T13—T8) 


A represents the contact area, 7'; is the temperature of the thermometer, 
Ts the temperature of the copper and f a constant. The experimental 
value of 6 was about 10? erg/s cm? deg?. Mendoza’s basic idea was 
developed and modified by several investigators [3], [4]. The copper 
vanes were replaced by copper wires and the contact area was largely 
increased. In an unpublished Leiden experiment by Wheatley, 500 wires 
with a total contact area of 100 cm? were used. In order to cool 10 grms 
of CrK-alum to 0.06° K with this apparatus a time of three quarters of 
an hour was needed. 

Better results were obtained by Roprnson in Oxford [5]. The sample 
was not compressed at all in his experiments, glycerol being used as a 
binding agent. With a very large contact area it was possible to cool 
a pill of cerium magnesium nitrate to 0.025° K within one hour. A cooling 
salt of this type was used in the experiments on nuclear demagnetization 
of Kurti e.a. [6]. A mixture of CrK-alum, glycerol and water cooled 
the metal nuclear sample, the total area of contact being 400 em?. In 
an experiment, where instead of the nuclear sample a CeMg-nitrate 
mixture was cooled, which was prepared in the same way as the cooling 
substance, a temperature of 0.011° K was measured 10 minutes after 
demagnetization, the heat-leak at the place of the thermometer being 
1 erg/min. The heat transport between the cerium salt and the CrK-alum 
mixture could be described by formula (3.1) when f/ is about 10° erg/s 
em? deg®. 

Experiments on the cooling of single crystals have been reported by 
WHEATLEY ¢.a. [7]. Single crystals of ferric ammonium alum were glued 
to a thermal link of copper or crystalline quartz. The experimental data 
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for the heat-flow, obtained at temperatures between 0.2 and 0.16° K 
gave much larger values than that measured in earlier experiments. 
Also, in this case formula (3.1) was in good agreement with the experiments, 
6 having a value of about 104 erg/s cm? deg®. 

All authors ascribe the thermal resistance found in the heat flow 
experiments to a contact resistance at the boundary between different 
materials. 


2. Recent Leiden experiments 

In 1956 experiments on heat transfer were started, making use of 
solidified solutions of paramagnetic ions in alcohol as a cooling salt. 
1-propanol was chosen for the solvent because of its low freezing point 
(146° K). The low freezing point is an advantage in experiments on 
thermal contact, the contact being made at a temperature below the 
normal freezing point and the shrinking of materials at lower temperatures 
being of no importance. As stated in the foregoing chapter, ethyl alcohol 
has the disadvantage that it may crystallize when passing slowly through 
the liquid-air temperature region. Nearly all nitrates and chlorides can 
be dissolved in alcohols. High values for the concentration of paramagnetic 
ions and hence large thermal capacities at low temperatures can be 
obtained with ferric chloride, chromic chloride and chromic nitrate 
(Table II, I). A solution of ferric chloride has not been used in the heat 
transfer experiments because it reacts chemically with copper. With a 
chromic nitrate solution, in combination with enamelled copper wires, 
a chemical reaction takes place only after the sample had been cooled 
down and warmed up again to room-temperature a few times. In principle 
the reactions might have been avoided if use was made of a copper 
chloride solution for the cooling salt or of gold wires instead of copper 
wires. 

In a first experiment both the cooling salt and the thermometer 
consisted of a solidified solution of chromium chloride in propanol. The 
solutions were contained in glass spheres, their volumes being 11 and 
3.5 cm3, respectively. Thermal contact was attained by means of a string 
of 1000 copper wires 0.05 mm in diameter. The area of contact between 
this heat-link and the alcoholic solutions was 30 em2 on each side, the 
average distance between the two samples being 7 em. Data about the 
heat transfer between 0.07 and 0.10° K have been obtained. 

The other experiments concern heat transport from paramagnetic 
single crystals to a solution of paramagnetic ions in 1-propanol. Fig. 3.1 
shows the apparatus. The cooling agent had a volume of 12 em? and 
consisted of a solution of chromium nitrate [Cr(NO3)3-9H20] in propyl 
alcohol, the concentration being 7.5 x 10-4 mole/em3. Thermal contact 
between the cooling solution and the copper heat-link was achieved by 
means of 130 copper wires 0.07 em in diameter and with a total area of 
70 cm? in contact with the solution. The wires were soldered in a copper 
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vase, which was connected to a wire of electrolytically pure copper with 
a diameter which was varied for several measuring days between 0.02 and 
0.09 cm. On the upper end the heat link ended in a circular copper plate 
0.02 cm thick and 1.2 cm in diameter. Two sets of experiments with single 
crystals were performed. In the first one two halves of a spherical single 
crystal of CrK-alum were glued to each side of the copper plate with 
Apiezon N grease. The crystals, which had a total volume of 0.75 cm’, 
were pressed against the copper plate by two german silver pinches 
(contact area 2 cm?). The heat flow from the CrK-alum crystals was 
measured between 0.6 and 0.04° K. The susceptibilities of the alum 
and of the cooling solution were used as thermometers. 
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Apparatus used for indirect cooling experiments. 
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CrK CrK-alum thermometer 
18 ) 1 2 
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C copper plate 
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H.L  heat-link 

C.W contact wires 
Gc F.T filling tube 
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In a second set of experiments the Crk-alum was replaced by small 
CeMg-nitrate crystals containing some radioactive ®°Co. The surface area 
in contact with the copper was about 1 cm? and the total volume of the 
crystals about 0.15 cm®. In this case the temperature was deduced from 


the y-ray anisotropy. 


3. Haperimental data and analysis 


During the first experiments on 2 chromium chloride samples a 
temperature difference was created by demagnetizing the two pills from 
different fields. The magnetic temperatures of both solutions were measured 
as functions of time with the a.c. bridge. A typical run which shows that 
it is possible to perform a temperature measurement, in less than five 
seconds, can be seen in fig. 3.2. During the first five minutes, when the 
temperature difference is large, it decreases rapidly, but it takes a long 
time before thermal equilibrium is reached. The temperatures have been 
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evaluated from the magnetic temperatures making use of gk=ta Kh 
as found for a chromium chloride solution prepared at 70° C. 

The 3.5 cm3 sample (thermometer) was demagnetized from the lower 
field and thus it has the higher temperature in fig. 3.2. During an 
experimental run it was possible to reverse the sign of the temperature 
difference by means of a magnetic field of 500 oersted, which was placed 
at the position of the smaller sample. Upon removing this field after 
several minutes, the smaller sample reached the lowest temperature, 
warming up rather quickly afterwards. 
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Fig. 3.2 Heat transfer between two alcoholic solutions of chromium chloride. 
O I/f versus time 3.5 cm? sample 
A 1/T versus time 11 cm’ sample 


Values of the heat-flow have been derived at the moments when 
relatively large temperature differences were present, the higher tempe- 
rature varying between 0.10 and 0.07° K. Assuming that the heat-flow 
could be described by: 


(3.2) Q= B(Ty®—T 2) 


where 7’, and 7’. are the temperatures of the “‘hot’’ and the cold sample, 
it was not possible to find an accurate value for n from the measurements. 
Agreement between formula and experiment could be obtained for n 
values between 2 and 3, the coefficient B being 54 x 102 erg/s deg? in the 
case where n=2 and 6.5 x 104 erg/s deg? for n=3. 

If the contact resistance between copper wires and paramagnetic 
solution dominates and if it corresponds to formula (3.2) with n=3 as 
suggested by Mendoza, Wheatley e.a. and Kurti e.a. the value of f 
equals 0.5 x 104 erg/s em? deg? since B= Af/2 and A=30 em2. In fact, 
this relatively high value for f will be still larger, because the thermal 
resistance of the copper heat-link is not negligible in this case. If 
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alternatively the total thermal resistance found is ascribed to the heat- 
link, a value for the thermal conductivity of the copper wires at 
0.1° K of K=0.4 7 watt/em deg is derived. This value is of the order of 
magnitude of the value estimated from the electrical conductivity of 
electrolytically pure copper using Wiedemann—Franz law which roughly 
equals the value for the thermal conductivity of copper reported by 
several investigators (average value 1.8 7’ watt/em deg, ref. [8]). The 
experiment with two alcoholic samples indicates that the thermal contact 
between an alcoholic solution and copper wires is good, the values of 
6 being not smaller than the values reported by Wheatley e.a. for a 
single crystal, glued very carefully to a copper plate. 

During the experiment in which CrK-alum crystals were used as a 
thermometer, the crystals were heated up by a 2 MHz oscillating field 
after demagnetization from fields between 5 and 22 kOe. According to 
AMBLER and Hupson [9] chromium alums have a strong absorption 
for frequencies of this order of magnitude up to temperatures of about 
0.1° K. At higher temperatures the heat will mainly be generated in 
the copper plate. Starting from 0.01 K° it was found to be possible to 
reach temperatures between 0.05 and 0.5° K in a few seconds. After 
each heating period the susceptibilities of the CrK-alum and the cooling 
solution were measured as functions of time by means of an a.c. measuring 
bridge. A frequency of 225 Hz and a field amplitude of about 0.1 oersted 
were used. 

From the susceptibilities the temperatures of the “hot” and the “cold” 
sample (7, and 7.) were obtained as functions of time. The relation 
between susceptibility and temperature for CrK-alum has been taken 
from measurements of BEuN e.a. [10], a correction for the value of the 
Curie-Weiss constant § found in the helium temperature range by 
Van Disk and Durtevx [11], being applied. Table III, 1 shows the 
relation between temperature, entropy and susceptibility. For the 
chromium solution this relation has been given in the foregoing chapter. 


TABLE III, I 


Thermal and magnetic properties of CrK-alum 
——_—_—_ 


EST | S/R 7/R Xx 108 
oe i Na a A 
0.40 1.335 5.23 
0.30 1.321 6.87 
0.25 (RP arat | 8.10 
0.20 1.233 10.0 
0.15 1.145 12.8 
0.10 TAO Zt | 
0.08 0.949 19.3 
0.06 0.866 21.1 
0.05 0.821 23.9 


0.04 0.774 25.9 
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A typical run is shown in fig. 3.3. 


Fig. 3.3 A typical run. Temperatures of CrK-alum crystal and cooling solution 
as functions of time. 
© Tp» temperature upper salt 


ol 


[]) Te temperature cooling agent 


Using the known specific heat of CrK-alum the heat flow can easily 
be calculated. The heat flow at various values of 7 and 7’c is given in 
table TIT, II. It varies from 580 erg/s for a 7, =0.55° K to 0.25 erg/s for 
T,=0.04° K, 

TABLE III, I 


The heat flow between a Crk-crystal and the cooling agent at several temperatures 


T Ts | Q @ Tn Tc Q @ 
OK oK erg/s calculated K OK | erg/a | sapaneie: 
erg/s | erg/s 
Pees Oe eee i Oe MOA teens | EES PH el, Bee. 
0.55 0.25 580 | 600 0.17 | 0.06 43 46 
0.55 | 0.385 | 450 440 0.10 0.06 | 6.7 7.2 
0.48 i O26 | 370 390 0.10 / 0.06 | 6.9 1.2 
0.40 0.17 320 | 321 0.10 | 0.04 10.2 | 9.0 
0.35 | 0.18 | 200 | 206 0.09 | 0.06 | 28 | 44 
O88 4 0:13 220 | 252 | 0.08 0.055 | 2.7 | 29 
0.30 0.115 230 | 177 ‘|! 0.08 | 0.04 | 3.2 3.2 
0.30 0.115 160 177 0.07 0.04 2.0 Ls 
0.29 0.10 165 174 0.06 0.049 1.0 1.0 
0.26 0.07 | 126 119 | (0.06 0.044 | 0.75 0.7 
0.26 0.07 108 119 0.054 | 0.04 0.44 0.49 
0.25 0.07 | 130 127 || «0.048 «=| (0.025 0.44 0.46 
0.21 0.075 | 69 | 80s || s0.04 | 0.02 0.29 0.28 
0.18 0.08 568 50 0.04 0.02 0.25 0.28 
O17 | O08 | sh] ag) | | 


In principle, in this experiment five thermal ‘‘resistances’’ occur. They 
are: the thermal resistance of the cooling solution, the contact resistance 
between the heat-link and aleohol, the resistance of the heat-link, the 
contact resistance between the crystal and copper and the thermal 
resistance of the CrK-alum crystal. The first two may be neglected since 
the contact resistance at the position of the cooling salt is much smaller 
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than at the position of the single crystals, as can be concluded from the 
two alcoholic pill experiments and also from an experiment in which 
the number of contact wires present in the aleohol was decreased by a 
factor of 20. It is only after absorption of considerable amounts of heat 
that the thermal resistance of the alcohol may become a limiting factor 
for the heat flow. 

Let us assume that the heat flow through each of the other 3 resistances 
can be described by a formula like formula (3.2) Q=B(Ty2—T") in 
which 7’) and 7, are now the temperatures on both sides of the resistance. 

On the basis of earlier data on thermal contact experiments one would 
expect that a resistance with n=3 is the most important. For this reason 
we analysed our measurements by comparing the values obtained for 
QO/(T'3 —T'3) at different values of 7y. If a formula with n=3 could 
describe the heat flow @/(7'n3—T7'3) would be independent of 7p. Figure 3.4 
shows that this is not at all the case. Q/(7n3—72) has a maximum for 
7, =0.2° K and thus a thermal resistance with n>3 has to be present 
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Fig. 3.4 Check of the formula Q = B(Tn3 — T.3). If the formula were valid the 
Q/(L'n3 — T°3) versus Tn curve should be a horizontal line. 


at temperatures lower than 0.2° K and a resistance with a corresponding 
n<3 is present at temperatures higher than 0.2° K. It was found to be 
possible to describe the heat flow by a combination of two thermal 
resistances, with the values n=2 and n=4, respectively. In table IIT, I 
the experimental values of Q@ are compared with values calculated by 


means of: 
(3.3) Q=12 x 104(7n4—T'p4) erg/s 
(3.4) Q = 26 x 102(7'p2—7'e2) erg/s 


In these equations 7’) represents a temperature between Th and ie. 
Within the experimental accuracy agreement has been obtained. A 
thermal resistance with n=3 may be present, but its coefficient B has 
to be larger than 5 x 104 erg/s deg’, which means that if we ascribe such 
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a term to a boundary resistance between the alum and the copper 
B>2.5 x 104 erg/s deg? cm?. 

The thermal resistance corresponding to formula (3.4) was obviously 
due to the heat-link, since its coefficient varies from 40x10? to 
2.7 102 erg/s deg? when the diameter of the copper was varied from 
0.9 to 0.2 mm. The corresponding values for the thermal conductivity 
of copper were K —0.327' and K=0.647 watt/em deg. These relatively 
low values [8] are not surprising as two solder joints were present in 
the heat-link. 

Menpvoza [2], WHEATLEY e.a. [7] and Kurri e.a. [6] conclude from 
their experiments that a boundary resistance is described by @/A= 
=f(Tr3—T.3). If the contact resistance is described by such a formula 
in our experiments, we must ascribe the resistance with n=4 to the 
CrK-alum crystals. This agrees with the fact that the coefficient of the 
thermal resistance with n=4 was a constant for four successive measuring 
days, while the crystal-glue-copper contact was renewed before each 
measuring day. 

It is possible to estimate the thermal conductivity of the CrK-alum 
crystals from equation (3.3). If the thermal conductivity follows a 7 
law, the coefficient 2; of A= A:T? will be of the order of 4 « 104 erg/s em deg4 
(see next section). This value is much smaller than that measured by 
GARRETT [12]: 4,=5.3x 105 erg/s cm deg?. 

The value for the thermal conductivity of CrK-alum given by Garrett 
corresponds to a value for the phonon mean free path / of about 2 mm, 
whereas our value for 4; gives / is about 0.2 mm. This rather small value 
of the phonon mean free path may be due to the fact that the CrK-alum 
crystals had been cooled four times from room temperature to very 
low temperatures before the experimental data were obtained. A similar 
discrepancy exists between the values of the thermal conductivity 
measured at higher temperatures by Bis [13] and by Van pen BRoEK 
e.a. [14]. Van den Broek’s values are a factor 10 smaller than those 
of Bil, but he also made measurements on a sample which had been 
cooled from room temperature to helium temperatures many times. 

In the experiments with CeMg-nitrate the crystals were heated up 
immediately after demagnetization in the same way as described for 
the CrK-alum crystals. The temperatures were derived from the normalized 
y-ray intensities, W, measured in the direction of the alignment axis 
and in the direction perpendicular to it. The relation between the 
anisotropy e, defined as [W(2/2)—W(0)]/W(a/2), and the temperature 
has been taken from data of AmBiErR e.a. [15]. Figure 3.5 shows the 
results of two identical runs taken together. After 6 minutes ¢ reached 
a value of about 0.08, which corresponds to a temperature of 0.015 
+ 0.002° K. In this experiment the contact area between the CeMg- 


nitrate crystals and the copper plate was about 1 cm? and the radio- 
active heat input was about 0.06 erg/s. 
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Fig. 3.5 y-anisotropy of ®°Co in CeMg-nitrate as a function of time. The crystals 
were heated up at t = 1’30”.e = 0.08 corresponds to 7’ = 0.015° K. 


4. Derivation of the value for the thermal conductivity of CrK-alum 
The heat flow differential equation has been solved for the CrK-crystals 
with the following approximations: 


a) The crystal has a cylindrical shape. 
b) The specific heat is constant. 
c) The thermal conductivity can be written as 1,7". 
d) The temperature of the cold side of the crystal 7'=0. 

In this case the heat flow problem becomes one-dimensional. The 
equation to be solved is: 

he. 4C, oT 

Se Be Le 
Assuming ¢>0, the solution can be found by the method of separating 
variables, as the relative temperature distribution will become independent 
of its actual (unknown) initial conditions. Hence we write 7’ in the form: 


(3.6) T'(x, t)=f(x) g(t) 


f(z) is found by solving numerically the equation: 


o2f4 WMOf 
(3.7) sa = eo fist 


The relative temperature distribution f(x) is shown by fig. 3.6 in terms 
of 74 as a function of x. The solution will also hold for values of 7) #40 
as long as 7'p1<74. It is possible to calculate from fig. 3.6 the relation 
between the heat transport through 1 cm? cross-section of the crystal 
and the average temperature. 


(3.8) Q/A =0.85 A1(T4—T'p4)/L 
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In our case approximations 6 and c will be roughly justified. The 
specific heat of CrK-alum is nearly constant in the temperature region 
between 0.1 and 0.04° K where the conductivity of the crystal limits 
the heat flow in our experiments. The 73 dependence of the thermal 
conductivity has been measured by Garrerr [12] and is confirmed by 
our formula (3.3). 


O5 


° x os 1 
i 


Fig. 3.6 Temperature distribution in a cylindrical CrK-alum crystal of length L, 
which is cooled at « = 0. The heat flow through each em® cross-section of the 
crystal is given in arbitrary units by the derivatives of the curve. 


When approximating formula (3.8) in the case of the two spherical 
crystals by taking L=0.6 cm, the result is: 


(3.9) Q=2.8 A\(T4 —T >?) 


The coefficient will increase if a correction is applied for the fact that 
the average susceptibility has been measured instead of the average 
temperature. The magnitude of this correction depends mainly on the 
actual value of the temperature of the cold side of the crystal and hence 
varies throughout the experiments. Assuming (3.9) to be valid, formula 
(3.3) gives a thermal conductivity of a CrK-alum crystal not higher 
than A=4 x 10478 erg/s em deg. 

5. Comparison with theory 


Debye’s theory gives the following formulae for the thermal conductivity 
and lattice specific heat of non-metals at very low temperatures 


(3.10) A=(1/3)leCy = A1T3 
(3.11) Cy = 0.4 kin273/13¢8 
in which 


/—phonon mean free path, 
Cy =lattice specific heat per unit volume, 
c=velocity of sound, 
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The heat flow through a unit area of a cylindrical crystal with length 
L is given by: 


(3.12) Q/A =A (Pn —T4)/4L = b4r2U( Trt —Te4)/30 RBc2L 


The maximum possible heat flow at a given value of 7) is obtained with 
(—ieand 7-0: 
Thus: 


(3.13) Q/A < k4x?T' 4/30 h3c2=1 x 1OU7,*/c? erg/s cm? 


In the case of CeMg-nitrate equation (3.13) gives @/A <2 x 1067'n4 erg/s em2 
the value of c being deduced from the lattice specific heat at 1° K [16]. 
In our experiments we approached this maximum value of Q/A rather 
closely. The heat transport of 0.06 erg/s with 7 =0.015° K corresponds 
to QY/A is about 10°74 erg/s cm?. 

The question arises why the theoretical value calculated for ideal 
thermal contact between a crystal and copper should not surpass the 
experimental value many times, two reflecting surfaces and a glue layer 
being present at the contact surface. There is a possibility to explain 
this, discussed hereafter. 

It may be expected that in the small CeMg-nitrate crystals the phonon 
mean free path is of the order of the crystal thickness (1 mm). In this 
case one may doubt whether it is justified to relate the flow of energy 
and the temperature, as derived from the observed y-anisotropy, by 
means of Debye’s theory. The temperature of the system of ®Co nuclei 
only depends on the frequency distribution in that part of the phonon 
spectrum with which the cobalt nuclear spins interact and not on that 
of other frequencies. It cannot be expected that in such small crystals, 
the thermal vibrations into which the kinetic energy of the emitted 
p particles is transformed attain a Boltzmann distribution. Since the 
kinetic energy of the / particles is extremely high if compared to the 
energy of thermal vibrations, one may expect that a relatively large 
number of high frequency phonons be generated. If these phonons leave 
the crystal immediately the energy thus transported adds itself to energy 
transported in accordance with Debye’s theory, in which the temperature 
of the cobalt nuclear spins and the part of the phonon spectrum that 
interacts with them, is used. This offers a possibility to explain the fact 
that the observed energy transport almost equals the value to be expected 
for ideal thermal contact. 

Other data on the heat flow from radioactive crystals come from 
experiments on the asymmetry of f-rays. In these experiments one uses 
radioactive layers on a directly demagnetized crystal. In experiments 
on the f-asymmetry of 58Co and ®Mn nuclei of Postma ¢.a. [17], where 
the radioactive heat input was of the order of 0.1 erg/s cm®, the layer 
warmed up quickly to a temperature which varied between 0.013 and 
0.020° K for different experiments, as has been derived from the observed 
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y-anisotropies. The energy transport was extremely large if compared to 
Debye’s theory which may suggest that also in these experiments there 
is no thermal equilibrium in the phonon system. 


6. The thermal resistence of a crystal-glue-copper contact 


It may be expected that apart from the radioactive heat input there 
is a heat leak in the CeMg-nitrate crystals of the order of 0.01 erg/s. This 
indicates that the resistance of a crystal-glue-copper contact is rather 
small. It will be shown in this section that the resistance of a glue layer 
may not be so important and that it is in principle even possible to have 
better thermal contact using a glue than in case of perfect mechanical 
contact without a glue layer. 

We shall first neglect the presence of the glue layer and consider the 
transmission coefficient for phonons at the surface between a crystal 
and copper. The transmission coefficient J’, is determined by the following 
formula given by RAayueiGH [18], phonons being in principle acoustical 
waves of extremely high frequencies. 


(3.14) Tix = 40102C1C2 cos 9; cos O2/(01¢1 COs 42+ o2¢2 Cos 4)? 


$2 is the angle with the normal on the reflecting surface at which a phonon 
with an angle of incidence 4; will be refracted. The angles 4; and 42 
are related by Snellius’ law of acoustic refraction. 

In order to know the total amount of refracted energy it will be 
necessary to integrate over all incident angles 4;. Calculations have been 
made by Lirrie [19], who described the thermal resistance at a surface 
between two different materials by 


(3.15) Q/A =5 x 1016[ I}/e;? + 21t/c1?](Tn* —TZ'c4) erg/s cm? 


I and J‘ are the transmission coefficients and c; and c; the propagation 
velocities of longitudinal and transverse phonons, respectively. The 
evaluation of J’; and /', is difficult because each wave, whether longitudinal 
or transverse, in general, breaks up into four waves, reflected transverse 
and longitudinal waves and refracted transverse and longitudinal waves 
as has been stated by Kousky [20]. Little calculated the transmission 
coefficients by numerical integration, using an electronic computer. In 
his paper the results are given in the form of parametric plots, where 
the values of 7 and J can be derived from the ratios of the densities 
and the velocities of sound in the two media. 

Using the values of 9 and ¢ as known from measurements at room 
temperature, Little derived that the reflection of phonons at a contact 
surface between a crystal (CrK-alum) and copper will correspond to a 
boundary thermal resistance of the kind given by formula (3.2) with 
n=4 and B/A=0.5 x 108 erg/s em? deg’. The values of o and c of CeMg- 
nitrate crystals are slightly different from the ones in CrK-alum; the 
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calculated value of B will be somewhat higher but it remains of the 
same order of magnitude, 106 erg/s em? deg’. 

If a glue layer is present there are two reflecting surfaces between the 
crystal and the copper plate. The heat flow from a crystal at temperature 
Tn, of density 91, and in which the velocity of sound is ¢1, through a glue 
layer of thickness L, density 92 and velocity of sound cz to a metal at 
temperature zero with quantities 93 and cg will be determined by the 
following equations: 


(3.16) (Q/A)s1= lise 1017(7"4 — 1734) Iy2/¢e;2 = Iss 10171(734 =1'3*)/¢22L = 
1 x 101" I537'34/co2_ erg/s cm? 


where 7, and 73 are the temperatures of the glue at the crystal-glue 
surface and the glue-copper surface, respectively, and / is the mean 
free path of the phonons in the glue. By eliminating 7, and 73 from 
the equations the ratio of the heat flow values using glue and without 
glue can be derived. The result is: 


(3.17) Qei/Q = [M12 23/I3] [Las + c22( Pes +1/L)P12/e12(l/L) | 


As an approximation formula (3.14) with 0;=62:=z2/2, will be used for 
the calculation of the transmission coefficient. The approximation gives 
reliable results as can be seen from the comparison of values of the trans- 
mission coefficients calculated by Little with the approximate ones. 
In figs. 3.7 and 3.8 the values of Q¢i/Q are given as a function of the 
density and the velocity of sound in the glue, using the following values 
for the other quantities: 91=2 g/cm’, 903=9 g/em?, c;=2.3 x 105 cm/s, 
¢3= 3.5 x 105 cm/s. When calculating the curves of fig. 3.7 the thermal 
resistance of the glue layer has been neglected (//L=1). Fig. 3.8 is con- 
structed using //L=0.1. Fig. 3.7 shows that in principle there is a 
possibility of increasing the transfer of heat when using a glue, if the 
glue has an acoustical impedance (gc) which does not differ much from 
the acoustical impedance of the crystal, but while it has a much lower 
velocity of sound. Comparison of the figures 3.7 and 3.8 shows that the 
thermal resistance of a glue layer (J/L<1) has no large influence if the 
velocity of sound has a low enough value. 

In our experiments with single crystals glued to a copper plate, the 
Apiezon glue can be expected to have a low ¢ value (¢<1 x 10° cm/s) 
which might explain why the thermal resistance of the glue layer is 
small, even if the mean free path of the phonons in the glue is several 
times smaller than the thickness of the glue layer (0.05 mm) at a 
temperature of 0.015° K. 

A tentative explanation for the 73 dependence of the heat flow through 
a erystal-glue-copper contact, as found at temperatures of about 0.2° K 
by WHEATLEY e.a. [7], might be given if the mean free path of a phonon 
in the glue increases with increasing wave length. Such a dependence 
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Fig. 3.7 The ratio of the heat flow from a perfect CeMg-nitrate crystal to copper, 
using a glue and without it, as a function of the density and the velocity of sound 
in the glue. The thermal resistance of the glue layer itself is neglected (l/L = 1). 


Fig. 3.8 The ratio of the energy transport from a CeMg-nitrate erystal to copper, 
using a glue and without it, as a funetion of the properties of the glue. The thermal 
resistance of the glue layer is described by l/L = 0.1. 


often has been found experimentally in the attenuation of normal acoustic 
waves in high polymers [20]. 

If the thermal resistance measured by Wheatley e.a. is only due to 
the resistance of the glue layer, we obtain 1/Z—0.02 at 0.2° K. (1 ~ 10-3mm). 


7. Conclusions and comparison of different indirect cooling methods 

The lowest temperatures obtainable by indirect cooling methods in 
which conductivity by phonons is essential, are of the order of 0.01° K. 
Formula (3.13) shows that the heat transportable at 0.01° K through 
1 cm? area is of the order of 0.01 erg/s. It would be difficult to reduce 
the heat leaks below this amount. A temperature of 0.01° K can only 
be reached in crystals in which the free path of the phonons is not much 
larger than the smallest dimension of the crystals, thus, only perfect or 
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very thin crystals can be used. In experiments with radioactive crystals 
no high activities can be used if «, - or f+ radiation is emitted. For a 
decay starting with K-capture no restriction is necessary. 

In experiments on, for example, two stage demagnetization, nuclear 
orientation in external magnetic fields or nuclear alignment in anti- 
ferromagnetics, where materials with large heat capacities are to be 
cooled, the time needed to obtain a considerable decrease in temperature 
will also be an important factor. For the cases where thermal conductivity 
by phonons is essential, the time factor can be estimated. Fig. 3.9 gives 
the temperature of a particular crystal as a function of time, the properties 
of the crystal being assumed as follows: the entropy of the crystal per 
mole can be written as 10-4R/T?, the thickness of the crystal is 1 mm, 
the molecular weight is 500 and the density is 2. The heat transport is 
assumed to be described by 10674 erg/s cm2. Fig. 3.9 shows that it takes 
3 hours in this case to reach a temperature of 0.01° K. For other values 
of the parameters the corresponding curve can be found easily, the times 
in fig. 3.9 being plotted on a logarithmic scale and a change of the 
parameters resulting in a constant shift of the whole curve. 
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Fig. 3.9 The temperature of a cylindrical crystal as a function of time, the 
properties of the crystal being assumed as follows. 
density = 2 g/cm? entropy = 10-4h/T? 
molar weight = 500 thickness = | mm 
The cooling rate is given by 10674 erg/s cm? 


The results obtained when cooling antiferromagnetics (chapter VI) 
can be compared with this calculation of the cooling time. A crystal of 
CoNH,-tutton salt of thickness 3 mm could be cooled to about 0.035° K 
in 1000s, the entropy of the nuclear spin system being given by 
age Ue eds ACE 

Surface resistances are not necessarily important at the lowest tempe- 
ratures if crystals are glued to flat plates. The thickness of the glue layer 
has to be of the order of 0.01 cm. The resistance of a heat-link with 
electron thermal conductivity can be neglected at the lowest temperatures 
if the diameter is not too small. With heat-links of crystalline quartz, 
as are used by Wheatley e¢.a., 0.01° K cannot be reached. The high 
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velocity of sound (5x 105 em/s) gives a maximum heat flow through 
a unit area of a quartz heat-link of about 3 x 1057 erg/s cm?. In practice 
this number will be reduced if the heat-link is longer than a few times 
the phonon mean free path. 

The cooling salts which have been used successfully by several inves- 
tigators can be divided into two groups. One group consists of the kind 
of cooling substances containing large numbers of copper wires in a 
paramagnetic salt and a binding agent, the latter being glycerol [6], 
Apiezon oil [21] or alcohol. The results have shown that the thermal 
contact using propanol is very good, the freezing point of the alcohol 
corresponding to a temperature below which the shrinking of materials 
is not so important. However, not only the thermal contact but also the 
heat capacity of the cooling salt present in the intimate surroundings 
of the copper wires should be taken into account. Magnetically more 
concentrated crystals can be added to an alcoholic solution to increase 
the heat capacity, but experimentally it turns out that not more than 
30 percent crystals (CrK-alum) can be allowed, the heat contact being 
reduced when using larger amounts of crystals. It appears that, when 
cooling large thermal capacities, a glycerol and CrK-alum mixture may 
be preferable to an alcoholic solution for the cooling salt, the heat capacity 
being larger, the contact resistance being nearly as good and its method 
of preparation at room temperature being handy. 

In other methods for indirect cooling use is made of single crystals. 
A thermal contact made by means of a glue and a copper plate is preferred 
above a contact made by single crystals grown around metal wires 
(ref. [22] and [23]), the copper wires shrinking more than the crystals 
when cooled below room temperature. : 

A difficulty encountered in using glued single crystals is the preparation 
of the crystals and the establishment of a large contact area. The thermal 
conductivity of crystals can be expected to be much better than that 
of mixtures, which means that if a large enough contact area can be 
attained a single crystal cooling salt can give the best results, especially 
if a large amount of heat must be absorbed by the cooling salt. Such 
apparatus has been used in the investigations on nuclear alignment 
(chapter VI), a contact area of about 60 cm? being obtained by using 
16 flat cylindrical plates 2 em? in cross-section and 3 mm thick eut from 
a large CrK-alum crystal. All crystal plates were glued to copper on 
both sides. In fact this cooling salt can be represented by one large 
eylinder 1.5 mm long with a cross-sectional area of 60 em2. The heat flow 
through such a CrK-alum cylinder calculated by means of formula (3.8) 
and the experimentally derived thermal conductivity is 12 x 106(7'y4—7',4) 
erg/s, which consequently indicates that its resistance can be neglected 
if the substance to be cooled has a contact area of only a few cm2. 


(To be continued ) 
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PHYSICS 


SOME EXPERIMENTS ON HEAT TRANSFER AND MAGNETISM 
BELOW 1°K. IV 
BY 


A. R. MIEDEMA 


(Communicated by Prof. C. J. Gorrer at the meeting of October 29, 1960) 


THE MAGNETIC BEHAVIOUR OF FOUR TUTTON SALTS 
BELOW 1° kK 


The investigated tutton salts have the general formula M+*+M2*(SOq4)o- 6H2O 
where M*++ is a divalent ion and M+ a monovalent ion. The crystals are 
monoclinic, their structure having been investigated by HorrMman [1]. 
The angle between the crystalline a and c axes was found to be 106°, the 
b axis being perpendicular to the ac plane. The dimensions of the unit 
cell in the directions of the crystalline a, b and ¢ axes are 9.2, 12.4 and 
6.2 A, respectively. 

There are two lattice positions for the magnetic ions, [0, 0,0] and 
[144, %, 0], the environment of the one in the ac plane being the mirror 
image of that of the other. The crystal field acting on the magnetic ions 
has roughly tetragonal symmetry and has the same character for the 
two lattice positions. Only the directions in space of the symmetry axes 
of the crystal field are different. Fig. 4.1 gives the positions of the 
tetragonal axes (T) and of the principal axes of magnetization (Ky, Ko, Ks). 
The parameters y and «, in the figure shown for the case of CoNH,4-tutton 
salt, are different for different tutton salts as can be seen in table IV, I. 
The data have been found by paramagnetic resonance methods. For 


c 


Ky 


RE OPVESe bK,- plane 


Fig. 4.1 Positions of magnetic susceptibility axes (K) and tetragonal field axes (T) 
with respect to the crystallographic axes in cobalt ammonium sulphate. 
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TABLE IV, I 
Parameters of the magnetic ions in four tutton salts. The values for Co and Mn ions 
are obtained on a ZnNH4-tutton salt crystal 
eee ee 


Tutton | | | 9) 

salt | © an a pe Tay eee ah i ry Ref, 
Se ee Chee et 
Cuk =| ¥ Wiee | 43 | 256 | ae |) oaa | Bemasey ec. [10] 
CuNH, | ¥% 65 | 39 | 2.45 | 2.06 2.12 | BLEANEY e.a. [10] 
CoNHa | % 137 34 6.45 3.06 3.06 VAN DEN Broek e.a. [11] 
MnNH, | 4/2 | 58 | 32 | D/k=0.038°K| g=2 | Breanny ea. [12] 


Co and Mn ions the parameters have been derived from measurements 
on a crystal which was magnetically diluted with Zn ions. Susceptibility 
measurements have shown that the parameters are roughly the same 
for concentrated CoNHy,-tutton salt [2], but for the concentrated 
manganese salt wy apparently is about 150°, as will be discussed in 
chap. VI, section 6. 

Experiments performed with single crystals of MnNHgy-, CoNH,-, 
CuNHy- and Cuk-tutton salt will be reported. The samples were used 
in the shapes in which they had been grown in order to prevent confusion 
as to the directions of the axes. The susceptibilities were corrected to 
the values for a spherical sample, although the investigated crystals 
did not have an ellipsoidal shape. Though a demagnetizing factor has 
only a well defined meaning for an ellipsoidally shaped crystal it is 
assumed that the susceptibilities can be corrected in first approximation 
by means of an average N value. This NV has been determined experimen- 
tally by measuring the susceptibility of a piece of Armco iron modelled 
to the shape of the crystal, as compared with that of an Armco iron 
sphere (see Hurtskamp e.a. [3]). The magnetic induction B in such a 
ferromagnetic sample is inversely proportional to 1/N for magnetic fields 
which are considerably smaller than the saturation value and considerably 
higher than the coercive field. The estimated N values are given in 
table IV, II. 

TAB LB ys ul 

Weights and demagnetizing factors (units 47/3) of the investigated tutton salt 

crystals 


Tutton salt | Weight g Mi N2 N3 


The ballistically measured susceptibilities have been corrected to the 
value for a spherical sample by means of: 


(4.1) £=Kel(1 + exe) 
in which ye is the experimental value of the susceptibility, ¢ = (42/3 —N )o/Mo, 
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o is the density and Mo the molecular weight. If strong relaxation effects 
occur the formula can not be used for susceptibilities measured with an 
alternating field. In that case the following formulae have been used: 


"” 


. Fs ae ee ee 
(4.2) x 1+ £7(%e? + Xe”) =f 2exe' 
and 

1 vo! te(ye'2+ ye"? 
(4.3) iy! Gee a ae 


~ TF e?(xe’2+ xe") + 2exe’ 


where ye’ and ye" are the experimental values of the susceptibility per mole. 


1. Manganese ammonium tutton salt 


Previously measurements on powdered crystals of MnNHy,-tutton salt 
have been performed by Cooker [4] and by STEENLAND e.a. [5], who 
found a pronounced maximum in the susceptibility versus entropy curve. 
Cooke measured the relation between temperature and entropy above 
the Néel temperature; Steenland e.a. did the same below that temperature. 

The present experiments were performed with a single crystal, the 
different magnetic axes being placed in the direction of the vertical 
measuring field on different measuring days. Fig. 4.2 shows the adiabatic 
susceptibilities in zero field. The measurements in the Ks direction were 
performed ballistically, while for the K; and Kg directions an a.c. bridge 
was used as the values of 7” proved to be negligible in these directions. 
At entropies above In 4 the susceptibilities are not very different for the 
three directions, the largest value being found in the K, and the lowest 
value in the Ke direction. At an entropy of about R In 4 the Kg suscepti- 
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Fig. 4.2. MnNH4-tutton salt. The susceptibility measured in zero field as a function 
of entropy for the three magnetic axes. 
A ky direction CL) Ke direction © Ks direction 
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bility increases sharply with decreasing entropy to a maximum value 
which if calculated per cm? equals the value 3/42 within the experimental 
accuracy of 2 %. At the entropy S=1.25 R the K, susceptibility passes 
through a pronounced maximum, whereas the K» susceptibility becomes 
independent of entropy. This entropy of 1.25 R will be considered as the 
entropy at the Néel temperature (Sy). At the lowest entropy reached 
(S=0.29 R), the Ki susceptibility has decreased to y=7x 10-8 R the 
Kg and K3 susceptibilities being 15 x 10-8 R and 63 x 10-8 R, respectively. 
Relaxation effects only occur in the K3 direction. The values of y’ and y" 
have been measured as functions of entropy. Both are strongly dependent 
on the amplitude and the frequency of the measuring field. Fig. 4.3 
shows the susceptibilities (yo, y’ and y”) in the Kg direction measured 
with an alternating field of 0.05 Oe amplitude at a frequency of 225 Hz. 
At relatively high entropies, above Sy, 7” is negligibly small and ,’ 
equals yo. Below Sy the value of 7’ decreases with decreasing entropy, 
whereas 7” reaches a maximum at S is about 0.8 R. The maximum value 
of y” corresponds to y"/yo=0.22. The susceptibilities in fig. 4.3 are 
corrected for the demagnetizing factor using formulae (4.2) and (4.3). 


In4 o=a 
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Fig. 4.3 MnNHa-tutton salt. yo, 77 and y” measured in the direction of the Ks 
axis as functions of entropy. The alternating field has an amplitude of 0.05 Oc 
and a frequency of 225 Hz. 


Bry Ks Xa 

No remanences have been observed in the three directions in which 
the susceptibility measurements were performed, using a measuring field 
of 1 oersted. By means of the relaxation heating method temperature 
determinations were carried out. As thermometric parameters the values 
of 7’ in the Kg direction in primary fields of 0.05 and 0.7 Oe were used. 
The amount of heat supplied was determined from the measured values 
of 7”. Using the field of 0.7 Oe amplitude for heating the sample, the 
heat supply was of the order of 1000 erg/s. The amplitude of the field 
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was chosen large enough that the 7" versus S curve did not show a 
maximum (see section 3). Thus y” decreases with rising temperature 
which keeps the temperature distribution in the salt more homogeneous. 
After a heating period of 100 seconds the field was reduced to 0.05 Oe 
and y’ was measured. Three such heating periods were performed after 
each demagnetization. The average temperature in one heating period 
could be calculated using the value of AS derived from the 0.05 Oe 
measurements. From the values of dy’/dt measured during a heating 
period absolute temperatures have also been derived with the aid of the 
relation 7’ = (dQ/dt)/(dS/dy’) (dy'/dt). The results are given in fig. 4.4 and 
in table IV, Ill. They can be compared with the relation between entropy 
and absolute temperature obtained by Cooke, which was obtained using 
y-ray heating at entropies above Sy (solid line in fig. 4.4). The agreement 
is good and connecting our experimental curve with that of Cooke, a 
value of 0.14° K for 7’y is found. Measurements of Steenland e.a. gave 
temperatures of the same order of magnitude but their measurements 
may have been less accurate, since in powdered crystals the relaxation 


TABLE IV, III 


Relation between entropy and absolute temperature in MnNH4-tutton salt 


| 
1.0 | 0.9 0.8 0.7 | 0.6 | 0.5 0.4 


S/R | 
0.117 | 0.109 | 0.102 | 0.093 | 0.086 | 0.080 | 0.074 


eS) 


i 
0.124 
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Fig. 4.4 MnNH4-tutton salt. The relation between absolute temperature and 

entropy. The solid line in the high entropy region has been derived from measure- 

ments of Cooke. The x’ versus S relations with two different field amplitudes ho 
are used as a thermometer. 


© ho = 0.05 Oc A ho = 0.7 Oc 
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50 


Fig. 4.5 MnNHgz-tutton salt. The susceptibility in the Ky direction as a function 
of a longitudinal field. 
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heating is very inhomogeneous due to the different orientations of the K3 
axes of the crystals with respect to the a.c. field. 

The susceptibilities in the K; and Kg directions have been measured 
adiabatically as functions of a longitudinal field and at several values of 
the entropy below Sy. For the K, directions and for three values of the 
entropy it is seen in fig. 4.5 that the susceptibility, as measured with the 
a.c. bridge, decreases slightly at first and reaches a minimum in fields 
of about 100 Oe. Upon increasing the field strength still further, the 
susceptibility shows a pronounced increase (for S=0.29 R it increases by 
a factor of 40 with respect to the minimum value) and reaches a maximum 
for field values which vary from 250 Oe for S=1.16 R to 460 Oe for 
S=0.29 R. The susceptibility curves show an antiferromagnetic character, 
the strong increase in fields of a few hundred oersted probably being due 
to the transition from the antiferromagnetic into the paramagnetic 
state. In the same way as described by Garrerr [2] for CoNH,-tutton 
salt the changes in temperature produced on isentropic magnetization 
have been derived for MnNH,-tutton salt. From the y,, versus H curves 
the adiabatic magnetization curves were deduced by graphical integration. 
Using the formula AT = —of#(IM/dS)ndH the temperature during isen- 
tropic magnetization has been calculated. The temperature reaches a 
minimum at a certain field strength and thereafter rises again. It will 
be assumed here, which may not be strictly correct, that the locus of 
the minima of the 7' versus H curves (fig. 4.6) corresponds to the boundary 
of the antiferromagnetic region. The dashed line connects these minima 
and the extrapolated value at zero temperature gives 540 Oe for the 
critical field in the K, direction. The Kg susceptibility was measured 
ballistically. It can be seen in fig. 4.7 that the magnetic behaviour in 
this direction is very different from that in the Ky, direction. The 
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Fig. 4.6 


The temperature variation during isentropic magnetization in antiferro- 
magnetic MnNH4-tutton salt. The dashed line represents the boundary of the 


magnetically ordered region. 
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Fig. 4.7 MnNH4,-tutton salt. The susceptibility in the Ks direction as a function 
of a longitudinal field. 
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susceptibility retains the very high value which it had at zero field up 
to fields which vary from 50 Oe for S=1.07 R to about 100 Oe for 
S=0.29 Rk. At higher field values the susceptibility drops to about 30 % 
of its initial value and becomes nearly independent of the magnetic field. 

The x3 versus H curves show a ferromagnetic character. The values of 
3/4 per cm® for the susceptibilities correspond to an infinite value if 
reduced to an infinitely long cylinder. The susceptibility of 1/N means, 
in fact, that the y value of the MnNH,-tutton salt crystal was equal to 
the susceptibility of the Armco iron sample of the same shape and 
volume which was used when determining the approximate demagnetizing 
factor. 

In order to derive a value for the ferromagnetic moment from the 
isentropic magnetization measurements, the experimental M values are 
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plotted versus H —(4/3)My in fig. 4.8 (My represents the magnetization 
moment per cm#). The values of the apparently ferromagnetic saturation 
moment, Mfe, are found at each entropy by extrapolating the linear 
parts of the magnetization curves to H—(42/3)M,=0. Mite increases 
with decreasing temperature as can be seen from fig. 4.9. The value 
of the ferromagnetic saturation moment at complete ordering of the 
spins is found by extrapolating the curve to T <7'y. This magnetization 
is given in table IV, IV together with the values of Ty, Sy and the 


TABLE IV, IV 


Comparison of the experimental results obtained with single crystals of CoNH4a- 
and MnNHa.-tutton salt 


Oar 


(degree) 


| Ty | Sy | Mee 


«x 106 
a0 | 


| 


CoNH: | 0.084 | 0.43 | 85 400 10 
MnNH, 0.14 | 125 | 94 540 16 
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Fig. 4.8 MnNHz-tutton salt. Isentropic magnetization curves for a field in the Ks 
direction. The external fields have been corrected with H = H. — (42/3)My. 
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Fig. 4.9 The experimental values of the ferromagnetic saturation moment, Mte, 
plotted versus temperature. 
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critical field in the K, direction. The extrapolated value of Mj. equals 
31% of the total saturation moment, the latter being calculated with 


S=5/. and g=2. 


2. Cobalt ammonium tutton salt 

A single crystal of CoNH,-tutton salt has been already investigated by 
Garrett [2] who measured pronounced anisotropies in the zero field 
susceptibilities. The K, susceptibility, which is the highest in the liquid 
helium temperature range, shows a pronounced maximum as a function 
of entropy. At the entropy corresponding to this maximum the Kg 
susceptibility increases sharply and reaches a value of about 3/47 per cm3. 
Garrett also determined the relation between entropy and temperature 
and derived a value of 400 Oe for the critical field in the Ky, direction. 

New experiments on a CoNHy,-tutton salt crystal were performed in 
order to make a complete comparison between CoN Hy- and MnNHy,-tutton 
salt possible. The ballistic measurements showed a constant susceptibility 
equal to 3/4 per cm? within the experimental accuracy (2 °,) in the Kg 
direction at entropies below Sy. The maximum in the yg versus S curve 
reported by Garrett, apparently is due to the fact that ys; was measured 
with a 40 Hz a.c. field (see section 3). 


80 
40 |— 4 
8X3 
10 R 
4 ee 
fe) H 150 Oe 300 


———_—_» 


Fig. 4.10 CoNH4,-tutton salt. The Ks susceptibility as a function of a longitudinal 
field. 


Oo S = 05060 4 A S = 0.262 R [1] S = 0.365 R 

Fig. 4.10 shows the susceptibility as a function of a longitudinal field 
in the Kg direction. The v3 versus H curves show the same character as 
the corresponding curves for the manganese salt. The susceptibility 
remains nearly constant at its large yo value in small fields and decreases 
sharply at a field the value of which is largest at the lowest entropy. 
In fields larger than 150 Oe yg becomes nearly constant again at a value 
of 10x 10-8 R. In the same way as described for MnNHy,-tutton salt 
values of the apparently ferromagnetic saturation magnetization have 
been derived at several entropies. Mt.increases with decreasing temperature 
as may be seen in fig. 4.9. The value of Mie at complete ordering of the 
spins is given in table IV, IV. It equals 65 °% of the saturation magneti- 
zation calculated with gg=4.40 and S= 1%, 
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3. Absorption and dispersion measurements 

The large specific heats of MnNHy- and CoNHy,-tutton salt below their 
Néel temperatures gave the opportunity of observing the relaxation 
effects occurring in the Kg direction at nearly constant temperatures. 
Using a primary field amplitude of 0.05 Oe and frequencies between 
3 and 1100 Hz!) the relaxation curves were measured at four values of 
the entropy. Fig. 4.11 shows the results for MnNH,-tutton salt. The 
scale factor of 7” is twice that for 7’. The measured y” values show a 
maximum within the range of applied frequencies for three of the 
entropies. From the shape of the curves, which is not a Debye shape, 
it may be seen that a well defined single relaxation time does not describe 


10r 7 
oO5 

x! 
Xo 

O° 0.25 

xl 

Xo 

e) 
° logy 2 4 
—_—s 


Fig. 4.11 MnNHz-tutton salt. Dispersion and absorption curves measured with a 
primary field amplitude of 0.05 Oe. 
O 8S O20 ais ies 0.83 
A S = 0.66 R vy S = 095 R 
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the phenomena satisfactorily, but an average value, Tay, may be derived. 
The values of tay at entropies of 0.95 R, 0.83 R, 0.66 R and 0.29 RF are 
0.2, 0.5, 2, and about 60 millisecond, respectively. The character of the 
relations between the entropy and both 7’ and y” measured for 225 Hz 
(fig. 4.3) are in agreement with these measurements, as the former may 
be considered as a dispersion and absorption curve measured at constant 
frequency with varying relaxation time. The relaxation time values proved 
to be strongly dependent on the amplitude of the measuring field 
(fig. 4.12). With increasing amplitude the relaxation time becomes shorter 
and the maximum in the 7” versus log v curve becomes higher, which 
indicates that the distribution of the relaxation times becomes narrower. 


1) The Hartshorn bridge is the one used by VAN DER Maren and VAN DEN 
Broek in paramagnetic relaxation experiments. The cooperation of Wike di, We a0), 
Broek during these experiments is gratefully acknowledged. 
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Fig. 4.12 MnNHg,-tutton salt. Dispersion and absorption curves measured at 
S = 0.29 R with several amplitudes of the primary field. 
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Fig. 4.138 CoNH4-tutton salt. Dispersion and absorption curves measured with a 
primary field amplitude of 0.05 Oe. 
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The derived values of tay at S=0.29 R are 60, 4 and 1.2 millisecond for 
a field amplitude of 0.05, 0.125 and 0.25 Oe, respectively. In fig. 4.12 it 
may be seen that at the high frequency side y’ approaches a constant 
value which corresponds to y'/yo ~ 0.30. 

The absorption and dispersion for CoNH,-tutton salt have been 
measured at three entropies using a measuring field amplitude of 0.05 Oe. 
The maxima in the y” versus log v curves are more pronounced (fig. 4.13) 
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than in the case of manganese and the relaxation times are better defined. 
From fig. 4.13 we derive tay=4, 1.3 and 0.6 millisecond for S=0.01 R, 
0.19 R and 0.25 R. respectively. At the lowest attainable entropy Tay has 
been measured for three values of the measuring field amplitude, (0.05 Oe; 
010"-Oc and! 0.24 Oe) giving for tay 4.0, 3.0, and 1.5 millisecond, 
respectively (fig. 4.14). 
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Fig. 4.14 CoNHz-tutton salt. Dispersion and absorption curves measured at 
S = 0.01 R& with several amplitudes of the primary field. 
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From fig. 4.13 and 4.14 it may be seen that 7’ tends to approach the 
same constant value of 7’/7o which is nearly equal to 0.15 for all measured 
curves at the high frequency side. 

At an entropy of 0.19 #& the influence of a magnetic field on the value 
of Tay has been studied. In small fields the values of tay were independent 
of the field within the measuring accuracy. At the field value where the 
ballistic susceptibility (fig. 4.10) drops, the relaxation effects vanish. The 
absorption and dispersion measurements for CoNHy,-tutton salt are 
satisfactorily described by one rather well defined relaxation time. For 
the manganese salt the absorption maximum is less pronounced and a 
broad interval of relaxation times seems to be present. The y’/yo value 
does not approach y’/yo=1 at low frequences. Combined with the 
dependence of the relaxation curves on the amplitude of the alternating 
field this suggests that some kind of a coercive field may be present. 
Such a coercive field would be of the order of a few hundredths of an 
oersted and therefore there is no disagreement with the absence of a 
remanence, measured ballistically with a measuring field of 1 Oe. 
STHENLAND e.a. [5] reported a small remanence of 0.10 + 0.05 gauss 
cm3/mole, whereas Kurrt [6] measured much higher values on a powdered 
sample of MnNHa,-tutton salt. 
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The differential susceptibility in a magnetic field, in which the ferro- 
magnetic saturation is reached, and the zero field susceptibility in 
alternating fields of about 1000 Hz are nearly equal for both salts. This 
suggests that the susceptibility in a longitudinal field in the Kg direction 
can be divided into two parts, a ferromagnetic part where the spin system 
shows long relaxation times and a quasi-paramagnetic part, where the 
relaxation effects have disappeared. 


4. Copper ammonium tutton salt 

Copper ammonium sulphate hexahydrate has been investigated at 
temperatures down to 0.05° K by Cooke e.a. [7]. The deviations from 
Curie’s law were measured and found to be described by 6=0.010° K. 

In the present investigation the adiabatic susceptibility has been 
measured in the directions of the three magnetic axes as a function of 
entropy in zero field and as a function of a longitudinal field at several 
entropies. 

The measurements in zero field (fig. 4.15) were performed with an 
alternating field of 225 Hz and 0.1 Oe, the differences from the ballistically 
measured susceptibility being negligible. The curves show a maximum 
for two magnetic axes (Ky, Ks) whereas the Ky susceptibility tends to 
approach a constant value at low entropies. The susceptibility maxima 
occur at entropies of 0.42 R and 0.36 R for the K, and the Kg curve, 
respectively. The Kj susceptibility, which is largest at high entropies 
due to the fact that for this direction the g factor has the largest value, 
becomes the smallest of the three susceptibilities at low entropies. 
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Fig. 4.15  CuNH4-tutton salt. The relation between susceptibility and entropy. 
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z” was smaller than 1 x 10-® R for all three measuring directions. The 
dependence of the adiabatically measured susceptibility on a longitudinal 
field can be seen in fig. 4.16 for the lowest entropy, S=0.2 R. The curves 
show a pronounced maximum in two cases (Ki, Ks) in fields of about 
150 Oe, the increase of the susceptibility being more than a factor two 
for the K, curve. The Ky curve shows an increase of only 7 % in the 
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Fig. 4.16 CuNHz-tutton salt. The adiabatic susceptibility as a function of a 
longitudinal magnetic field at S = 0.2 R. 
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same region. The susceptibility decreases sharply at fields between 150 
and 180 Oe for all three directions. 

Fig. 4.17 gives the dependence on entropy of the susceptibility in a 
longitudinal field for the Ky direction. The maximum in the y,, versus H 
curve becomes less pronounced at higher entropies but it is still present 
at an entropy of 0.38 R. The results obtained with CuNH4.-tutton salt 
suggest the occurrence of antiferromagnetism. Maxima in the x versus S 
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Fig. 4.17 CuNHa-tutton salt. The Ki susceptibility as a function of a longitudinal 
field for several entropies. 
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and the y versus H curves have been found for two magnetic axes Ky 
and K3. This may be explained if the antiferromagnetic ordering takes 
place in the K;Kg plane, while the alignment direction(s) do not coincide 


with one of the magnetic axes. 


5. Copper potassium tutton salt 

Previous demagnetization experiments with powdered crystals of 
copper potassium sulphate hexahydrate were performed by Garrett [8] 
and by STEENLAND e.a. [9], [5]. The zero field susceptibility showed a 
maximum as a function of entropy. At entropies above this maximum 
Garrett derived the absolute temperature versus entropy relation from 
his magnetic data. He gives a value of 0.034° K for the Curie-Weiss 
constant, the sign being positive. Steenland e.a. evaluated the 7’ versus S 
relation in the entropy region below the susceptibility maximum. The 
temperature of the Curie point was found to be about 0.04° K. 

The present investigation of a single crystal has been performed in 
the directions of the 3 magnetic axes. From fig. 4.18 it can be seen that 
the zero field susceptibility, measured with an alternating field of 225 Hz 
and 0.1 Oe, has a maximum for the K; and Kg directions. Its values are 
much larger than those found for the CuNH4-salt, the x scale of fig. 4.15 
being two times larger. The Ke susceptibility becomes nearly constant 
at entropies below 0.4 & at a value which differs not more than 5 °% from 
3/42 per cm®. The anisotropy of the susceptibility at low entropies differs 
completely from the one in the high entropy region where the Ky, 
susceptibility is largest, corresponding to the large g value in that direction. 
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Fig. 4.18 CukK-tutton salt. The relation between susceptibility and entropy. 
A Ky direction [| Ke direction © Ks direction 
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x” is measurably large below S=0.4 R. Its values are considerably 
smaller than those found in the CoNHy- and MnNHy,-tutton salts, ra 
being plotted on a tenfold scale in fig. 4.18. The values of x" in the Ke 
direction are higher than those found in the other two directions. In 
agreement with the small values of y” present in CuK-tutton salt the 
ballistically measured susceptibility differs only slightly from the y’ values 
given in fig. 4.18 (<2 %). 
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In small longitudinal fields the susceptibility keeps its high value but 
it decreases sharply to a few percent of its initial value in fields of about 
100 Oe. The curves obtained at S=0.2 R are shown in fig. 4.19. At higher 
entropies, but below S=0.35 R, the behaviour remains the same, the 
steep fall of the susceptibility occurring in lower magnetic fields. The 
remarkable fact that in small fields the highest susceptibility is reached 
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Fig. 4.19 Cuk-tutton salt. The susceptibility as a function of a longitudinal field 
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Fig. 4.20 CuK-tutton salt. The relation between absolute temperature and entropy. 


The solid line, given in the high temperature region, is derived from data of Garrett. 
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in the direction with the smallest g value can be checked easily. It is 
possible to derive the magnetization as a function of the magnetic field 
by graphical integration of the curves in fig. 4.19. In agreement with 
the anisotropy of the g values the magnetization in a field of 250 Oe 
is largest in the K; direction and smallest in the Ke direction. 

Using the 7” heating method, absolute temperatures have been derived. 
Both ' and 7”, measured in the Kg direction, were used as thermometers. 
The alternating field had an amplitude of 0.5 Oe which supplied heat 
to the sample of the order of 10 erg/s. The results can be seen in fig. 4.20, 
experimental points of Steenland e.a. being also given. The solid line 
drawn in the high temperature region represents data of Garrett. When 
connecting the two experimental curves a pronounced minimum of the 
specific heat is found at S=0.34 R, two maxima for the specific heat 
occurring near 0.05 and 0.007° K. 

The magnetic behaviour of CuK-tutton salt seems to be ferromagnetic 
in the Ke direction. The value of the ferromagnetic saturation moment 
obtained by extrapolation to 7’=0 differs less than 8 % from the total 
saturation moment, calculated with g=2.04 and S=}/9. 


(To be continued ) 
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PHYSICS 


THEORY OF RELAXATION PHENOMENA IN A 
MONATOMIC GAS. I 


By 


M. J. OFFERHAUS 


(Communicated by Prof. J. pp Borr at the meeting of October 29, 1960) 


1. Introduction and survey 1) 


In this paper we are going to describe some aspects of the development 
of the state of a gas consisting of identical, spherically symmetric 
molecules, starting from a completely arbitrary initial state. One can 
distinguish between the three stages defined by BocotuBov [1]: 

I. The stage of initial mixing, taking about the time of a collision, 
i.e. the average time during which two colliding molecules are in each 
other’s field of force: 


(1.1a) t ~ ofv 


(o=diameter of a molecule; v=(k7'/m)? ~ mean thermal velocity of the 
molecules; ~ means: of the order of magnitude of, or: approximately 
equal to). 

II. When stage I has finished, the development of the gas is henceforth 
described by the velocity-space-distribution function for a single particle, 
f(¢; r,t), which, under certain conditions, changes according to the 
Boltzmann equation (BE). This is called the kinetic stage, which is of a 
duration of the average time between successive collisions of the same 
molecule: 


(1.1b) ty Uv 


(/=mean free path). 

Ill. After stage Il, the maximum of information on the gas is 
contained in the first five moments (i.e. ¢-integrals) of / : 9 (mass density), 
co(‘“mass-averaged” velocity) and 7’ (temperature), which are functions 
of r and ¢; together we call them y!)(r, ¢) (A=1, 2, 3). The changes in ¢ 
of these y are given by the hydrodynamical equations (HDE). These, 
however, are not self-sufficient, like the BE; in order to solve them, 


1) Some abbreviations used throughout parts I-IV : 


BE Boltzmann equation; HDE hydrodynamical equations; HD hydrodynamical ; 
MB Maxwell—Boltzmann; ChE Chapman—Enskog ; LHS left hand side; 
RHS right hand side; ~ of the order-of-magnitude of, or: approximately equal to. 


Lists of symbols and references are added to part IV. 
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one has to watch the distribution function f, which by now has become 
a special solution of the BE closely linked with the moments y. This is 
the hydrodynamical (HD) stage, of a duration 


(1.1) try > dv 


(d—dimension of the vessel). In the absence of external forces, the HD 
stage ends in the equilibrium distribution of Maxwell-Boltzmann (MB). 

In a gas of moderate density, o <1 <d: the 3 stages will be well 
separated. 

It is evident that even a schematic description of stage I meets with 
ereat difficulties, and it will not concern us here. If, at the beginning of 
stage II, the BE for the single-particle distribution is valid—we shall 
not dwell on the arguments and restrictive conditions for its validity — this 
provides a starting point for a treatment of stages II as well as III. It 
appears, however, that what is called the kinetic theory of gases mostly 
concentrates on stage III and on the special solution of the BE, valid 
in this stage, developed by CHapman and ENnskoG (ChE solution; see, 
e.g. [2], passim; [3], chs. 7 and 8; [4]). 

Ch. 2: Now this paper will try to give a description of the kinetic 
stage II, during which many connections between stages II and III will 
become apparent. This is why we shall first give a survey of some important 
features of the ChE solution of the BE. 

Ch. 3: If, at the beginning of the kinetic stage, the distribution is 
independent of r:/=f/(¢;t), the kinetic process IL leads straight to the 
equilibrium (MB) distribution /, and the HD stage III is cut short. 
This special process, to be numbered IIA, first came up for discussion 
in connection with Boltzmann’s #-theorem. A solution for the linearized 
case, valid when f(¢; #=0)—/ is small, was given by WaLDMANN ([4], 
p. 367); the essential features of stage ILA for arbitrary initial f were 
demonstrated by Kac [5], who incidentally used a simplified one- 
dimensional model of the collision mechanism. We shall present a detailed 
treatment of this spatially uniform process; the deviation /—/ is 
subject to a decay of exponential type; this is a typical relaxation process. 

Ch. 4: The special role played by the MB distribution / in the spatially 
uniform kinetic process is, in a non-uniform gas, taken over by the 
normal ChE distribution, /°"; again f—/° dies out by relaxation. This 
process, to which, as far as we know, little reference has been made in 
literature previously, constitutes the general behaviour of a gas during 
the kinetic stage IT. 

The general solutions of the BE for stage II, suggested in chs. 3 and 4, 
have in common with the normal ChE solution for stage III, that for 
the distribution function series expansions are given, of which only the 
first few terms are elaborated, whereas the general character of the left 
out terms can be made plausible. 

Ch. 5: Next the local production of entropy during the kinetic stage 
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is discussed, for both the uniform and the non-uniform kinetic processes, 
as well as that of the HD stage. As far as possible, connections have been 
established with the terminology of the thermodynamics of irreversible 
processes. 

Ch. 6: During stage II, there is a development of the transport 
quantities (HD pressure tensor and heat current vector) towards the HD 
values of stage IIT. Collisions try to break down this transport continually 
throughout stages I] and III. This phenomenon, for which we suggest 
the term quasi-relaxation, was first mentioned, though with an inter- 
pretation we do not share, by KoHuer [6]. 

Ch. 7: The relaxation processes constituting the kinetic stage (chs. 3 
and 4), as well as the quasi-relaxation of the transport currents (ch. 6) 
are characterized by relaxation times which are inverse eigenvalues of 
the so-called linearized collision operator 2. WALDMANN ([4], p. 366 sqq.) 
as well as WanG CHANG and UHLENBECK [7] have given the complete 
eigenvalue theory of this 2 for the special case of a gas with Maxwellian 
intermolecular potential (7~4-repulsion), as well as some properties of Q 
for a gas with general, spherically symmetric force law. Accordingly, we 
shall present a qualitative evaluation of the various relaxation processes 
for the Maxwell gas, followed by an approximative treatment for a 
non-Maxwell intermolecular potential. The main relaxation times of 
stage II, which are of the order of the time between collisions, also play 
the role of coefficients in the standard expression for the ChE distribution 
function of stage III, thus providing one more link between stages IT 
and III. 

Ch. 8: Finally, in order to make the kinetic stage stand out against 
the following HD stage, we sum up the time scales for the duration of 
the various HD processes. As stated in (1.1c), the fastest of these takes 
a time d/v in comparison to the time //v needed for the kinetic stage. 


2. Velocity distribution in the hydrodynamical stage 

We shall first recall some important features of the ChE method for 
solving the BE (see ch. 1 for references). This equation, which gives 
the time derivative of the molecular distribution function f(¢ ; r,#) in 
a one-component dilute gas, reads as follows: 


(2.1) wf=—e- fF 4S th) = rf +0cf. 


The first two terms of the right hand side (RHS) represent changes in / 
due to streaming and to the influence of an external force (per mass unit) F: 


def of Oa 
(2.2) sf =— €-5,— Se 
the third one accounts for collisions: 


(2.3) dof SI(ff) = Sde fdz sin x [ders (Je—aa), 1) -{ fle’) er’) — fle) lea) 5. 


118 


Integrations are over the azimuth e, over the angle of deflection y and 
over the velocity ¢; of a particle colliding with one of velocity ¢. ¢’ and =i 
denote velocities-after-collision of two particles with initial ¢, ¢,; they 
depend on y, ¢ as well as on ¢, ¢;. The differential cross-section « depends 
on the absolute velocity of relative motion, g=|¢—¢1| and on 7; the 
functional relation «(g, 7) is the point where the intermolecular potential 
V enters. This is supposed to have spherical symmetry: V = V(r), r being 
the distance between the centres of two molecules. The «-notation, 
originally used in the quantum-mechanical theory, equally applies in a 
classical treatment. Most of our discussions will not distinguish between 
the two. 

There are 5 independent collision-invariant functions of ¢, characterized 
by the property 


(2.4) ue) + 4(e1)=4(€') + (41). 

Obviously (2.4) is satisfied by the following 7, indicated by a super- 
seript [0A]: 

(2.5) yl] =m; x02] —=me; 703] = dmc?. 

Integration with / yields the following moments of / : 

(2.6) yA = fde fy”, 


of which y!!] represents the zero moment of /; pl! the first moment, 
and yl], a contraction of the second. They can be expressed in terms 
of the 5 parameters: 9, ¢) and 7 (see ch. 1), and reversely: 


2 hb [3] 2)" 

(2.7) o=yl; oe ==; Sn on 5 ar: 
»kT'/m is really the moment {de /-)C2, and its use indicates that moments 
defined with the ‘peculiar’ or relative velocity C= ¢— eo are preferable 
over those containing c. 

In the following, the collective symbol wy will denote either y'7! 
(A=1, 2, 3) or e, c and 7’, 

Kquations for the change of y'’! in time are found by integrating y'°?! 
times the BE (2.1) over e: 


(2.8) dept = fe z™ dof 
+ a collision term that vanishes: 
(2.9) fde fAI(ffi)=0 — (A=1, 2, 3), 


this being a property of the J operator and of the collisional invariants 
74) for all functions f(¢). The streaming contribution (2.8) contains an 
F term (cf. eq. 2.2) which is treated with partial integration; thus we 
get the hydrodynamical equations (HDE): 


(2.10) p= —P- p++ Fp,  (4=1,2, 3) 
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in which 


(2.10a) piai+ = fde fey" 
and 
(2.10b) pUl- = fde (df/de) 7, 


It appears that the ¢-derivative of a moment y"! is expressed in terms 
of the values and r-derivatives of moments higher (P+) and lower (-) 
than y itself. Among all moments p+ we find the 5 yp! themselves, as 
well as the full second moment of / and a contraction of the third one; 
as in the case of yl], we prefer a formulation in terms of moments 
containing C, viz. the pressure tensor P and the heat current vector q, 
of which the first is ‘‘made traceless” by subtracting |x the trace from 
the diagonal elements, which is indicated by the symbol °; the remaining 
tensor is that of the HD pressure: 


(2.1a,b) P= P— str P= | de fmCC: q= | def; mO2C. 


The RHS, then, of the 5 eqs. (2.10) contain values and gradients of the 
y themselves and of P and q. Both these higher moments enter linearly. 
We turn first to the case of a gas without external forces, in which 
the initial distribution is spatially uniform. [(¢;t) now follows from a 
simplified BE: 
(2.12) def =S (ffir). 
The first 5 moments y”! of f are, of course, uniform too; they are, 
moreover, constant in time, for, on account of (2.9), integration of 7! 
times (2.12) simply results in 
(2.13) dy =0. 
Next we introduce the idea, characteristic of the ChE approach, that 
f shall depend on r and ¢ through the 5 moments yw: 
(2.14) fle; r, t)=f(e; p(r, 4), 


meaning that, out of all possible f having the values of the 5 integrals 
y4 in common, our f is a special one which is known once the » are given. 
Now the » do not change in ¢, and f has to follow the changes in y which 


are none: 


r) ‘ 
(2.15) def = ey =0. 
A 


Thus, the f-y-relationship (2.14) necessarily implies equilibrium and f 
follows from (2.12), which yields: 

(2.16) =O 1, 

the MB equilibrium distribution, containing the given values of p (1.€. @, €o 


and 7’) as parameters. 
The non-uniform case (f depending on r as well) is, in the ChE scheme, 
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based on the wider hypothesis of a distribution function depending on 
rand ¢ through the moments y and their r-derivatives of various order, 
so that f(r) is connected with the values of the y ina neighbourhood of r. 


One can express this by writing / as a functional of y: 
(2.17) fOr (¢; r,t) = f& (e¢lp(r, f)). 


Again, f has to follow the changes of y: 
218 of [A] 
(2.18) = 2 Sok Fil uy A 


Here 6f/dy! is the functional derivative of f with respect to y'”. In 
our case, one determines it by using the requirement that / has to be a 
function of the independent variables y”!, dy/dx, ..., °y/d2, ..., for which 
we will use the collective symbol D,y! (D for differentiation, « a variable 


index); so 


: of + th 
(2.19) ef = 2 eS (Daylily de(D, yl”), 


in which d; and the r-differentiation D, may change place: 
990 (a) 
(2.20) uf = >. > 3D.yy cD, a (dy), 


so that the functional derivative has the form of a differential operator 
in r-space: 

sf df 
2.21 Rs In i PR Pi 
( ) dpi ~ 2D, pia Do 
The second idea underlying the ChE theory is that of a series develop- 
ment of f: 


(2.22)  fch=fO4fM4...=fO14+4H442+...)=f(1+¢4), 


with terms /® normalized in such a way that, no matter where the 
series } / is cut off, one obtaines the correct values y(r) on which the 
method is based: 


(2.23) pit — yl; yD — 9 (yA = fdefirrylory, 


Unlike these p, higher moments do contain contributions from / other 
than pes in fact, because / will appear to be a local MB distribution, 


q and P are 0 when calculated on the base of /, the main contributions 
coming from /@): 


(2.24) gq —0; P 0, 


Now let us see how the RHS of the HDE (2:8 or 10) for dyp read if 
we substitute the terms of the / series (2.22) one by one. We start with 
/): if the distribution did not consist of more than a local MB {©, y™ 
would vary (cf. 2.8) according to: 


(2,25a) 94 pl — JdcylM, 7, 
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The elaborated form of this RHS (see 2.10) would show the moments 


y,P and q based on / alone, ie., the exact yOl—w but P20 
(cf. 2:23, 24) (Euler equations). Now add 7: the RHS of (2.8) are 
increased by 

(2.25b) de Vy = fdey Oy, f) 


. . . iS 

in the formulation of (2.10), terms with P@® and q are added on the 
RHS (Navier-Stokes equations), and so on. Each new term { gives 
an additive correction 


(2.25c) dep ae fdey1y.f 


on the Euler values 9;y™”; in all, 


co 
(2.26) Wry = Fy, yp 
s=0 
and, consequently, 
9 97 Of”) : 
(2.27) vf” = > Spi Dd a7) yt = > 0,47), 
A s $ 
if we define 
292 def Ou 
(2.28) de) uel) = > 5m Spin oF (op. 


Thus the series development (2.22) of a distribution function chosen to 
be a functional of » (2.17) leads to an analogous development (2.26) of 
the time derivative. Insert both series (2:22, 26) into the BE (2.1): 


(2.29) S Of — 0e 5 fo+ > 8) (f fr). 
£—=Ohs—0) r=0 r=0 3= 
This expression is divided into equations for separate determination of 
f,/®,... by introducing a parameter 4 as follows: 
(2.30a) YAO > (1/O)fO+FM+ f+... 
(2.30b) Y d¢6) > 9, + Hd¢M + 62,2) + 


Introducing these series into (2.29) and putting terms containing equal 
powers of 6 separately=0, we get: 


[1/67] Ces i2) 
(2.31) \ [1/0] wOFO =f O +I(f OAM) + S(f Of) 
[0°] Jf D + 9,Hf O) = 
] = def D+ I(fFOA2) + S(fOAM) +S(fOAM), ete. 


Of course, 9 really=1 and the §-expansion does not in itself imply a 
development of / in a series of decreasing terms. Still, one would like the 
treatment to converge rapidly, because the [0°]-equation is already very 
cumbersome. This means that the local MB f, which dutifully emerges 
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as “zero” approximation from the [1/62]-equation, should be the main 


contribution to f. In fact, if the next term satisfies /@ </, the total 
collision term J (ffi) is effectively approximated by 


\ ORE ea) 
(2.32) » = f(c)-fdefdy sin xfdaf(er){pM(e' + 6™(er’) )— dM (e)— $4 (e1)} 
ae stich 


which serves as a definition of the ‘‘linearized collision operator” 2 acting 
on any function of ¢ instead of 4. Using this 2, the [1/6 ]-equation can 
be written as 


(2.33) | FO=- 240 |. 

in which 

(2.34a) FO) = (1/f) (gf © — df) 

can be determined from / alone. Similarly the general [0"~? |-equation 
from (2.31) takes the form 

(2.33b) FO=240|, 


in which 


(2.34b) FP = (1/f){dsf"-? — > of alba tide = JT(fefi"-°)} 


contains f,/@,...,f@-D, 

The linear operator 2 defined by (2.32) if of the greatest importance, 
for the ChE theory as well as for the extensions we are going to suggest. 
For the moment, we need the following properties of 2, which we state 
without proof: 

a) A scalar product of functions of ¢ be defined by 


(2.35a) (F,@)=(G, F)S fde fO FG 


[If F and G are, instead of scalars, vector or tensor functions, we define 
correspondingly : 


(2.35b) (6 G) = pie Gy) ; = 22(F xp» @ 
Then 
(2.36a) (F, 2G)=(G, QF) 


or, introducing “‘square bracket” notation: [F,@]2(F, 2G): 
(2.36b) [f, Gj=(6,.7 
Q is a symmetric operator. 


[It is necessary to mention the following discrepancy between our 
square brackets and those used by Chapman and Cowling: 


(2.37) (Ff, QG) = ([F, G4] = vlF, Gh (n=0/m).] 
(our notation) (Ch. and C.) 
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b) 2 is a “non-negative” operator: 
(2.38) (f, QF) = (Ff, F]>0. 
c) The eigenvalue equation: 
(2.39) yd = gyyie" 


has a discrete spectrum of real eigenvalues , each with eigenfunctions 
74°" numbered by A; on account of the symmetry of Q (see (a)), 


(2.40) (io), yi’41) — 0 if He 


whereas different eigenfunctions with the same w can of course be made 
orthogonal. These 7, which we also suppose to be normalized: 


(2.41) Cee ses 5 


form a complete set of functions in ¢-space. 
d) On account of (b), 2 has non-negative eigenvalues only: 


(2.42) w>0; 


as in (2.38), the =—sign holds only if the corresponding eigenfunction 
is a linear combination of the 5 collisional invariants 7”. In fact, it is 
clear that 2 as defined by (2.32) annihilates the functions 7 characterized 
by (2.4): 

(2.43) QyA=0, 


and it is equally possible to prove (2.4) for any x satisfying (2.43). This 
explains the notation 7°” used for the collision-invariant functions from 
the beginning. [For use in the following theory, 7"! as defined by (2.5) 
have of course still to be normalized and made mutually orthogonal: 


(2.44) y/o = — ; (1021 = YZ Ce; 7/1081 = y= (cre - 5) . 
in which 

det ym 
(2.45) . or SY FC.) 


We now return to the set of equations (2.33). If we formally develop 
the LHS, F™, in a series of eigenfunctions of 2: 


(2.46) Fo = > oe fOrleodl yloal with frtea = (F, lol), 
o +d 


the solution 4 can, again formally, be written down: 
i ; 

i= r)[WA] ,,Lwa . (r)[mA] _ = f(r)[wA] 

(2.47) pee vi yi) with g” mal , 


In both (2.46) and (2.47), w=0 is excluded from the summation: 
> , for F™ as well as ¢™ are orthogonal to all collisional invariants: 


o>0 


(2.48) 0; 
(2.49) (©, 7°") =0. 
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Here (2.49) is just the normalization of 4” introduced in eq. (2.23). 
ee (2.48) is a solubility condition for eqs. (2.33), for if it did not hold, 
Q applied to 6 as of (2.47) could never yield #™. The proof that this 
condition is actually satisfied, runs as follows: substitute the definition 
of F®, (2.43b): 


( (FP, x a 2 fede yh fO pe) — = fdcx™d, ft- a)) = 


ri 


/ = ~ J de x1 94! r—1- a | (0) 3 5 [acy J (fe fi "= ey. 
The last sum is a collision contribution which vanishes separately, because 
of (2.9) and the more general property 

(2.51) fde 74 I (fg) +I (gf) } = 9, (A=4 2s 8) 


which again holds for all functions /,g. The other terms in the RHS 
of (2.50) cancel, as is proved as follows: using (2.28), 


SD) [0A] > ,(r—1—e fi = 1 (47) {0A} fle) 

(2.52) iS. | dex ry; este 'y sa | [ « le aie peal 
o=0 A 

in which the remaining integral=y'!d,. because of (2.23), and the 


functional derivation (2.21) means simply: 


< —_— erie 
(2.53) sari (YO) = Ox28.0' 
o (2.50) becomes 
(2.54) (F”, ; iol) =fde {oa Ye fT—Y — dlr-D mi a: 


for this vanishes according to (2.25c). 

(2.48) and (2.49) are apparently connected, for the normalization 
(2:23, 49) of f,...,f@-) is essential for the step (2.53) in the proof of 
(2.49), i.e., for the solubility of the equation determining /“. 

Two steps in the actual computation of the complete distribution 
f° are very difficult: 

a) successive computation of the PF (ef. 2.34b). Only F® is well 
known; from r=2 onwards, cumbersome collision terms like J(f@/,) 
occur. 

b) solving the eigenvalue problem of the 2 operator, in order to 
develop all #™ in series of y'°”! (eq. 2.46). This has only been accomplished 
for a Maxwellian intermolecular potential, V(r)=-const./r4. For some 
details, we refer to ch. 7. 

About the convergence of the basic series (2.22) something will be 
said in ch. 8. Essentially, we have a development in powers of J —=mean 
free path x gradient, which implies that convergence is ensured if the 
gradients of the moments (and gradients of these gradients, etc.) cause 
little relative change in w over a distance 1. 


(To be continued ) 
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3. Velocity distribution in the kinetic stage: uniform case 

We are now going to outline the general solution to the BE for a 
spatially uniform gas. References to existing literature dealing with this 
case were already given in ch. 1. 

Consider a uniform distribution f independent of r which, at t=0, 
has the completely arbitrary shape 


(3.1) (0) = fle:t=0). 


Kqs. (2.12), (2.13) give the simplified forms of the BE for f and of the 
HDE for the moments y of f: 


df JS (ffi); dey" —0. 


We now introduce a MB distribution f(c; y) based, at each ¢, on the 
moments y of the real f. This is a mathematical concept and _ its 
t-derivative is, of course, not given by a BE but by eq. (2.15): 


‘ (0) om) Wal = % 
(3.2) oj > i : 
4 


So we have a real distribution f changing in ¢, with constant moments y, 
and a time-independent f, equally with moments y; the difference 


~ ~ 


(3.3) jf fle; t) = f (e)$les t) 
is a function of ¢ and ¢ with, as its first 5 moments, 
(3.4) pl = de fy = (p, 7) =0 
and an initial value: 

(3.5) fb (0) S f(e;t=0) = f(0) —f. 
if is found by substituting (3.3) into the BE (2.12): 
(3.6) rWfO+A=ILO+f hO +h) 


or, with (2.16) and (3.2), 
(3.7) f=I(f OA) + IFA) +I (f Of) +I (FA). 
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We now proceed to solve this equation, in close analogy to the 
ChE method. Again, the total / is developed in a series of terms, of which 
f® is the: first: 


(3.8) f=fO+fHfO4fO+fo+...=f/O1 44M +g +...) =fO(1+49), 
with the normalization: 
(3.9) pao yd; paw — ...=0 (pi ) = [de fir yy, 


in accordance with (3.4). Introduce the expansion (3.8) into (3.7): 


(3.10) ¥ 2: f= > FFON) (FO) =f) 


r=] r=0,1,... ¢=0,1,. 


and insert a constant =1 as follows: 


) + (1/0)f +fM + Of@+...; 


“HR? 


oo 
(3.1 1a) bs 
te 7 - , ¥ 
(3.11b) > d¢ f — (1/6) def +d¢ f@ + Ode f@ + 
or, using (3.2): 
(3.11¢) SY 46) +2 dM + 1 GS + 621 AO +... 
r=1 


Putting these series in (3.10) and separating terms according to powers 
of 0, we find: 


\ [1/02] 0=JS(f f,) 
(3.12) [1/0] df Y= fF OPH —J(f OFM) + I(F Mf, ) 
/ [0°] df 2 =f OPA) =S(fOFD) + IFOAM) +I(FOAM), ete. 


The choice of the 0-developments (3.11) is such that the [1/62]-equation 
yields {, as required; using the definition of 2 (2.32), we write the 
[1/6 ]-equation for f™ as follows: 


(3.13) | YOY = — OGW |. 


The [0"*]-equation expresses fj in terms of /, 7, a fe); 


(3.13b) wd” = FO OPM |, (r>2) 


in which 
(3.14) Fr) = (1/f ) {7( ja fit- 1)) JF, a Se J(feDf, (1), 


Kgs. (3.13) form a hierarchy of equations for successive terms of the 
general distribution / in our special, r-independent case, just like eqs. 
(2.33) did for a special f in the general, r-dependent case, This time, 
f enters explicitly and accordingly we have to take the initial condition 
(3.5) into account. The initial deviation (0) can be spread at will over 
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the terms 4; it is logical, and will certainly promote rapid convergence, 
to assign all of it to po: 


(3.15) $%(¢;1=0)=4(0); d@(c; #=0)=...=0. 


4 * 2 « « eee: : Poy Pe 
Solution of eqs. (3.13) again involves developments of # and ¢@ in terms 
of 74, from which 7 will be absent: 


(3.16) (FO, vio) — 9. | 
(3.17) (A, v4) — 0: > 2=1,2,3. 
(3.18) ig (0),7°7) = 0. 


(3.16) is a solubility condition which, like (2.48), is easily proved to be 
satisfied with the use of (2.9) and (2.51); (3.17) is the normalizing condition 
put in (3.9), and (3.18) is simply (3.4) at ¢=0. 

As in the ChE case, two difficulties bar the way towards a complete 
solution of the hierarchy of equations (3.13): solution of the Q-eigenvalue 
problem and successive computation of /'”) from (3.14). We can, however, 
give a formal treatment, from which the time-dependence of the complete 
¢ becomes clear. 

The first equation ( (3.13a) for r=1 is formally solved by developing 


ty 


dM (t), as well as AM (t=0)=4(0), in eigenfunctions of Q: 


(3.19) b (e;t) 2 > > puter (t) ion (c) 


and 


~~ ~~ 


(3.20) $@(c;t= 0) = ¢4(0) = ¥ Yo y4 (c) with po” = (4(0), 7). 
o A 


g4) is the initial value of 9 !4: 


(3.21) G0 (4 (¢ = 0) = GM, 
Substituting (3.19) into (3.13a), we have for each mA: 
(3.22) (24+) GVA(t) = 0, 

with the initial condition (3.21); the solution reads: 
(3.23) PVlOA(t) = God eo 


and the full oa (3.19) is correspondingly given by 


(3.24) = ae ge wy [wal ( => 2 yu )aol ( e-ot, 


wo 


in which the W-notation is used in order to group the terms of o according 
to their exp(— wt)-factors. This simply means that each m/-mode of the 
initial 4(0) decays with 1 [wm as its relaxation time. 

Next we at determine ¢@). Substitute 7@ = fd}, as given . by (3.24), 
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back into the definition of #®@ (3.14, r=2); we find: 


(3.25) 


Fo 1 J (fF? fi?) — > G2) [a] p—ot. 
o>0 


in which the sums of any two w occurring in (3.24), which are of course 
all +0, have been arranged in order of increasing size and called ©. 
Now write out eq. (3.13b) for r=2 

(3.26) (4 4 Q) go 2) p G2) [1 e- ot _ (0). 


o 


Develop both ¢” and G®! in y-series: 


(3.27) G2)() (¢) —= mY >) g (2) (ww) yor (c) 
wo A 

and 

(3.28) db? (e; t) => d¢ pp'2) (mal / )y (c c): 


initially (ef. 3.15): 
(3.29) p41 (t= 0)=0. 


Substitute (3:27, 28) into (3.26); for each mA, 


(3.30) (d1 A. w) pi) taal (¢ )- g@itwodl get = ge gaiceh y= == {). 

@+o 
Here we have supposed that the eigenvalue w is found again among 
the sums-of-eigenvalues ® occurring in the development (3.25) of F®), 
Of course, if this is not true, we simply omit g®!°°!, The solution of 
(3.30) is easily seen to have the form: 
(3.31) pi?) wad (¢ ) = (2) foal fp’ 2)loAlt)e oti Sy (2) weal ot 


— 


Oro 


In order to determine the coefficients y, we substitute this back into 
(3.30) and put coefficients of each exp (— ft) separately =0: 


9) OF "(2 A 2 A . 
(3.32a) y'| ) [w. 1_ gi ) (eneal = () ; 
(3.32b) (co — &)p® loa __g(2) (oh) _ ), 
while the initial condition (3.29) serves to fix the one remaining coefficient: 


(3.32c) (2) Coal 4. > 2) worl — 0), 
Ot 

or, using (3.32b), 

(3.32d) prion gi?) low} 


= o—o 
oro 
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Putting (3.32a, b, d) in (3.31) and going back to (3.28), we find: 


\ be Nees {== 2, aH y, cov (c) (9 2) [ood] fo —at 1 SS g{?) Lowa) a 
(3 33) ae o—w 
|] — > > iA) (¢) PS tis g ) [awa] 1 Flee ( \. 


\ wo 4 
\ 


in which we have introduced: 


\ (®4o) Blve! (t) — Flo) (t) et a : 
(3.34) d o—w 
Eo! (t) = lim Hl! (t) = te, 


As for the time-dependent factors of the various modes, t exp (—«) 
starts from 0 at ¢=0 at a rate 1, to reach its maximum value I/we at 
t=1/w; from that time, it decreases rapidly. H(t) behaves very much 
like texp (—of): it also starts from 0 at a rate 1; it stays in between 
texp (—ot) and texp (—@f), reaching a maximum in between those of 
texp (—ot) and texp(—ot) at a time between 1/m and 1/6. This is 
how all modes behave. 

Finally, we rewrite (3.33) in a notation arranging the terms of p® 
according to their t-dependent factors exp (—t) and ¢t exp (—@t): 


(3.35) P(6 51) (PON) Beene chi) en 
Having dealt with $@ in such detail, we can be short in discussing 
higher terms. Next comes 


\ F@)(c; t) = (1/f O)LT(FOA,@) +I(fOf My} = 
(3.36) / = > (Gel 4 GB lal) eat 
as follows on substituting /® and f®, as given by (3: 24, 35), into 
(3.14, r=3); this time @ represents triple sums of eigenvalues m. The 
solution $®) of eq. (3.13b) for r=3 will read: 


(3.37) b(c; t)= Y (Yc) 4 YO") (0) (ep 4 Yo") (01 (¢) £2) e7 ot, 
Continuing in this way, one finds that the general F) (6) consist of 
terms: exp (— Gt) times a polynomial in ¢ of maximal degree r-- 2 (r—1). 
Of course, whether these degrees are actually reached, depends on the 
recurrence of the eigenvalues » among the @, which are r-fold sums of 
w values. 

Thus, ¢ goes to zero for f -> co, because all terms in the developments 
of J, é2) etc. contain decreasing powers exp (—t), where the @ are 
those w which occur in the y-development (3.20) of the initial (0), and 
sums @ of these w. The smallest of all possible @ is m1, the first positive 
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eigenvalue of 2, and its inverse is the longest of the decay times of the 
terms in d, which is the characteristic relaxation time 7 for the kinetic 
stage: 

(3.38) T= 1/oK, 


Of course, the effective relaxation time is somewhat smaller if one or 
more of the lowest @ are a from (3.20). A discussion of the time 
scale t will be given in ch. 

Finally, we should say sinatties on the conv ergence of the series 
(3.8). The general character of the terms dO, éd®,... is such that 
bV(e; t) initially equals (0) and then decays, for all ¢, with a relaxation 
time  T; p® and higher terms start at 0 for all ¢; in a time of order r, 
they build up some maximal value (as functions of f, for all ¢) and then 
they decay exponentially like é) does. Now in order to compare the 
terms /,/@,/@,..., which are functions of ¢ and t, we consider in 
particular a time t ~ 1 (because, if f > rt, all terms following / vanish 
anyway) and a velocity ¢ ~ ¢o (because, in calculating physical properties 
like entropy production and heat or momentum transport, one or two 
factors ¢ are put in a scalar product, i.e., integrated after multiplication 
by exp {—m/(¢—c¢o)?/2k7}). First, we substitute these values in the series 
expansion of $) (3.24); supposing this does contain the first positive 
eigenvalue m,, and neglecting higher modes m2,... on account of the 
accompanying factors exp (— wot), ..., we get: 

(3.39a) b(eo; T) & (Lfe)pleiy'4) (eg) & /e)b(eo; t=0) )=@. 
Likewise, it appears that 
(3.39b) b)(€o; tT) ef 


etc.: l/e times the initial value of the deviation 4 at ¢=c¢o serves as a 
parameter setting the pace for the convergence of the ¢ series. 


4. Velocity distribution in the kinetic stage: non-uniform case 

After the special (ChE) solution of the r-dependent case (ch. 2) and 
the general solution of the spatially uniform case (ch. 3), we shall now 
determine the general non-uniform velocity-space distribution f(¢; r, f) 
from the BE, which reads (2.1): 


(4.1) de f(€s r,t) =dsf +I (ff). 


Many steps in this solution find their parallels in ch. 2, in ch. 3, or in both. 

In the ChE method, the solution / had the special characteristic of 
being a functional of the py; this, together with (4.1), proved sufficient 
to determine /. Without such restriction, we are free to impose an initial 
condition: 


(4.2) f(¢; r; t=0)=f(0). 
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Let the solution f(¢; r,t) have y“(r,t) for its first 5 moments. On the 
base of these y, one can determine a complete ChE distribution 7°, from 
which / differs by an amount 


(4.3) f-fc=fle; r, t)=f Od, 


/® denoting the MB distribution based, like f, on yp (in fact, this is 
the first term in the series which constitutes /C®), We do not require the 
difference # to be small. At t=0, it has a given value 


(4.4) j(e; r; t=0) =f —fen(0) —f ©(0)4(0). 

Both f and /“ have y"™! for their first 5 moments, so that those of } vanish : 
(4.5) wr, t) = fde fy = (d. wie) = (): 

in particular, 

(4.6) pr; =0) =(4(0), x) =0. 


It is this ¢ (or f) for which we are going to draw up an equation. Because 
of the definition (4.3), this imphes that we should know how f and fC 
change in time. f obeys to the BE (4.1); 0; f" does not follow from a BE, 
because it is not a real distribution but an implicite function of r and t 
through the y, so that we should start from the HDE for y (see 2.10). 
The RHS of these equations contain the 5 moments y™”! themselves, as 
well as the higher moments P and q, the latter ones linearly. Now we 
first determine 0; y computing these moments y, P,q on the base of 


{Ch only: this means inserting the exact y but ChE values PC, q°h, 


ie) ~ 3 
Then we only have to add terms with P and q arising from the correction 
Ld a . =< Dos . 
/, because f does not change the w values (see 4.5). In all, this gives 


~ 


(4.7) ay yl — fde 7D. (fch4 “f= ,Ch pil + wp, 


in which the second equality, holding term by term, serves to define the 
‘partial’ t-derivatives 2,°2 and 0, of the y”. The change in f°, a 
functional of the y, can likewise be divided into two parts: 


6fch : ; re 
(4.8) ny fen = Soo ny yi = 9,Ch fen L Dy i Ch, 
2 


if the obvious definitions hold: 


ou \ . 
(4.9) ¢°8 w(cly'(r, t)) = 2 owl) oe, 


(4.9b) de w(e|pl"(r, £)) = > oa ey. 


Eq. (4.8) gives a comparison between the actual time derivative of f° 
and the one we would find if the distribution were really {© and no 


L32 
more. The latter follows from the BE as it reads in the ChE case: 
(4.10) DeCh fCh= do fCh+J(f Caf, C2), 


whereas, for the full distribution /, we have the BE (4.1) which reads, 
introducing f by means of (4.3): 


(4.11) (f+ fy =I (fM+fy+I(fr+f, e+ fr). 


Combination of eqs. (4: 8, 10, 11) gives the equation for d, / and, conse- 


~ 


quently, that for 0; ¢: 


~ 


(4.12) — uf=fO dip + def © + db dif © =dsf +I (fA) +I (fA) +I (ff). 


The next two steps are, of course, series development of the various 
quantities and insertion of powers of 6=1, in order to get a hierarchy 
of equations to be solved in succession. This time, we might as well 
combine these two steps and write immediately: 


(4.13a) foha(1/6) f+f+...=(1/6) f (14+ 046 +...); 
(4.13b) an ee f (06M + A262) + ...), 
with the normalization: 

(4.14a) pA m= fde fry — yl5,9; 

(4.14b) Pa = fde fryr—Q (ry > 1). 


Higher moments of /, like P and q, also consist of contributions from the 
terms of (4.138a, b), and so do the RHS of the HDE: 


(4.130) dg hyll = fede yl (Ff 4 f+ ...) = (94 + O94 +...) 
(4.13d) De pl = fde AQ (fO 4 f24 y= (O34) + 6342+... yh, 
as well as the “partial” ¢-derivatives of any functional wu of yp: 
(4.13e) O42 w= (de + O9¢M +... ue; 
(4.13f) of u= ne + 62342) + sas JU, 


if, like in (4.9), we define: 


(4. 15a) de) u(e|p) = > a dy"); 
(4.15b) 4) w (¢|w) = SP ea dep, 


As regards d4 f = d4( (fd), because /© is t- -dependent (unlike the uniform 
/), we must assign the factor 0" either to d;/ +) or to d, é”. The latter 
choice seems preferable: 


(4.13g) i b=d1 69 4 O04 G2 +... 
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Substitute (4.13a, b, e, f, g) into (4.12) and put coefficients of each 


power 6” separately =0; for the J terms, we use an obvious shorthand 
notation : 


\ [1/0] fr dM =J( / (0)7,@) 4 J (Ff, ) = J(O1; 10) 
[0°] WFO +f OIG (2) +40), (0) =>. f@ 4+J(02:11;11; 20; TM) 
de DFM +9,@f O + cae +409 “ 0) + B20) O) 4 

1 


P( 


+$03,MF =a. F@ +5(03; 12; 21; 12; 91; 30; 13; 31), ete. 


There is a complete analogy to the comparable equations in the uniform 
case (3.12) in that operations on the unknown 4 include the linearized 
collision operator in ¢ space and differentiation 2;. This leads to an 
abbreviated notation identical with (3.13a, b): 


~ 


(4.17a) 2 6M = — 240; 


~ 


(4.17b) 4 OM = FO ~ 04 


in which, this time, both #™, which has the new meaning: 


aS = fio) fa 
FO) =(1/f){9.f ¢-) — 0-1-2) f (2) — (2), 7—-1-e) (0) — 
( if yi sl] t t 
(4.18) a ae or se 
Se Ore eh DIG hi OO = 0s 0.70) 
o=1 1 


oe 
and 4) depend parametrically on r. The initial condition (4.4) can again 
be imposed on ¢ alone: 


(4.19) AM(c; r; t=0)=4(0)dr0. 
Normalization of 4 must take place according to (4.14b): 
(4.20) (A, yi) =0. 


If we now develop all functions occurring in the general eq. (4.17b) in 
eigenfunctions y!] of 2, d% é” will not contain collisional invariants 
741 because, by (4.20), ), 6 does not at any time, and 246 is of course 


ae to 7%], so that again 
(4.21) (FO), (041) = 0 


is a condition of solubility. As in the ChE case, (4.21) for a given r is 
proved by the aid of (4.20) for 1, 2,...,7—1, as follows: 


Cited ie 
(FO, '04]) = fde 104] Of) = > fde 04] o47-1-2) f © fod 
o=0 
—1 75, pw) 
FE Ss fde 10) fh (2) 9,7 -1-8) f (0) == ys fde plore) 9 g(r —1-0) ¢ (0) 4 


(4.22) = 


— 


os 2! [de 71. J(o, r—@; 7-0; & —0)= 
1) + (IL) + (11) +(1V) +(V), 


my iss} 
a 
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in which we work out the single terms as follows: 


(I) =d:0-Dyl4) by (4.13d); 


r=2 


Ce cee Bye ty a fide lon by (4.15b) 
= 5 Dap (see 2: a 53); 
493) | HD=- by Ze spi Jade Goin © by (4.13¢) 
=0 by (4.14b); 
(IV) =- 55 ay -ehylt) og Jed ier zivrif © by (4.13d) 


=(0, again by (4.14b); 
(V) =0, by (2: 9, 51). 


So finally (1)+(11)=0, which proves (4.21). 

We conclude that the equations (4.17), with the normalization (4.20), 
the solubility conditions (4.21) and the initial condition (4.19) — normalized 
by (4.6) —have exactly the same character as the corresponding expressions 
for the r-independent case (3: 13, 17, 16, 15, 18), the only differences being 
(a) that all functions depend parametrically on r, and (b) that F© in the 
present case (4.18) is still more cumbersome to compute than the spatially 
uniform F”) defined in (3.14). In fact, F as determined from f, f®, 

.,/¢-) now contains terms linear in the /@ as well as quadratic, so 
that already the expansion of F@) in decreasing e-powers (3.25): 

FO) = ¥ GO) eo 
will contain -values which are simple eigenvalues @, as well as binary 
sums of these m. However, the occurrence of simple w values among the 
sums @ has already been taken into account in the treatment of ch. 38, 


so that the results obtained there fully apply to the present, r-dependent 
case. 


The terms of the deviation of / from a normal distribution of ChE type: 
j/=f—f%, thus decay according to equations like: 


(3.24) for ¢@; 
(3: 28, 33, 35) for J@); 
(3.37) for 6), etc. 
Like in the spatially uniform case (cf. eq. 3.39a), the first-order deviation 


in r at t=1; c=cp is of the order-of-magnitude of the initial (t=0) 
deviation in r at c=oq: 


(4.24) pO(r; €o; t) (Le) hr; e; t=0)=9(r), 


This time, analogous equations (like 3.39b) for ¢@ ete. only hold if the 
terms of each “‘partial BE” (4.16) are of the same order-of-magnitude, 
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so that the ‘‘@-procedure” (eqs. 4.13) is not merely formal. Only then 
the local 9(r)-value serves again as a convergence parameter. This require- 
ment is met if not only the initial 4(r) is comparable to the local Chapman 
é(r), but also similar relations hold between successive gradients (d/dr) 
of ¢ and ¢. This of course limits the range of initial values ¢(r; t=0), for 
which the given procedure results in a quickly converging expansion of 
the distribution function at later times. Physically, this means that, even 
if the initial ¢ is <1 everywhere, a smooth approach to the ChE 
distribution for all r is only guaranteed if restrictions are applied to the 
gradients of ¢ (t=0), whereas, in other cases, local disturbances may 
travel through the gas before they are eventually damped out. Such 
processes, though formally covered by our treatment, will be treated 
with better hope of quick convergence by choosing a different “0-grouping”’ 
method. 

As a conclusion to this treatment of the general non-uniform kinetic 
process, we may say that this theory provides an argument in favour 
of statements like that of CHapMAN and CowLIna ([2], p. 117), who said 
that the velocity distribution in a gas, starting at an arbitrary initial 
value, moves towards a normal value; and that, on physical grounds, 
this normal distribution has to be the one of ChE, which is fully determined 
once one knows the five moments y!’] as functions of r. 


(To be continued) 
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5. Entropy production 


In order to discuss the entropy density and its development in time, 
we start from the following definitions: 

a) u-function: as usual in statistical mechanics, we instead of ¢ take 
the momentum 


(5.1) p=me 


as a variable. In a space of momentum p and position r, we consider a 
dimensionless distribution f, defined by 


(5.2) F (pi r,t) dp ae h3f(c;r,t)de, or f(p) = (h/m)3f(c) 


(h=Planck’s constant); in terms of this f, we define the #-function, 
following PRIGOGINE [8], by 


(5.3) u(r)= —h-3 f[dpf(1—Inf), 
or, back in velocity space: 
(5.4) a(r)= —Jf de f{1+ 3 In (m/h)—In Ff}. 


The integral f{ dru is dimensionless, so 7 can be looked upon as the local 
density of a dimensionless probability function. 
b) S-density: the entropy per unit volume is defined as 


(5.5) s(r)= —ku(r)=k J de f{1+3 In (m/h) —In f}. 


This definition differs from the usual one (with —/ In f in the integrand) 
by the addition of 1+3 In (m/h) to —In f. This choice has the advantage 
of leading back, in the case of equilibrium, to the Sackur—Tetrode formula 
for the entropy content of an ideal gas in a volume V, as derived by 
statistical mechanics: 


(5.6) S=s0V=knV In ( 1/n) amy” e'n| : 


h2 
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We now take a non-equilibrium distribution 
(5.7) fe; r, )=f© (c) (1+¢4(¢; r, t)); 


for the time being, ¢ will denote an arbitrary deviation (of ChE type or 
not) from a local MB distribution. With 


(5.8) In f© =In {n(m/2akT)"2} —C*2 = y—C*2 
and the expansion 
(5.9) In (1+¢)=¢—4¢?+ -..., 


we get the following expression for s(r), in which the integrand has been 
arranged according to powers of ¢: 


(s=k fde f[{1+3 In (m/h) —y +O0*2}4 


D1 
(5.10) ( + 4{3 In (m/h) —y +C*?}—3¢?+...]. 


Now because f and f® must have the same first 5 moments yl], ¢ is 
normalized by 


(5.11) [de fOd=f de fdC*2=0, 


so that the contribution linear in ¢ is zero: 


(s=k f de f[{1+3 In (m/h) —y+C**} -3P'+...] = 


(5.12) ) == (0) 1 52) 4 0)4 


Here, s® represents s as found by using only the MB term f of f; s® 
forms a first, essentially non-positive correction. It is easily seen that the 
insertion of a development for ¢, like 


(5.13) f=fO(14+46M+44@4...) 


of terms decreasing in order of magnitude leads to the analogous result: 


(5.14) s=s Lk fde fg +higher terms, 


the contributions of 42, 4° etc. being confined to these higher terms. 
The time derivative of s reads: 


(5.15) d48=k f de d% f{3 In (m/h) —In fF}. 
By inserting the BE: 
0 0 
(516) m%f=r%mftrof af = —< Sh —F. Shs tof = JIA) 


we divide d;s into a streaming part: 


(5.17) Iss=k f de wf{3 In (m/h) —In fF}, 
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which can have either sign, and a collision part: 
(5.18) Ios =k f de J(ffi){3 In (m/h) —In Ff}; 
the constant 3 In (m/h) can be omitted on account of (2.9): 
(5.19) I W8= —k f de I (ffi) Inf. 


In the usual proof of the #-theorem (see, ¢.g. [2], p. 69), it is shown that 
this integral is essentially <0, so that 


(5.20) do8 > 0. 


We shall now write d,s, as given by (5.20), in terms of ¢, like we did 
to s itself in (5.10). Again by (2.9), we diminish In f by In /, which is a 
collisional invariant (ef. 5.8); J(f/;) is of course zero; develop In (1+ ¢) 
in a power series and introduce 2 (ef. 2.32): 


(5.21) ds=k f de{fO2d—JS(f Md, f:1di)}(6—4$4? + $6? — +..-). 
For small ¢, 


(5.22) 


(bracket notation introduced in (2.37)), which is >0, just like the exact 
d,s. In these formulae, ¢ always represents the relative difference between 
the actual, arbitrarily chosen / and the corresponding local MB distribution 
(see 5.7). This quantity we called ¢ in the ChE case (ch. 2); é in the 
general uniform case (ch. 3) and $+¢ in the general non-uniform case 
(ch. 4). For each of these cases, numbered III, ILA, Il according to the 
classification of ch. 1, we shall now derive approximative expressions 
showing how changes of s by streaming and by collisions (5: 17, 18) 
influence the terms of the series (5.12). To this end, we write down each 
time: 

a) the full time derivatives of the first 2 terms of (5.12), s© and s@); 

b) the s- and c-parts of the t-derivative of the full s. 
Of course, 


(5.23) de(8 + 8@ 4 s@) +...) =dIs8+ dIo8 


and we want to compare the separate terms on both sides. For shortness, 
we will use the abbreviation 


(5.24) 3 In (m/h) —In f© = K(C). 


K is a dimensionless function of the length C of the vector of relative 
velocity C, of order-of-magnitude 1, reaching a maximum of about 10 


for C=0 in a typical case (He at NTP); K is a linear combination of 
collisional invariants 710), 
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As to the first remark in column. 5, 
as Of (9) AG 0) — ofl P 4 
(5.26) yl Mee 2 5ym dep]; 9, FO = > 5 yi d2Oypl4), 


and all 5 d//dpl4l are some collisional invariant 7! times /. All 
approximations in column 4 are properties of the ChE correction ¢, valid 
in particular if the terms of 

(5.27) fon=fO(1+4¢)=f/%1+46%+¢42 +...) 

converge rapidly. If this is actually the case, we may also suppose ¢ 
to approximately satisfy the equation for 4%, (2.33a): 

(5.28) 9, fF) zw Isf O —f OO. 

Now it is evident from column 3 that the approximations given in (5.25) 
for the terms of (5.23) balance. In fact, we conclude: 


(5.29a) d¢ 8 w& Iss+ C8; 
(5.29b) d¢ S@) w 488)... & 0, 
with, in particular, 

(5.29c) Ics ~ k[d, d]: 


at all times, s exceeds s®@ by two orders-of-magnitude; changes in s 
both by streaming and by collision are mainly used to adapt s to the 
changes in /® during the HD development. 

As we already stated, the sign of dss depends on the case considered. 
In a typical transport experiment, during which a gradient is maintained, 
it will even be negative, because streaming has to make up for the 
constant increase of s by collisions (5.29¢c). As an example, in the case 
of heat transport q down a constant temperature gradient 7’, the values 
of column 3 read: 


(5.30) 48 =d48@) = 0; ss= —T 19g -VI'= —A(V Nn T)?; dos = +A(V NM TP 


(A= coefficient of heat conductivity). 

For a detailed discussion of the entropy production in the ChE case, 
and of the restrictions on the validity of the approximative result (5.29a), 
see PRIGOGINE [8]. 


II A. General uniform distribution: f (14-4) 


(1) (2) (3) (4) (5) 
term explicit expression ee i giro ain remarl 
wO=kf dere f-K = 0 = 4 /=0 
d¢8) = —kfde($r./ -$2+/OSard) =—kJde/Mbred . df =0 
Iss =k fdedsf {K—In(1+¢)} = 0 asi te he0 


~~ 


des =k fdeJ(ffr) {K —In (14+4)} wk fde/-0d aes a ae 
- —s) 1 
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The approximations in column 4 rest on ¢ <1 only. Again, if d is expanded 
in a convergent series: 


(5.32) o=b%4+46@+..., 
d will approximately satisfy the equation for its first term ¢@ (3.13a): 
(5.33) 4d ~ —Qd, 


so that the balance (5.23) now reads, in detail: 
(5.34) Ot g (0) — OF OSS OF o¢ 82) x dCs, 
with 

dcs  k[d, 4]. 
This time s does not change by streaming. Collisions, however, cause 
changes in s: not in s which is tied to the constant /, but in s@) which, 
containing ¢@), vanishes as ¢ dies out by relaxation. (Of course, the same 
holds for s®) ete., which contain higher powers of ¢). s®@ dies out twice 
as fast as ¢: if the longest-living mode of ¢ is (cf. 3.24) 


(5.35a) bY) ~ WAll(e)e-t, 
then 

(5.35b) s2) ~ g@)[l(e)e-2t, 

in which 

(5.35c) o2)(4— —1k( YO), Yale), 


(~ means: asymptotically equal to). 

The formulae of this section may serve to illustrate the usual proofs 
of the #-theorem (like in [2], p. 69), by giving the behaviour of ¢, 7 and s 
in detail, for the case of small initial ¢. 


II. General non-uniform distribution: f(1+¢+¢) 
(5.36) (see following page) 


This time, the approximations in column 3 satisfy the balance (5.23) 
if both ¢ and d approximately fit the equations for their first terms 
(5: 28, 33), which is the case if the series developments (5: 27, 32) for fC 
and f both converge reasonably .In detail, we thus find: 


d¢ 5) 
O¢ g (2) 


~ 


Is s+kl[o, o+¢] 
kd, b 5 $| 


@ 


(5.37) ee Asn F/O, 
with (d) ws kidt+¢,o+4l. 


Conclusions, in this approximation, are: s changes in t owing to both 
streaming and collisions; s®) is only influenced by collisions (changes in 
s®) and higher terms, although of course present, are not noticeable on 
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the present scale). When d has eventually died out, as described in ch. 4 
the ordinary ChE balance (5.29) remains. 

Comparing the three cases: IIT, 6; ITA, 6; H, 6+ ¢, we conclude: 

1) In the cases of s® changing in time (III, II), ¢s contains a 
“streaming” contribution dss which can have either sign, 

2) The collisional contribution des, which is essentially positive, finds 
expression in: 

a) an increase k[¢?, ¢]>0 in s® if 640; 

b) an increase k[d, o]>0 in s® if +0, to which should be added 

c) equal changes k[¢, 6]=k[¢, 6] (ef. 2.36b) in both s and s®), if 
d+0 as well as $40. 

We shall next look at the irreversible entropy production dss from the 
viewpoint of the thermodynamics of irreversible processes. 


) 


ITI. ChE distribution f(1+¢): we introduce transport ‘“‘forces’ X 
and corresponding “currents” J as follows: for the viscosity, 
ie) 


(5.38) Xy= —eo and i= (tensor of HD pressure) 


(the symbol — serves to symmetrize a tensor: To h(T g's.) and 
for the heat conductivity, 


(5.39) X,=—VInT and J,=q (vector of heat transport). 


These forces occur in the LHS of the integral equation (2.33a) for 4: 


(5.40) (§—0*2)C.7 In T—2€*C*: Pep= 24 


and, consequently, in the solution 
oO 


(5.41) dU =A-V InT+B: eo; 


equalizing coefficients of each separate gradient, we have: 
(5.42) (§—O0*2)C=QA; —2C*C* =2B. 


The LHS functions in these integral equations happen to return in the 
definitions of the currents J; therefore 


(5.43a) 1 =P= kT, —2C*C*) = — kT, Bi; 
(5.43b) h=q= —kT(4, (§-072)C)= — kT Id, Al. 


Multiply each equation by 1/7’ times the corresponding X and add the 
results, supposing the ¢-equation (5.41) to approximately hold for ¢: 


(5.44) TUS: X,+Jy-%) © kd, $1 © 208. 


Here the first ~ sign comes in because of the approximation $ ~ pf) ; 
the second one is introduced by the linearization of the collision operator, 
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which presupposes ¢ < 1. These neglections are both allowed if the series 
(5.27) is duly convergent. 

The irreversible s-production has thus been written as a sum of products 
of streams and corresponding forces. The balance (5.29) can accordingly 
also be formulated as: 


\ 048) we Iss ei fiee d+ IX 


5,45 
ee) (d8@ w... w 0. 


Using the phenomenological relations, which connect the currents J to 
the forces by means of the transport coefficients (~ of viscosity, A of 
heat conductivity) : 

( d= 2X, = Ly Xp; 
~ > t i ih pu 7] 
(5.46) ; 
(J,=ATX,= LyX; 
one can also write dJcs in terms of the phenomenological coefficients L: 


(5.47) do8 wT (LyX? + LyX"). 


II. General r-dependent distribution f(1+¢4 6): As in the ChE case 
(III), we consider the sum (5.44) 


Ce LP XK, +Ja: X,;). 


in which the transport currents now read: 


eas (J, =P= —kP(6+¢, —2C*C*) = — [6 +4, _ 
a (S=q= —kT($ +4, (§—C*)C) = —kT[p +4, A 
which leads to 

(5.49) TASIX ~ kid+¢, 6] 


instead of (5.44). Comparing this result with the balance (5.37), we see 
that, like in the ChE case (5.45), 7-1! ¥ JX renders the contribution of 
collisions to the change in 8: 


(5.50) 0,8) ~~ ss -4 J iis! > JX. 
This latter contribution can now have either sign; together with the 
collisional change in s@), of course, it forms des which is always positive. 


Finally, we point out the following difference between the ChE (4) 
case and the general (¢+¢4) one: in the ¢ case (5: 44, 47), 


(5.51) cs=T-1 ¥ JX=T- ¥ LXX; 

in the ¢+¢ case, however, these three quantities all differ: (5: 37d, 49): 
(5.524) dcs=k[b+4, +4] 

(5.52b) TS IX=k[$+4, 4] 


(5.52c) T } LXX=k{4, 4], 
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(ec) holding because the relation between the ChE-¢ and the gradients 
is the one (5.41) valid in the case of 4 alone. Looking at eqs. (5.52a, b, c), 
one is tempted to suggest that des could be identified with 7 eae, 
(L being the inverse of the LZ matrix), but this is not correct. 


6. Transport currents; quasi-relaxation 

We shall now describe the behaviour of the transport currents JP 
and J,=q during the kinetic stage, like we did for the entropy density 
in the preceding chapter. The treatment will correspond to section II of 
ch. 5, for we shall immediately consider the general non-uniform distribu- 
tion ((14+¢+¢4), of which the ¢- and ¢-distributions are, of course, 
special cases. 

We introduce functions G and H as follows: 


(6.1) SIC Ch=G (0 )G=c, 
(6.2) H, = —kTG, = —mCC; H,= —kTG,= (4m02 — 3k TIC. 


Using these H, we can write the transport currents J as follows (cf. 5.48): 


\ P=J,,=f de fH, =(H,, 6+4) 


(6.3) ( q=J,=J de fH,=(H;, +9), 


in which we have omitted terms 


(6.4) (H,.1)=0; (Hy, 1)=0 


which vanish because the MB term f/) does not contribute to P and q 
(cf. 2.24). 

Apparently, both currents J consist of the ChE contribution (//, ¢) + 
a deviation (H, d) which, like ¢ itself, dies out during the kinetic stage. 
A formal expression for the t-dependence of J will be given later 
(see 6.14a). 

Next, we make a division of the derivative 0,J, in the way we did for 
Ds ia (62-17, 18); 


(6.5) Del =04 if de ji =f de (dsf =i dcf \H +f de i 0-H =ded + ded + ond 5 


here dsJ derives from changes in f by streaming; dcJ from changes in 
f by collisions; 047 comes in because H,, and H, are functions of ¢ depending 
parametrically on the moments y (from (6.2) it is seen that H, contains 
co, H, contains both ¢9 and 7’); hence: 

0 
(6.6) dpi (es yl) — > att dpi! 


We now make approximations for the terms of (6.5), just like we did 
for d¢s in (5.36): 
10 Series B 


146 


(1) (2) (3) (4) 
term ee approximate value approximations used 
(0) 
ud =fded (fH) ~ J de fOH{ud b+d K1;fObD Ku fs de ie re 940) 
wd =f[dersf-H ~ffdersf/- HA Is f Ys fO 
deT =IdeI(fp)H ~ fde {—/Q(b+9)} H Tih) ~ —fO2(b+¢) 
iJ=fdefurH wfdefd,0 A 4+¢<1;4H ~ >, H. 


Here 9, applied to a function u(¢; y) has the meaning given in (2.28), 
where d;) applied to a functional u(¢|y) was defined; in this simple case, 


> ou . 
(6.8) D4) u(c; y) — z Syl D4) yl), 
A 


If 6+¢<1 as supposed, the differences yy—dyp and dxu—1w are 
relatively small for any wv. This and the other approximations of column 4 
lead to the result for d;/ given in column 3 of (6.7) in the following way: 


~ 


: (dF =f de[H{(1 + $+) uf O+f Up+fd 62 bias 
cata ww f def{d,(f OH) + (OH dd} d(H, 1) +(H, did), 


of which the first term vanishes by (6.4). 
Resuming (6.7) and combining the ds and dq contributions, which are 
both, in a sense, streaming effects, we find: 


(9s +dx)J a A(dsf © —¢f ©) =(H, FO) 
(6.10) dct «=o pw —(H, Q[6+4]) 
J «=e (A, a 


Here F® is the LHS (2.34a) of the integral equation (2.33a) for 4; 
like the d,s-balance (5.36), the present one is in order on condition that 
¢ approximately satisfies the ¢“)-equation and ¢ the one for 6 (eqs. 
5: 28, 33), ie. on condition that the series developments of 4 and ¢ are 
both reasonably convergent. This leads to the following survey, listing 
the terms of (6.10) for each of the 3 distributions treated in chs. 2, 3, 4: 


—(H,¢+¢]| —(H,4] 
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We draw the following conclusions with regard to the development of 
J =(H,4+4) in time (stages indicated by roman numerals as in chs. 1 
and 5): 

II A. Kinetic stage in a spatially uniform gas, ¢ only: J is not changed 
by streaming, but collisions cause a change in J to the amount 


(6.12a) de (H, $) ~ —[H, 4] . 


This equation indicates that J dies out due to a relaxation process, to 
which we shall return presently. 

II. Kinetic process in a non-uniform gas, ¢+4¢: streaming changes 
tend to build up J during the kinetic stage at a constant rate depending 
on the existing y gradients: 


(6.1 1a) QQst+dn) J w~ (H, F®) w [H, 6] 


but collisions try to annihilate J at a rate increasing together with the 
value of J itself: 


(6.12b) 


Together, 
(6.11b) ud ~ —[H, d]. 


¢ does not change but ¢ dies out rapidly. 
III. As the kinetic stage comes to an end and the next (HD) stage 
sets in, 


(6.12c) ded —[H, P| » —(dst+du)J: 


streaming and collision contributions cancel. As f reaches the ChE value, 
the transport J reaches the HD value (/, 4). 

In each of these 3 cases, the collisional derivatives of J are related to J 
itself by equations of the same type (6.12a, b,c). Only in case II A 
collisions are not counteracted by streaming; in case II there is such a 
counteraction, in case III there is even a balance. We would like therefore 
to indicate the phenomenon of (6.12) as quasi-relaxation. 

The implications of eq. (6.12b) will become clearer if, using the 
approximations dx?) and ¢~ 4 also in J itself, we develop 
6,6 and H in Q-eigenfunctions y. As for H, 


(6.13) H(c)= Pp: Allyl with hol —(H, zt), 


in which, as usual, w>0, because all H are orthogonal to the collision- 
invariant functions 7!) of (2.44). The developments of 6M and od” are 


given in (2.47) and (3.24): 
(= ()[oa]yloA]. f(t) — plorle—otyloa] ; 
d >> (1/a) f get; D2 e yt 
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they contain the 7-series coefficients of F“) and of d(t=0). Substitution 
of these series into (6.3) and (6.12b) gives: 


(6. 14a) i ~ Sy, { ( | /@) f G)[4] } qiorle—ot} | = bey J () ’ 


wo as os 


(6.14b) ded & eae (f Dla] 4 Glolme—ot) Alea] == ape, dod lo); 


wo a bd 


here, J and dcJ have been written as sums of ,,modes’’, of which the one 
numbered wA decays according to the relaxation equation 


(6.15) (doJ lo = — wJ lI, 


Again, 1/@ assumes the role of a relaxation time for the mA mode. In 
the next chapter we will give one-term approximations to the develop- 
ments (6.13) of both H,, and H;, which of course greatly simplify matters, 
leaving us with two quasi-relaxation times 1/@ only, one for P and one 
for q (t, and 1,). 

Mention of something like this quasi-relaxation was first made by 
KouLER [6], who dealt with the change-by-collisions of the 3 quantities 


(6.16) 6,,=(m/k) f de f (C2 —4C?) (x=2, y, 2) 


which he characterized.as 7',—7’, the excess of the “‘temperature”’ 7’, of 
the «-component of heat motion over 7=} 57, during the HD stage. 
Clearly : 


(6.17) KT —T,)=—P 


Xa? 


the diagonal element of the HD pressure; therefore Kohler’s result is 
contained in the present treatment. From his formula, reading in our 
notation: 


(6.18) dc0, = —0,/t 


x P 
Kohler concludes: 

a) Tt, measures the speed of mutual exchange of energy between the 
degrees of freedom a; 

b) the viscosity m arises as a consequence of the finite time of 
establishment of the ‘‘one-direction-temperatures” ,. 

Instead, we would prefer the following statements: 


a) The quasi-relaxation phenomenon involving the second-degree 


moment P has to do with the transport of a first-degree function of C, 
i.e. with the transport of momentum. 7, measures the speed of the exchange 
of momentum in one direction « bekinnen colliding molecules, like t,, 
measures the transport of kinetic energy; 


b) the viscosity j results from the equilibrium between increase 


(by streaming) and breakdown (by collisions) of the pressure tensor P. 


Incidentally, Kohler complicates his argument unnecessarily by the 
use of so-called main axes (a system in which the eo tensor is diagonal). 


(To be continued) 


PELYSUGS 


THEORY OF RELAXATION PHENOMENA IN A 
MONATOMIC GAS. IV 


BY 


M. J. OFFERHAUS 


(Communicated by Prof. J. p—E Borr at the meeting of October 29, 1960) 


7. Collision operator; time scale of the kinetic stage 

Numerical description of the relaxation processes discussed in the 
preceding chapters requires solving the eigenvalue problem of the Q 
operator (2.39). The results will largely depend on the nature of the 
(spherically symmetric) intermolecular potential field, V(r), and, for 
general V, no complete solutions are known. There is one property of 
this general solution, however, which follows from the invariance of Q 
with respect to rotations in ¢ space about the point ¢o. It is easier to 
consider (2 as an integral operator in C space, invariant under rotations about 
the origin. Choosing polar coordinates C,6,q, it is easy to prove the 
following statements: 

1) The matrix of 2 in a representation in spherical harmonics Y7m(0¢): 


(7.1) f sin 049 f dpYi,A0G)QY (Op) 


is diagonal with respect to J; the diagonal elements, 2”, which are 
integral operators in one-dimensional C’ space, do not depend on m. 
2) If 2 has eigenfunctions &”) with corresponding eigenvalues w(): 


(7.2) QM) Elin) = yr) Er) 


Q has these w(”) as its eigenvalues with, as corresponding eigenfunctions, 
the products &Y: 


(7.3) QlLEM(C)¥im(O~)]=oOMEC(C)Yim(Op) (m=—l, —1+1,..., +). 


Thus the 2 eigenvalue problem in C space is reduced to the Q- 
problem in C space, which does not make things easier; more 
important, the angle-dependent part of each eigenfunction y appears 
to be a spherical harmonic. Each w(") is (2/-+ 1)-fold degenerate, to begin 
with; further degeneration is still possible, because w with different (Ir) 
may coincide. 

This is illustrated by the one case in which the 2-problem has been 
solved: that of a Maawellian gas (subscript M), in which the molecules 


repel each other according to 


(7.4) Vuln) =x/r*. 
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The solution for the Maxwell case, which was first indicated by MaxwELL 
himself, has been elaborated, using different methods, by WALDMANN ({4], 
p. 370) and Wana Cuana & UHLENBECK [7]. The general eigenfunction 
of Qy™ reads: 


(7.5) Ey) =CO*1S 1), (C*?) 


(the S(x) are Sonine polynoms, a special choice of associated Laguerre 
polynoms). We quote the lowest 3 eigenvalues, with the corresponding 
eigenfunctions, in the following table (neither } nor y=éY are 
normalized): 

TABLE 1 


ae SUES EEEEEEEEEEESEEEEERE 


Se aad a : degen- 
OM lr Ey) Yim eracy ? 
\ 00 I Yoo 
‘ 3 a 7 
wo = 0 Ol 3—C*™ Y 00 19 =D 
/ 10 ] Y11,10,1-1 
5 5 i es . 
(oz 30% 450% Yu 
om=1.374 nl 2x%/m x : y=4 
fii c+ 0+?) Ye) eee 
Te eb ‘tee : . 
(0s Soe Tom Lott Tua | 
wom = 2.062 nl 2x/m oe . ° ; 
“ 35 i yo = ) 
3 12 on( _(i*2 50%) ) a 2 © 
=5 01M / 5 5 11,10,1—1 \ 


The mo=0 entry with v9=5 is, of course, common to all potentials V(r). 
The degeneracies of the next two eigenvalues mim and moy: »1=4 and 
v2=9, probably are special characteristics of the Maxwell gas. We find, 
in particular, that the 3 components of the G, or H, vectors (see 6: 1, 2) 
are linear combinations of the 3 functions y@" ¥y», whereas the 9 com- 
ponents of the G, or H, tensors (see again 6: 1,2), of which 5 are 
independent, are all combinations of the 5 functions £2 Yo,. Thus: 


(7.6) QyG,=wimG;; QmuG,,=wemG,,. 
Owing to these properties of the Maxwell gas, many expressions given 


on the preceding pages become especially simple and elegant; e.g., 
1) the integral equations (5.42) 


G,=2mMA; G,=2mB 
are satisfied by 
Cia) A= (1/@im)G;,; B=(1/m2m)G,,; 


hence, the ChE distribution, first approximation (5.41) reads: 
° 


(7.8) $Y) = (1/oim)G,-V In 7 + (1/w2m)G,,: Veo. 
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2) The transport currents in the HD stage read (5.43): 


ge =O (He Oy" (G, Gay in T= 


W1M 


kT 5 nk?T 


nina GaG \eo= 2 7 ee 
301M BES) 2 moi VE AV’; 
(7.9a) ee - 0 
Q)=J @= D)\ eee Pica 
P Ju ae pi )) jah Mom (Ge G, c V co) a 
° ° 6 
aT es Ue 
= — / =e) een : 
5@om V co(G,, Si) 2 Wom YD ames TCO) 


apparently, the coefficients of heat conductivity and viscosity are 
expressed in terms of the eigenvalues wy by 


5 nk?T nkT 
, We 


2 mom ’ om 


(7.10) he 


In (7.9), use has been made of the following values of scalar products: 


(7.11) (G,, G,)=- 


3) Each component of both H, and H,, being an eigenfunction of 
Qy, has a y-development (6.13) consisting of one term only; the same 
is true for the developments (6.14a, b) of J and dc¢J. Consequently the 
quasi-relaxation of J, and J,, due to collisions is described by equations of 
the type (6.15): 

(7.12) dIcJ,= —omJ;; Icd,. = — wom ,,. 

So we can interpret 1/m:m and 1/m2m as the relaxation times connected 
with the transport of heat and that of momentum: 

(7.18) l/om=T%; 1/mem =T,,( = 27). 

Using these relaxation times 7 as well as the transport forces X (see 
5: 38, 39), we can write eq. (5.40) and its solution (7.8) as follows: 
(7.14) —G,-X,—G,: X,=2ud™; 

(7.15) jl) = —1,G,-X,—1,G,: “Ge 


Finally we introduce 1, and t, in the transport coefficients (7.10): 


nk? T 


ta) USN <r. 
™m ao b Le 


(3) 
(7.16) A=5 


4) As we have remarked at the end of ch. 3, the inverse of the lowest 
positive eigenvalue w; has the meaning of a relaxation time for the decay 
of ¢ during the kinetic stage. In a Maxwell gas, 

1 1 m 


(7al7)) ie ®O1M 1.374 2x reas 
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This time t ~ the mean time between successive collisions. To recognize 
this, consider a pair of molecules moving towards each other: relative 
distance—r; reduced mass governing relative motion=4m; initially 
(roo), kinetic energy of relative motion=e and velocity of this motion: 
v=(4e/m)!. They can get as near as r=o=(x/e)! (see 7.4). Apparently 
for particles with thermal velocity (for pairs or single particles) ~ 2, 
the potential field acts as if they were hard spheres of diameter o, for 
which, as is well known, the mean free path 1 ~ 1/no?. The mean time 
between successive collisions is estimated by 


ime f 
‘4 m 1,/m 


de” ont 4x’ 


(7.17a) Tor = wo) 
which is of the same order-of-magnitude as t (7.17). Incidentally, te does 
not depend on e: for a Maxwell gas, both t and te are independent of 
temperature. 

We now return to the 2 problem for the non-Maxwellian general 
intermolecular potential. Several properties of 2, of its eigenvalues w and 
-functions y have already been mentioned in ch. 2 (eqs. 2: 35-45) and at 
the beginning of the present chapter (eqs. 7: 1-38). We do not know 
and ¥ but we do have at our disposal some scalar products 


(ym), Qyylo]) = [ym ll, yy lo’), 


matrix elements of the general 2 in a representation of functions ym, 
the eigenfunctions of the Maxwellian 2 (see e.g. [2], p. 161; [3], p. 511), 
which have been calculated for use in the determination of the transport 
coefficients. The number of available elements, however, is far too small 
to make possible a diagonalization of the 2 matrix. Of course, this 
matrix has a first row and a first column (@=0; @'’=0) consisting of 
zero elements; among the first non-zero elements are the following two: 


OR) 


(7.18) 


m 


( (G,, QG,) = [G,, G,J=16 n2 22.2), 


Here 2©.°) is a temperature-dependent average “‘cross-section”’, according 
to the definition of 25) in [2], p. 157 (Q(s) in the notation of this 
reference; for a discrepancy between definitions of square bracket expres- 


sions, see (2.37)), of dimension 73/f; the matrix elements of the normalized 
functions are, of course (cf. 7.11): 


fGen eee [G., Ga] "3 3 
T19 C ee Sl = — n0(2,2) +g) == bow et _ 2 02.2) — 2 
( ) (Gi, Gi) 15 oe (ag Go he QA 


The matrix elements @, and w,, Of dimension 1/t, are entitled to take 
the places held by wim and wm in the Maxwellian gas, as is seen from 


the following remarks, which we have arranged in 4 groups, like the 
analogous discussion for the Maxwell case. 
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1) The integral equations (5.42) for the coefficients found in pl): 
(7.20) Gi= DA; G,=2B, 
are usually solved by the Enskog method of successive approximation 
to A and B. We shall discuss the first approximation only, which reads: 
(ie 1) Aw~aG,; B + 0G, 
with @ and 6 constant; with (7.20), 
(7.22) G,—2(aG,) ~ 0; G,—Q(bG,) ~ 0. 
a and 6 are now fixed by requiring that the first of the LHS in (7.22) 
shall be | G,, the second | G,: 
(7.23) (G,, G,—2[aG,])=0; (G,, G,—Q[bG,,])=0, 
or, using (7.19): 
(7.24) a=1/a,; 6=1/a,. 
So the distribution ¢ reads, in first Enskog approximation: 
(7.25) [PY i = (1/m,)G,-V In T+ (1/w,)G,,: Veo. 

2) From the formulae (5.43) for the transport currents in a Ch 
distribution 4: 

(4) =J,0 = REG , $) = 4G, AT 
(7.26) i — 
(PM =J@ = —kT(G,, 6M) = —1kT(G,, B)Veo, 

we derive the following expressions for the coefficients of transport, 
which are linear in A and B: 
(7.27) A=$k(G,, A); w= kT (G,, B). 


Using the integral equations (7.20), we find alternative expressions 


quadratic in A and B: 
(7.28) A=2k[A, A]; w=kT[B, B]. 


These formulae are still exact. Now we substitute (7.21) in order to 
determine the first Enskog approximations to 2 and y; clearly the 
requirement (7.23) ensures that (7.27) and (7.28) lead to the same result, viz. 
5 nk? T 
| [2h =4 a(G, G,)—3 


(7.29) | “e 
[whi = kT L(G, G,)= = ; 


Oy 
3) The collisional time derivative of the currents J is given by the 

approximative formula (6.12b), which we modify using the symmetry 

Obata: 

(7.30) IJ =do(H, $+4) ~ —[H, 6+4¢]=—[4+4, H1= —(4+4, OH). 
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First Enskog approximations to 2H follow from eqs. (7.22) which, of 
course, hold for H as well as for G: 


(ical) H, ~ Q(aH,); H, » 2(6H,). 


lt 


Substituting this into (7.30), we get the following equations for quasi- 
relaxation : 


\ dof, & —C/a)\fi= —oh, 


7.32 | 
( dc, FS —(1/b)J,= = 4S, 
Apparently, 

a6 15 anda a 5 = 
(7.33) Lio, =0 = TanDea —% and 1l/w,=o= aaa = t 


are the relaxation times pertinent to the quasi-relaxation of heat transport 
and of momentum transport, respectively. Formulae (7: 14, 15) giving the 
integral equation for ¢@ and its solution in terms of X and _t, and (7.16) 
expressing 2 and yw by means of these 1, and 1, all remain valid (approxi- 
mately of course) in the non-Maxwell case. 

4) To give an approximate value for the lowest positive eigenvalue 
oi, we try a variation method, using, as trial functions, three of the 
eigenfunctions ym of 2m belonging to the lowest Maxwell eigenvalue 
wim, Viz. G;: 


(7.34) 
Again, 
(7.35) 


which approximately gives the relaxation time of ¢, i.e. the time scale 
of the kinetic stage. The following argument will show that this t is once 
more ~ the average time between collisions. We apply 2 (see definition, 
eq. 2.32) to G,; over all values of c,¢1,e and x, (G,+G,'—G,, —G;,) and «, 


the differential cross-section, will have average values G and 4, and, 
looking at orders-of-magnitude : 


ap z. G7, G nG2 ] 
(7.36) QG, = naG; tT ®& {Gi Ga) yy —= = —. 
(Gj, 2G) naG2 Ne 
An estimate for & is provided by vo?, v being the mean thermal velocity 
and o? the cross-section of the molecules. So 


i l 
7.37 fe i 
( ) é ono? Vv 


(/=mean free path), in which tc is some average time between successive 
collisions. 


8. Time scale of the hydrodynamical stage 


The development in time of the basic moments y (=e, ¢0,7'), in the 
absence of external forces, is determined by the HDE (2.10), which we 


quote in their usual form: 


i 0 =—y:Vo—oV -eo 
on co = arta 


dT = — cy VT 7g (P: Veo+V-9). 


— 
oe) 
— 

ces 


They follow from the BE by integration and are therefore valid during 
both the kinetic and the HD stage. During the kinetic stage, the rapid 
changes which the distribution function f undergoes will also be reflected 
in the ¢-derivatives of its moments y; when / has reached its normal 
ChE value f©2 and the HD stage has set in, further changes of the 
moments y will be much slower, as we shall now demonstrate. 

The ChE distribution /“ has the form of a “‘Hilbert expansion” (2.22): 


(8.2) fon FOL 4 60 4 $+. 


with successive terms of decreasing size; in fact, introducing a small 
positive number gy, we may suppose that 


(8.3) YD) wo; d@) w op; ... (O<@ < 1) 


If we restrict ourselves to the first correction, 6“), the transport quantities 
P and q in (8.1) are given by the phenomenological equations: 


(8.4) P—pU—2uiFeq and q=—A/T, 
with 

(8.5) p=koT/m 

and 

(8.6) Aas. oy. 


We substitute (8: 4, 5, 6) into (8.1), with the result: 


ot 0 =— ° 
o¢ Co = 2 (u/o)V - (Veo) 
ip 
= (uJ) \Ac+ 3 V+ (Veo) 
a +4mp (3ko)+ (Vco)2 


(3a) (3b) 


(The grouping of the RHS terms will be discussed in connection with 
eqs. (8.13)). 

The magnitude of the gradients of 7’ and ¢o occurring in (8. Hy) ae 
restricted by the estimates (8.3). To begin with, the estimate for 4 
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determines the magnitude of the first-order gradients; for, using (7.15), 
we have 


te} 


° — 
(8.8) pW) 7,(3—C*2)C-V nT — Ep Veo = 6,04 fv) A P. 


t,—the kinetic relaxation time Tt; t,=j7;3—C*? and C*C™ are of order 
unity (in the velocity range of interest) and C is comparable to the mean 
thermal velocity v=(kT'/m)!. Hence, 


(8.9a) d,0 wx w(VT/T) y= ~ <1; 
(8.9b) o,Y) w~ Wo ry <i. 


tu=lis a length ~ the mean free path; so (8.9a) illustrates the fact that, 
in a development of f/f in powers of lV, 4 is the linear term, which 
is small if the relative change over a distance / of the moments y (e.g., 7’) 
is <1. Likewise, (8.9b) restricts the change along / of ¢) to an amount 
small compared to //r=v. Apparently the thermal velocity v serves as 
the proper velocity of reference for co. 

Although ¢@ does not contain a g-term comparable to the V7'- and 

° 


’co-terms, it stands to reason that a condition similar to (8.9a, b) is 
imposed on |/o: 


(8.9¢) (Volo) > @ <1. 


For the gradient of the hydrostatic pressure p = ov?, (8.9a, ¢) imply a 
restriction 


(8.10) (l/ov?) Vp =p <1. 


In a similar way, we use the estimate of ¢@) in (8.3) in order to determine 
the magnitude of the second-order gradients of 7’ and ¢ 9 occurring in 
(8.7). This second correction consists of terms reading 2VV ; VV symbolizes 


products of first gradients (2 : +) as well as second-order gradients 


(like Ay/p). By virtue of (8.3), all these terms will be ~ ?; e.g., 


(8. 11a) b,a® w (P/T) AT » @; 


(8.11b) bi aw (P/v) V (Veo) © @. 


Finally, the magnitude of the viscosity coefficient j~ can be inferred 
from the connection (7.35) between Q@2) and the kinetic relaxation 
time Tt: 

SkT mv 


Sele L FS fi — ae ~S 
( ) wwe [uh 30 


a) oe 
2,2 


w ovl. 


We are now in a position to replace each term in the RHS of the 
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HDE (8.7) by its order-of-magnitude; using the estimates (8: 9, 10, 11, 12) 
we get: | 


All terms inside the brackets at the RHS are <1, as we may safely 
assume Co not to exceed v. Each term corresponds to a contribution 
Td:y/y, and therefore they indicate the local changes in 9, ¢o, 7’ relative 
to 9, v, T in a time t which it has previously taken the gas to pass through 
the kinetic stage. Apparently these changes are all <1, so that the HD 
development is certainly slow in comparison to the kinetic process. 

The grouping of terms in eqs. (8.7) and (8.13) has been made in order 
to distinguish between HD processes of a different nature; for each 
process we quote a characteristic time, viz. the time after which, 
according to (8.13), the process would have caused a relative change in wp 
or order unity: 

1) local changes in all three y due to bulk movement with velocity co: 


(8. 14a) tor,1 = tv Cop. 


2) remaining first order effects (decrease of o and 7’ if co is locally 
divergent; acceleration of the gas down a hydrostatic pressure gradient) : 


(8.14b) tir,2 — t/@. 


3a) second degree, first order process (dissipation of kinetic energy 
of bulk movement by friction), and 

3b) second order processes (local yield of transport of momentum and 
heat); for both (a) and (b), 


(8.14c) tr1,3 = t/q. 


Estimates of the y-gradients more specific than those in terms of the 
development parameter g are provided by the following suppositions: 

a) the local gradient /y, if continued along the vessel dimension d, 
would correspond to a relative p difference between the extremities of 
a ot order 1: 


(8.15a) UV p/p) ~ Id; 
b) in analogy, the gradient of Vy is restricted by 


(8.15b) (I/p) VV) ~ Vea. 
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Using these estimates instead of (8: 9, 10, 11), we obtain the following 
HD times for the various processes mentioned: 


(8.16) ir. d/co; tir,2= dv; in,3 = d2/lv. 


All these times are S> rt. Which of the corresponding processes predomi- 
nates, depends on which is the shortest. In general co<v, so the first 
order processes (2) will in general prevail. In their absence, streaming (1) 
and 2nd order or 2nd degree processes (3) can act. Which of these is first, 


depends on the ratio: 
(8.17) a 


This appears to be Reynolds’s number, and we can distinguish between 


the following 3 cases: 

a) Re <1, laminar flow: t; > tg: processes of 2nd order and 2nd 
degree are preponderant. 

b) Re x 1, flow still laminar: t; ~ ts: both groups of terms contribute 
to the HD development. 

c) Re>+1000, turbulent flow: first order terms dominate. 
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List of symbols 


(Missing symbols are explained near to the place where they are used. 
Superscripts and subscripts are mentioned separately at the bottom. In 
the equations quoted, the symbol is defined, first mentioned, or best 


elucidated.) 


= A) 
5.41 ) 
c 
(2.7) Co 
Cc 
(2.45) Gs 
(2: 33, 34; 3: 18, 14; lp 
4217, 18) ) 
i 
(Ge) seed 
(5.3) H 
(23) J (ffi) 
(5: 38, 39) tba 
l 
m 
n 
: Ve 
A 
(5.43) iq 
(5.5) 8 
4 
Cis) tt,1,101 
Dv 
(5: 38, 39) Kui Ay 
(2.00; a1; 415) 6 
= (A 
(7.9) 
(2.7) 0 
(ied) T 
(2.22; 3.8; 4.3) d 
(8.3; 3.39) Pp 
(2,392.5) yl] 
(2.6) yl 
(2.32) Q 
(2.39) o 
(3.25, etc.) o 
ot 
(22d 0.7) ds 
(273 3.5,138%6.7) dc 
(6.7) oH 
(2: 22, 23, 26, etc.) (r), (s) 
Ch 
(3.3; 4.3) ~ 


coefficients in 6 


velocity of a molecule 
local average velocity 
C— Co 

reduced value of C 


LHS of integral equations 


velocity-space distribution 

functions of C 

H-function (used in ch. 5 only) 

collision operator 

transport currents 

Boltzmann's constant 

mean free path 

mass of a molecule 

local number density 

pressure tensor 

heat current 

local entropy density 

temperature 

characteristic times of stages I, II, II 

mean thermal velocity (kT'/m)? 

transport “‘forces”’ 

series development parameter (= 1) 

coefficient of viscosity 

coefficient of heat conductivity 

mass density 

kinetic relaxation time 

deviation of f relative to f® or feb 

series development parameter 

eigenfunction of Q 

moment of 7/01 

linearized collision operator 

eigenvalue of 2 

sum of m values 

partial derivative with respect to ¢ 

contribution to oy by “streaming” 

contribution to dy by collisions 

contribution to % by change of H 

e.g.,/, pu, d,6), PO: term index 
in series of Hilbert type 

e.g., foo: ChE value 

ee. i,¢é, F: deviation from MB 
value (ch. 3) or from ChE value 
(ch. 4) 


2: 39, 46, etc.) 
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[wa | e.g 
[@ | © 
A, bl e 
ei. 
( ) e 
[ ] e 


, 741, pl: indices of 2-eigen- 
aon y and of y-development 
coefficients m; A counts degenerate 


levels 


.o., Wil; index of development 
coefficients in e-“!-series 
.2., Tj, T,: quantities connected with 
transport of heat/momentum 
=, : 
2 ail ! T px) 
ae Ls 
( Tr = $2 ro 
.2., (¢, 7): scalar re 
.2., [A, i matrix element of 2 


ORGANIC CHEMISTRY 


DISSOCIATION CONSTANTS OF 
ALKYL-SUBSTITUTED CYCLOHEXANECARBOXYLIC ACIDS 
AND CYCLOHEXANEACETIC ACIDS 
(preliminary communication) 

Bie 


H. VAN BEKKUM 1), P. E. VERKADE, anp B. M. WEPSTER 


(Communicated at the meeting of November 26, 1960) 


During the realization of a research program concerning synthesis, 
configuration, and conformation of cyclohexane derivatives a fairly large 
number of pure alkyl-substituted cyclohexanecarboxylic acids and 
cyclohexaneacetic acids became available. The apparent dissociation 
constants of these acids in 50 vol.°% ethanol-water at 25° were determined 
potentiometrically. Considering the present-day interest in cyclohexane 
derivatives, and in particular also in the acids in question, it would seem 
desirable to us to record these data in a preliminary communication. For 
additional information reference is made to the thesis of one of us (v. B.) 
shortly to appear and to papers to be published elsewhere at a later date. 

The cis- or trans-configuration of nearly all the 4-alkyleyclohexane- 
carboxylic acids involved in the present work has been established by us, 
as was that of the 2- and 3-methylcyclohexanecarboxylic acids by others, 
by a chemical method, viz. by correlating these compounds with the 
corresponding cyclohexanedicarboxylic acids, the configuration of which 
had been put beyond doubt. The assignment of the configuration to the 
other alkyl-substituted cyclohexanecarboxylic acids took place on the 
basis of the method of preparation and the stability of these compounds, 
while making use, wherever necessary and possible, of data in the literature 
of other nature. Finally the alkyl-substituted cyclohexaneacetic acids were 
correlated chemically with the corresponding alkyl-substituted cyclohexane- 
carboxylic acids. We feel sure that the configuration of the acids thus 
deduced by us is the correct one. 

For the determination of the apparent dissociation constants we used 
a Beckman Model G pH-meter with a glass electrode and a calomel 
reference electrode. Each compound was examined in three separate 
solutions at 1/3, 1/2, and 2/3 neutralization. The mean of the three values 
for the pK, thus found, which were invariably in excellent agreement, 
is given in the Tables below. 


1) Present address: Koninklijke/Shell-Laboratorium, Badhuisweg, Amsterdam-N. 
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In addition the melting points of the acids in question, as found by us, 
are listed in the Tables. Finally —invariably assuming the cyclohexane ring 
to have the rigid chair form, which is undoubtedly reasonable —for all 
the compounds examined the predominating conformation of the sub- 
stituents has been indicated with the aid of the letters e (for equatorial) 


Alkyl-substituted cyclohexanecarboxylic acids 
oe OE Eee 


gilevl opewote) | cis-isomer “) trans-isomer *) 

ie ela m.p. conf. pKa m.p. conf. pKa 
H | |  30-31° € 6.30 
4-methyl 30.5-31.5° ae 6.53 110.5-111.5° — ee 6.26 
4-ethyl |  22.5-23.5° ae 6.57 |  49.5-50° ee 6.28 
4-isopropyl | 40—41° ae 6.66 94-94.5° ee 6.26 
4-tert-butyl «118-1185 ae 6.78 | 175-176" ee 6.28 
4-neopenty! 64-64.5° ae 6.57 120.5—-121° ee 6.28 
3-methyl] b) ee 6.31 
3-isopropyl i) ee 6.41 
3-tert. butyl 94.5-95.5° ee 6.42 119-119.5° ae 6.95 
2-methyl 20—21° ae 6.52 53.5-54.5° ee 6.24 
2-ethyl 37.5-38° ae 6.67 75—76° ee 6.38 
2-isopropy] 35.5-36.5° ae 7.16 46.5—47.5° ee 6.47 
2-tert. butyl 64-65° ae 7.45 58.5—59.5° ee 6.65 
2,6-dimethyl | 87—88° aee 7.21 103—104° eee 6.34 
2,4,6-trimethyl 105.5-106° aeee —- 7.30 96—96.5° eeee =: 688 
3,5-dimethyl | 69-69.5° eee 6.35 
3,5-di-tert. butyl 192.5-193.5° eee 6.56 
1-methyl 37.5-38.5° ae 6.79 


“) In the acids with two or three alkyl groups examined by us these groups 
all occupy the same position relative to the ring; the acids in question are 
meso-acids. Cis and trans here indicate the position of the carboxyl group 
relative to the alkyl groups. 

6) Liquid; n7> 1.4548. 

¢) Liquid; ny 1.4648. 


Alkyl-substituted cyclohexaneacetic acids 
re ee en ee 


alkyl group(s) cis-isomer trans-isomer 

| m.p. cont.» “pity | m.p. conf. pK, 
H | 30-31° e 6.18 
4-methyl | 43-44° ea 6.24 73-73.5° ee 6.20 
4-isopropyl 40—41° ae 6.26 78-78.5° * ee 6.21 
4-tert. butyl | 98-99° ae 6.32 96.5-97° ee 6.23 
3-tert. butyl | 34-35° ee 6.27 
2-tert. butyl | 111-112° ae 6.44 
3,5-di-tert. butyl 149—150° ece 6.38 
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and a (for axial); the letter coming first refers to the conformation of 
the carboxyl group. 


In the following paragraphs we shall confine ourselves to a few brief 
remarks with reference to the pKa values listed in the Tables. For a 
fuller discussion of these numerical data reference may be made to the 
thesis of one of us (v. B.) shortly to appear and to papers to be published 
elsewhere at a later date. 

(1) Applying the basic idea of the work of WrysrEern and Honnsss 2) 
to certain compounds examined by us, we assume that the dissociation 
constants found for cis- and trans-4-tert.butyleyclohexanecarboxylic acids 
are those of the conformationally pure compounds (ae and ee respectively). 
On the basis of these two dissociation constants and of that found for 
cyclohexanecarboxylic acid, it is then possible to calculate, in the way 
introduced for reaction rates by the investigators just mentioned and 
applied by Srontow 3) and by Ticuy, JonAs, and Srcuer *) to dissociation 
constants, the free-energy difference between the two conformations of 
cyclohexanecarboxylic acid and of cyclohexane carboxylate ion, with the 
substituent in the axial and the equatorial position respectively; this 
calculation yields the following results: 44 Fcoon=1.6 kcal/mole and 
AAFcoo- =2.3 kcal/mole. 

(2) Considering the very small difference between the dissociation 
constants of trans-4-tert.butylcyclohexanecarboxylic acid and cyclohexane- 
carboxylic acid, it is no doubt more advisable to base the calculation of 
the two above-mentioned as well as any further 44/-values on the 
generally accepted free-energy difference between the two conformations 
of methyleyclohexane, 44 F'cu,= 1.8 keal/mole ®). We then follow the 
procedure applied by Sronow 3) and by Ticny, JonAs and Srcumr 4). 
On the basis of the dissociation constants of cis- and trans-4-tert.butyl- 
cyclohexanecarboxylic acids, of those of the other czs-4-alkyleyclohexane- 
carboxylic acids examined by us, and of the said value of A4Fcx,, the 
following free-energy differences have been calculated: 


AAFcoou = 1.5 kcal/mole 


AALReao- = 2:1 
AAF xt we 70) 
AAFipr = 2.5 


AA F Neop — 748) 
(3) The cis-3-alkyl- and cis-3,5-dialkyleyclohexanecarboxylic acids 
examined by us, with practically complete equatorial position of all the 


S. Wrinsrern and N. J. Houness, J. Am. Chem. Soc. 77, 5562 (1955). 

R. D. Srotow, J. Am. Chem. Soc. 81, 5806 (1959). 

M. Ticxy, J. JonAS, and J. Stcumr, Collection Czechosl. Chem. Commun. 24, 

3434 (1959). 
5) C. W. Becxert, K. 8. PirzEr, and R. Sprirzer, J. Am. Chem. Soc. 69, 2488 
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substituents, have a higher pKa value than cyclohexanecarboxylic acid; 
thus in these cases a small but marked acidity-weakening effect of the 
alkyl group(s) is unexpectedly found to occur, 

(4) In the 2-alkyl- and 2,6-dialkyl-substituted cyclohexanecarboxylic 
acids steric hindrance of the carboxyl group and the carboxylate ion by 
the alkyl group(s) can no doubt play an important role. In the cis-isomers 
in question, with the carboxyl group largely in axial position, in addition 
conformational hindrance by the axial hydrogen atoms in positions 
3 and 5 will occur. It is no doubt partly due to this that in these cases 
the difference between the dissociation constant of a cis-isomer and that 
of the corresponding trans-isomer increases considerably with increasing 
size of the alkyl group or with increasing number of these groups. The 
cis-acids of the type under consideration include remarkably weak acids. 

(5) In the alkyl-substituted cyclohexaneacetic acids axial or equatorial 
position of the —CH2.COOH group is found to give rise to only small 
differences of acidity. This was to be expected. 


For some of the acids examined by us, the dissociation constant had 
already been determined by other workers, though in other media #)*)®)7)8), 


Laboratory of Organic Chemistry, 
Technical University, 
Delft, Netherlands 


6) N. Zevrysky and N. IzGaryscuew, J. Russ. Phys. Chem. Soe. 40, 1379 (1908); 
Chem. Zentralbl. I, 531 (1909). 

*) J.F.J. Dippy, 8. R. C. Hugues and J. W. Laxton, J. Chem. Soc. 4102 (1954). 

8) P. F. Sommer, V. P. Arya, and W. Sron, Tetrahedron Letters 20, 18 (1960). 
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SYNTHESIS OF UNSATURATED HYDROCARBONS BY THE 
WITTIG REACTION 
(PRELIMINARY COMMUNICATION) 


BY 


H. HEITMAN, J. H. SPERNA WEILAND anv H. 0. HUISMAN 


(Communicated by Prof. J. F. ARENs at the meeting of September 24, 1960) 


A recent publication by McDonatp and Campsety [1] in which the 
attempted preparation of 1,6-diphenyl-1, 3,5-hexatriene (VI) was described 
by means of a Wirtic reaction with 1,4-bis(triphenylphosphonium)- 
2-butene dichloride and benzaldehyde, prompt us to communicate our 
own results in the same field at an early stage. 

1,4-bis(triphenylphosphonium)-2-butene dibromide (V) was prepared 
by selective partial hydrogenation of 2-butyn—1,4-diol (I) to 2-butene— 
1,4-diol (II) followed by treatment with dry hydrogen bromide and 
conversion of the so formed 1,4-dibromo—2-butene (III) into the corre- 
sponding bistriphenylphosphonium compound (V) by a modification of 
the method used by FRrepRIcH and HEennin@ [2] for the preparation of 
1,3-bis(triphenylphosphonium)propane dibromide using dimethylforma- 
mide as the solvent. 


; Lindlar HBr seh WC 
HOCH.2C = CCH20H ————_~+ HOCH2CH = CHCH2OH — rae BrCH2CH = CHCH2Br 


He iv 
if 1B III trans 
(m.p. 55.5—56°) 
DMF Br- Br- 
BrCHsCH = CHCH2Br + (CeHs)3 P —————> (Ce6Hs)3 P+ — CH2CH = CHCH2 — P+(CeHs)s 
EE Vi Vi 


Contrary to the results of McDonaLp and CampBe.n [1] the compound 
VI could be obtained in good yield by a Wrrric reaction carried out in 


alcoholic solution and with a lithiumalkoxide as the base. 
Bre Br= nase Sw fo 
(CeHs)3 P+ —CH2CH = CHCH2 — P+(C6Hs)s + 2 CeHsCHO > (0) 
Vv VI 


(m.p. 200.5-201.5°) 
max: 335, 351 and 370 my (ethanol) 
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The reaction was carricd out also with anisic aldehyde (VII), cinnamic 
aldehyde (VIII), f-ionylidene acetaldehyde (IX) and furfural (X), all 
leading to the corresponding unsaturated hydrocarbons (XI-XIV). 


S OCHs 
_ base Vw “& 
2 CHs0 —(O )— CHO + ¥ (0) 
CHs30 
VII XI 
(m.p. 250-251°) 
Amax: 346, 362 and 380 my (ethanol) 


_ base Yee we “DS e 
2CsHs;sCH = CHCHO + V = 


Vill XII 
(m.p. 252—254°) 
Amax: 380, 402 and 426 my (benzene) 


H. 
ye eg a base Be a ea ee 
2, LV 
IX 


XIIT 
(m.p. 160-161°) 
Amax: 412 my (cyclohexane) 


x XIV 
(m.p. 154-155°) 
Amax: 348, 366 and 388 my (cyclohexane) 


When the reaction was carried out with ketones we were unable to 
obtain the hydrocarbons to any substantial amount, although there 
were indications that traces of the desired compounds were formed. 

By taking a trialkylphosphine instead of triphenylphosphine (IV) the 
reaction took place much faster and the yields were considerably better. 

More detailed studies about the influence of the substituent in the 
aromatic nucleus of the aldehyde on the yield of the reaction are in progress. 

There are some indications leading to the assumption that in the 
cases described above one of the steps of the Wrrria reaction might 
proceed by a free radical mechanism. 

Full details of this work will be published in the thesis of Mr H. Herrman 


as well as in forthcoming TS. 
ac A Laboratory for Organic Chemistry 


of the University of Amsterdam 
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